a=dindgen(101)/5.
c=1 
;window,xs=800,ys=800 
set_plot,'ps'
device,file='rpsf.ps',/landscape
plot,a,(exp(-a^2/(2*c^2)^2)/(c*sqrt(2*!pi))),/ylog,yr=[1e-8,1],chars=1.5,xtit='Radius (arcsec)',ytit='Radial Point Spread Function'
c=1.1  
oplot,a,(exp(-a^2/(2*c^2)^2)/(c*sqrt(2*!pi))),linestyle=2
c=1.3
oplot,a,(exp(-a^2/(2*c^2)^2)/(c*sqrt(2*!pi))),linestyle=3
device,/close
set_plot,'x'
end