Trayectoria del electrón

n = 92;
Load('EHField_3D_5to15fs.mat');
XX (1) = 0;
YY(1)=y(40);
ZZ (1) = 332 * 10 ^ -9;
VX(1) = 0;
VY(1) = 0;
VZ(1) = 10 ^ 6;

para i = 2: n/c

velocidad = vt(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),i-1);
VX(i) = speed(1,1);
VY(i) = speed(1,2);
VZ(i) = speed(1,3);
posición = xz(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),i-1);
XX(i) = position(1,1);
YY(i) = position(1,2);
ZZ(i) = position(1,3);

plot3(XX,YY,ZZ)
Espera
final

Y ahora, tengo que simular desde muchos tiempo inicial diferente. Había escrito este código

borrar todo;
n = 92;
Load('EHField_3D_5to15fs.mat');
XX (1) = 0;
YY(1)=y(40);
ZZ (1) = 332 * 10 ^ -9;
VX(1) = 0;
VY(1) = 0;
VZ(1) = 10 ^ 6;

para a = 90
para j = a: n
para i = 2: n/c-a + 2

velocidad = vt(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),j-1);
VX(i) = speed(1,1);
VY(i) = speed(1,2);
VZ(i) = speed(1,3);
posición = xz(xx(i-1),yy(i-1),zz(i-1),vx(i-1),vy(i-1),vz(i-1),j-1);
XX(i) = position(1,1);
YY(i) = position(1,2);
ZZ(i) = position(1,3);
plot3(XX,YY,ZZ)
Espera
final
final
final