Originally Posted by

**Rorshach** Hey,

I have been trying( for the last several hours...) to make this horrid thing work(**you don't have to read the whole script, the main point is the several loops and the plot**):

I(1)=10^2;

II(1)=10^2;

III(1)=10^2;

IV(1)=10^2;

c=4;

s=0.5;

for k=1:2:9

T=0.1*k;

for j=1:2:9

f=0.1*j;

for h=1:5

M=0.5+0.5*h;

N(1)=4*10^2;

for i=1:20

I(i+1)=s*c*(2/(N(i)-1))*(I(i)*((1-T)*((N(i)-1)/2*c)+I(i)-1+(f/2)*II(i)+III(i)/2+(f/4)*IV(i)+M*T*(2*(I(i)-1)+f*II(i)+III(i)/2+f/4*IV(i)))+(f/4)*II(i)*III(i)*(1+M*T));

II(i+1)=s*c*(2/(N(i)-1))*(II(i)*((1-T)*((N(i)-1)/2*c)+(f/2)*I(i)+(f^2)*(II(i)-1)+(f/4)*III(i)+(f^2/2)*IV(i)+M*T*(2*(f^2)*(II(i)-1)+f*I(i)+(f/4)*III(i)+f^2/2*IV(i)))+f/4*I(i)*IV(i)*(1+M*T));

III(i+1)=s*c*(2/(N(i)-1))*(III(i)*(((N(i)-1)/2*c)+(1/2)*I(i)+(f/4)*II(i)+III(i)-1+(f/2)*IV(i)+M*T*((1/2)*I(i)+(f/4)*II(i)))+(f/4)*I(i)*IV(i)*(1+M*T));

IV(i+1)=s*c*(2/(N(i)-1))*(II(i)*(((N(i)-1)/2*c)+(f/4)*I(i)+(f^2/2)*(II(i))+(f/2)*III(i)+(f^2)*(IV(i)-1)+M*T*((f^2/2)*(II(i))+(f/4)*I(i)))+f/4*II(i)*III(i)*(1+M*T));

N(i+1)=I(i+1)+II(i+1)+III(i+1)+IV(i+1);

end

x=(1:i+1);

plot(x,I,x,II,x,III,x,IV);

end

end

end

The thing is that I am trying to some how save all of the plots so I can view them later. But since they are inside a loop,and there's a whole lot of them, simple thing like 'subplot' will not do.

Anyone knows how do I solve this issue?

Thanks in advanced.