k4=0.5*0.5*16.895979
rvals=(indgen(100)*0.66/99.+0.34)
rms=fltarr(100)
for i=0,99 do begin
   r=rvals(i)
   eps=0.33
   in=(r^indgen(11)-eps^indgen(11))/indgen(11)
   a=-(in(6)*in(4)-in(8)*in(2))/(in(4)^2-in(6)*in(2))
   b=(in(6)^2-in(8)*in(4))/(in(4)^2-in(6)*in(2))
   RMS(i)=k4*sqrt((in(10)+2*a*in(8)+(2*b+a^2)*in(6)+2*a*b*in(4)+b^2*in(2))/in(2))
endfor
set_plot,'ps'
device,filename='weiler.ps'
!mtitle='Apodized pupil'
!xtitle='Outer obscuration ratio'
!ytitle='RMS wavefront error'
plot,rvals,rms*2
!ytitle='Percent light loss - FOC'
plot,rvals,100.-(rvals^2-eps^2)*100./(1.-eps^2)
!ytitle='Percent light loss - WFPC'
plot,rvals,100.-(rvals^2-0.42^2)*100./(1.-0.42^2)
!ytitle='Percent Resolution limit loss'
plot,rvals,100.-rvals*100.
device,/close
set_plot,'X'
stop
end