sum=0; step=.1; for i=1:10 x(i)=i/10.; sum=sum+x(i)*step; y(i)=sum; end plot(x,x,'b'); hold on plot(x,y,'r'); hold on plot(x,y,'rx'); hold on plot(x,(.1)./y,'g'); y %print -dpng myplot.png; Copyright (C) 2009 Konstantin Kirillov