;$Id: pmat.pro,v 1.4 2016/10/26 22:06:37 brandenb Exp $
;
;  The eigenfunctions are plotted afterwards using ".r pevec.pro"
;  But before that, run on the command line ".r rayleighII.pro".
;  Report problems to Axel Brandenburg <Axel.Brandenburg@Colorado.edu>
;
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=8,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv -i idl.ps ~/tex/teach/ASTR_5410/Homeworks/fig/pmat.ps
;
!p.charsize=1.0
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
!p.multi=[0,2,1]
!x.title='!6'
!y.title='!6'
;
nscl=20
loadct,6
xxx=findgen(nscl*21)/nscl
contour,congrid(m,nscl*21,nscl*21),/fill,nlev=40,xxx,xxx
contour,congrid(s,nscl*21,nscl*21),/fill,nlev=40,xxx,xxx
loadct,5
;
fo="(11f7.3)"
print,m(0:10,0:10),fo=fo
print
print,s(0:10,0:10),fo=fo
END