common comp,$
pupil,r,x1,y,t,cx,cy,ax,ay,fo,sa,clx,cly,psfm,nx,ny,scl,sax,say,ashx,ashy,tx,ty
common done,done,crve
common screen,iscr
iscr=0 & done=0 & crve=1 
;case 2
a=[-2.92788,3.20713,6.06396,-2.99452,3.77422,0.113594,0.602555,-4.38223]
; case 3
a=[-2.80042,2.29146,6.30618,-1.81163,3.76222,0.876647,1.48907,-5.04596,$
 -0.309251,-4.998904e-02,-1.02750,-1.65498,-2.624347e-02,-5.816614e-02]
; case 4
a=[-3.37238,1.86156,-6.99383,-1.49684,5.65420,-1.15969,-1.40797,6.78794,$
-0.380177,-4.747110e-02,0.923253,1.50126,2.724005e-02, 8.259022e-02]
; case 1
a=[0.538994,0.704896,2.79230,-0.230412,0.580961]
convert,a,ac
ac=[-2.6453,0.08973,0.52430,-0.10045,0.0789783,-0.029463,-0.1011,0.1038]
crve=1
reads,'psfs',psfm,nx,ny
set_plot,'PS'
fun,x,ac,f
stop
end