function E = electric(t,y) % ELECTRIC Calculates electric field E(t,y) % Input values are scalar time t and vector y(6)= (x,y,z,vx,vy,vz) % Output values are column vector E(3)=(Ex,Ey,Ez) % E=zeros(3,1); % Electric and Magnetic Fields E(1)=0. ; E(2)=0.; E(3)=0. ;