pcheck_sin4.py

import matplotlib.pyplot as plt
import numpy as np
#plt.ion only works in interactive mode
plt.ion()
plt.rc('text',usetex=True)

#plt.rc('font',**{'family':'sans-serif','sans-serif':['Helvetica']})
## for Palatino and other serif fonts use:
plt.rc('font',**{'family':'serif','serif':['Times New Roman']})



x=np.logspace(-4,3*np.pi,100)
y=np.abs(np.sin(x))*np.exp(-x)
x2=np.logspace(-4,3*np.pi,40)
y2=np.abs(np.sin(x2))*np.exp(-x2)
y3=np.abs(np.sin(x+.4))*np.exp(-x)

plt.subplot(211)
plt.title('Height of ocean swells'r'~\Gamma',fontsize=24,color='b')
plt.yscale('log')
al=plt.plot(x,y,'k--',label='theory',linewidth=2.)
bl=plt.errorbar(x2,y2,xerr=.0005,yerr=.02,fmt='r^',label='data')
plt.xscale('log')
plt.ylabel(r'\Gamma',fontsize=24)

#plt.setp(al,linewidth=4.)
plt.legend(loc=2)
plt.xlim([0,10])
plt.ylim([1.e-3,1.3])

plt.subplot(212)
al=plt.plot(x,y3,'k--',label='bad theory',linewidth=2.)
#plt.setp(al,linewidth=4.)

bl=plt.errorbar(x2,y2,xerr=.0005,yerr=.02,fmt='r^',markersize=8,label='data')

plt.legend(loc=2)
plt.yscale('log')
plt.xscale('log')

plt.xlim([0,10])

#plt.ylim([0,1.3])
plt.ylim([1.e-3,1.3])

plt.ylabel(r'\Gamma',fontsize=24)
plt.xlabel(r'\omega t',fontsize=16)
plt.show()
plt.savefig('ocean.pdf',format='pdf',dpi=1000)