R = -0.1:0.01:0.05; for i=1:length(R) S=expGrowth(1000, R(i) ,30); FINAL(i,:) = S ; end surf(1:30 , R , FINAL)