How to store a vector under a for loop

Hi, so I want to store theta as a vector under my for loop. But only the last number of theta saves, does anyone know how to get that as a vector?

clear all;

v = 5.8;
u = 1/v;
Psph = [509.35 507.06 504.77 398.20 394.19 389.61 268.71 266.99 265.85];
E = [];
r1= 6731;
r2= 6731;
pi = 3.1416;
theta=[];

for delta = [29 30 31 59 60 61 89 90 91]

      theta = r1*sin(delta)*u;

end

   E = Psph/((r1.^2)*(r2.^2)*4*pi*(u.^2)*sin(delta)*(cos(theta)).^2);
      A = sqrt(E);

Answer 1

[James Tursa]

theta(delta) = r1*sin(delta)*u; % <-- Added the (delta) indexing

You will also need to change some of your matrix operators to element-wise operators. E.g.,

E = Psph./((r1.^2)*(r2.^2)*4*pi*(u.^2)*sin(delta).*(cos(theta)).^2);

But note that this entire loop:

for delta = [29 30 31 59 60 61 89 90 91]
      theta(delta) = r1*sin(delta)*u;
end

can be replaced with these lines:

      delta = [29 30 31 59 60 61 89 90 91]
      theta = r1*sin(delta)*u;

[Maria Gonzalez]

Thank you very much!

Answer 2

[Andrei Bobrov]

v = 5.8;
u = 1/v;
Psph = [509.35 507.06 504.77 398.20 394.19 389.61 268.71 266.99 265.85];
r1= 6731;
r2= 6731;
pi = 3.1416;
delta = [29 30 31 59 60 61 89 90 91];
theta = r1*sin(delta)*u;
E = Psph./((r1.^2)*(r2.^2)*4*pi*(u.^2)*sin(delta).*(cos(theta)).^2);
A = sqrt(E);

[Maria Gonzalez]

Thanks!