zemax源碼：
This DLL models an anamorphic aspheric surface.
This surface is essentially an even aspheric surface with different terms for
the x and y directions.
The sag is given by:
Z = ((CX*x*x)+(CY*y*y)) / (1 + sqrt(1-((1+KX)*CX*CX*x*x)-((1+KY)*CY*CY*y*y)))
+ AR*( (1 - AP)*x*x + (1 + AP)*y*y )^2
+ BR*( (1 - BP)*x*x + (1 + BP)*y*y )^3
+ CR*( (1 - CP)*x*x + (1 + CP)*y*y )^4
+ DR*( (1 - DP)*x*x + (1 + DP)*y*y )^5
Note the terms AR, BR, CR, and DR ... have units of length to the -3, -5, -7, and -9 power.
The terms AP, BP, CP, and DP are dimensionless.
The surface is rotationally symmetric only if AP = BP = CP = DP == 0 and CX = CY and KX = KY.
zemax源碼：
This DLL models a standard ZEMAX surface type, either plane, sphere, or conic
The surface also demonstrates a user-defined apodization filter
The filter is defined as part of the real ray trace, case 5
The filter can be used at the stop to produce x-y Gaussian apodization similar to the Gaussian pupil apodization in ZEMAX but separate in x and y.
The amplitude apodization is of the form EXP[-(Gx(x/R)^2 + Gy(y/R)^2)]
The transmission is of the form EXP[-2(Gx(x/R)^2 + Gy(y/R)^2)]
where
x^2 + y^2 = r^2
R = semi-diameter
The tranmitted intensity is maximum in the center.
T is set to 0 if semi-diameter < 1e-10 to avoid division by zero.
