figure() subplot(2,1,1) t = linspace(0,3,51) y1 = f1(t) y2 = f2(t) plot(t,y1,'b-',t,y2,'bo') hold = 'on' xlabel('t') ylabel('y') axis([t[0],t[-1],min(y2)-0.05,max(y2)+0.5]) legend(['t^2 * exp(-t^2)','t^4*exp(-t^2)']) title('Top Figure') t3 = t[::4] y3 = f2(t3) plot(t,y1,'b-',t3,y3,'ys') xlabel('t') ylabel('y') axis([0,4,-0.2,0.6]) legend(['t^2 *exp(-t^2)','t^4*exp(-t^2)']) savefig('tmp4.pdf') show()