&input nt=1000, dt=1, neqn=6, abserr=1.e-5, relerr=1.e-5, eq=1.602, em=9.1, oq=1., om=1., beta=0.1, e=0.,0.,0.001, b=0.,0.,0.1, y=1.,0.,0.,0.,1.e-3,0. /