elektron.mws

>    deq1:=m*diff(x(t),t,t)=q*B/(1e16*x(t)^2+1)*diff(y(t),t); deq2:=m*diff(y(t),t,t)=-q*B/(1e16*x(t)^2+1)*diff(x(t),t);

deq1 := m*diff(x(t),`$`(t,2)) = q*B/(.1e17*x(t)^2+1)*diff(y(t),t)

deq2 := m*diff(y(t),`$`(t,2)) = -q*B/(.1e17*x(t)^2+1)*diff(x(t),t)

>    with(DEtools):

>    m:=9e-31; q:=1.6e-19; B:=0.01;

m := .9e-30

q := .16e-18

B := .1e-1

>    DEplot3d([deq1,deq2],[x(t),y(t)],t=0..10e-9,
[[y(0)=0,x(0)=0,D(y)(0)=1,D(x)(0)=-20]],scene=[t,x(t),y(t)],linecolor=t,stepsize=5e-11,orientation=[0,90]);

[Maple Plot]

>