_ti='Fm250ki 10 in f9.63'
fo=479.772; 479.736
pt=5.33929; 5.34
st=-17.3471; -17.35
tt=-13.0171; -13.01
kp=-0.928398; -0.9055
rq=4370.63; 4369.949
;
rp=-4935
ps=1200
rs=-6350
ss=800
rt=rs
ts=-312.5
ks=0
kt=0
kq=0
qt=-45-pt-st-tt
;
print,'qt=',qt
print,'kp [tct]=',kp*kp
;
objp,-74*12.17e6,0*11.80e6,206264.812e6,1,6328
pupca,125,8,16
disp,0,0,ps
rot,0,pt
conic,kp,rp
mirror
rot,0,pt
disp,0,0,-ps
rot,0,st
conic,ks,rs
mirror
rot,0,st
disp,0,0,ss
rot,0,tt
conic,kt,rt
mirror
rot,0,tt
disp,0,0,ts
rot,0,qt
conic,kq,rq
mirror
rot,0,qt
disp,0,0,fo
flat
detector
spot ;
m1=rms(0)

objp,74*12.17e6,0*11.80e6,206264.812e6,1,6328
pupca,125,8,16
disp,0,0,ps
rot,0,pt
conic,kp,rp
mirror
rot,0,pt
disp,0,0,-ps
rot,0,st
conic,ks,rs
mirror
rot,0,st
disp,0,0,ss
rot,0,tt
conic,kt,rt
mirror
rot,0,tt
disp,0,0,ts
rot,0,qt
conic,kq,rq
mirror
rot,0,qt
disp,0,0,fo
flat
detector
spot ;
m3=rms(0)

objp,0*12.17e6,74*11.80e6,206264.812e6,1,6328
pupca,125,8,16
disp,0,0,ps
rot,0,pt
conic,kp,rp
mirror
rot,0,pt
disp,0,0,-ps
rot,0,st
conic,ks,rs
mirror
rot,0,st
disp,0,0,ss
rot,0,tt
conic,kt,rt
mirror
rot,0,tt
disp,0,0,ts
rot,0,qt
conic,kq,rq
mirror
rot,0,qt
disp,0,0,fo
flat
detector
spot ;
m4=rms(0)

objp,0*12.17e6,0*11.80e6,206264.812e6,1,6328
pupca,125,13,47
disp,0,0,ps
rot,0,pt
conic,kp,rp
mirror
rot,0,pt
disp,0,0,-ps
rot,0,st
conic,ks,rs
mirror
rot,0,st
disp,0,0,ss
rot,0,tt
conic,kt,rt
mirror
rot,0,tt
disp,0,0,ts
rot,0,qt
conic,kq,rq
mirror
rot,0,qt
disp,0,0,fo
flat
detector
spot ;
m2=rms(0)
print,'m1=',m1,', m2=',m2, ', m3=',m3,', m4=',m4

merit=m1*m1+1000*m2*m2+m3*m3+2*m4*m4
merit=merit/13
print,'merit=',merit
end