> | 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); |
> | with(DEtools): |
> | m:=9e-31; q:=1.6e-19; B:=0.01; |
> | 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]); |
> |