from pylab import *

Xmin = 0.0;  Xmax = 5.0;  Npoints = 500
DelX = (Xmax - Xmin)/Npoints
x = arange(Xmin, Xmax, DelX)
y = sin(30*x) + sin(33*x)

plot(x,y)

show()