Witam, mam problem z pętlą for.Chciałbym napisać pętle liczącą taką tabelę

screenshot-20200402190028.png

Nie wiem jak to zapisać aby w iteracji zerowej wartość ekw = 0.01, a w następnych iteracjach liczyło ją z iteracji poprzedniej.

a = 6378137;
GM = 3986005*10^8;
J2 = 108263*10^(-8);
omega = 7292115*10^(-11);

ekw = 0.01;
e = sqrt(ekw);
ep = e/sqrt(1-ekw);
z = (1+(3/(ep^2)))*atan(ep)-(3/ep);
ekw(1) = (3*J2)+((4*(omega^2)*(a^3)*(e^3))/(15*GM*z));

for i = 1:6
     el(i) = sqrt(ekw(i));
     epl(i) = el(i)/sqrt(1-ekw(i));
     zl(i) = (1+(3/epl(i)^2))*atan(epl(i))-(3*epl(i));
     ekwl(i+1) = (3*J2)+((4*(omega^2)*(a^3)*(e^3))/(15*GM*zl(i));
 end
``