CLOSE WINDOW j3 d=O!  
DELETEMFO ALL M,1Yce%+}  
W = 5 Pe@M_ r  
K= 20 R:SFj!W1  
sample = 80 
#W`
>vd}  
SYSP 11, W*6   |0:&dw?*!  
SYSP 12, 1 ;EstUs3  
apod_factor =1/POWER((1/3),1) pVe@HJy6G  
SYSP 13, apod_factor )%p.v P'p  
FOR I, 1, sample+1, 1 "-JJ6Bk  
INSERTMFO I J0xOB;rd  
SETOPERAND I, 11, “REAY” Isv@V.  
SETOPERAND I, 9, 1 xzF@v>2S+  
norm_pupil_coord = I/sample ypEcjVPD  
X= norm_ pupil_coord*W*3 iyNyj44 H  
# Un-normilized pupil coordinate JzH\_,,  
S= K*SQRT(1-EXPE(-2*POWER(X,2)/POWER(W,2))) T)Z2=5V  
#same as the final equation in the article A
P.WTFf  
SETOPERAND I, 8, -S !FOPFPn  
# set target. Positive input #coordinatehas negative ^yDCX  
# output coordinate (positive lens) thusthe negative sign Sd?:+\bS;  
SETOPERAND I, 7, norm_pupil_coord *k'D%}N:  
# Set Py column of the operand 0] e=  
SETOPERAND I, 2, NSUR() De`)`\U  
# specify the image surface for the REAYvalue calculation \IY)2C<e  
NEXT lVY`^pw?  
OPTIMIZE ($:s}_<>s  
UPDATE ALL # Update all windows !J ")TP=  
END