##^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^Lance Simms, Stanford University 2009
from pylab import *
import pyfits

a=[1,2]
b=[1,2]

plot(a,b)
xlabel(r'$\lambda$','hello',size=12)
ylabel('hello in units of'+r'$\lambda (\angstrom)$',size=14)
show()