from scitools.all import * # import numpy and plotting t = linspace(0, 3, 51) # 51 points between 0 and 3 y = t**2*exp(-t**2) # vectorized expression plot(t, y) hardcopy('tmp1.eps') # make PostScript image for reports hardcopy('tmp1.png') # make PNG image for web pages