function ydot=emode(t,yn,E,B)
%global Bx,By,Bz,Ex,Ey,Ez
global m q
ydot=zeros(6,1);
    
%yn=[x;y;z;vx;vy;vz];
B=Bfield(yn(1),yn(2),yn(3));
E=Efield(yn(1),yn(2),yn(3));
vx=yn(4);
vy=yn(5);
vz=yn(6);

ax=q/m*((vy*B(3)-vz*B(2))+E(1));
ay=q/m*((vz*B(1)-vx*B(3))+E(2));
az=q/m*((vx*B(2)-vy*B(1))+E(3));
f=[vx;vy;vz;ax;ay;az];

ydot=f;