Thanks Madhan Ravi for the idea,
You could apppend the values of Sbar1 at each iteration to another variable, say result like this
Nn=4;
result = double.empty(0,3);
for i=1:1:Nn
for X1=0:45:360
c= cos(X1*pi/180);
s= sin(X1*pi/180);
R= [1 0 0; 0 1 0; 0 0 2];
A1= [c^2 s^2 2*s*c; s^2 c^2 -2*s*c; -s*c s*c c^2-s^2];
Sbar1= [R*A1^-1*R^-1*s*A1];
result = [result,Sbar1]
end
end
Here first 3 columns of result will be the Sbar1 of first iteration, next 3 columns will be Sbar1 of second iteration, and so on.
