cadekd

I hope is not too late for an answer here... You have two options for plot what you want: 1. Add an index to the "I" in to the loop: i=1; for T = 3000:1000:8000 [i(:,i)] = intensity(L,T,c,h,k); i=i+1; end 2. just add the plot inside the loop but you need to include the "hold on" instruction to keep the previous plot. for T = 3000:1000:8000 = intensity(L,T,c,h,k); plot(L,) hold on end grid on