from __future__ import division
from scipy import *              # Load SciPy
import pylab                     # Load Pylab

x = linspace(-pi, pi, 200)       # Generate a list of 200 points from -pi to pi
y = sin(x) + 2*cos(3*x)          # Compute a function of x

pylab.plot(x, y)                 # Plot y against x
pylab.show()