PST自动计算程序

很多人比较关心怎样自动计算PST，提供两种方法（另一种明天提供），目前提供的方法是直接计算PST曲线，程序如下，直接打开ASAP的代码窗口，复制进行即可运行，如果需要计算其它系统的PST，替换雷竞技下载找ray666点vip系统部分即可。

SYSTEM NEW
RESET

UNITS M
WAVELENGTH 10 UM
MEDIA; (N=2.4) =ZnSe
MODELS
  LAMBERTIAN .0001
  LAMBERTIAN .01
  LAMBERTIAN .0001*N^2
EDGES
  RECTANGLE Z -5 2@.001 !primary image
  RECTANGLE Z .5 2@.004 !final image
ENT OBJECTS
  OPTICAL Z -4.05 ELLIP 1.04 =Window; INTERFACE 0 1 0 "ZnSe"
    SCATTER MODEL 3; TOWARDS Z 1 (.001/5)/N
  OPTICAL Z -3.95 ELLIP 1.04 =Window.2; INTERFACE 0 1 "ZnSe" 0
    SCATTER MODEL 1; TOWARDS Z 1 (.001/5)
  OPTICAL Z 0 -10 -1.0353 ELLIP 2@1 .31914 =Primary Mirror
    INTERFACE 1
    SCATTER MODEL 1; TOWARDS 1
  OPTICAL Z -3.9 -2.9333 -3.1576 ELLIP .24044 =Secondary Mirror
    INTERFACE 1
    SCATTER MODEL 1; TOWARDS 2
  PLANE Z .5 ELLIP .1 =Image
  TUBE Z -5 2@1.05 -.05 2@1 =Main Baffle Tube
    SCATTER MODEL 2; TOWARDS 1
  TUBE Z -3.9098 2@.24044 -3.4231 2@.33623 =InnerSecondaryBaffle
    SCATTER MODEL 2; TOWARDS 1; TOWARDS 2
  TUBE Z -2.0214 2@.18019 -.0051 2@.31914 =Outer Conical Baffle
    SCATTER MODEL 2; TOWARDS 1
  TUBE Z -2.0214 2@.18019 .5 2@.1 =Inner Conical Baffle
    SCATTER MODEL 2; TOWARDS 2; TOWARDS 1
RETURN
PROFILE
BEAMS INC GEO; HALT 1.E-19
M=10000
$ITER ANGLE 1 29 -29 PSNIT
{ GRID ELLIP Z -5 -4@1.05 2@SQRT(4*M/3.1416)
    SOURCE DIR Z (ANGLE) 90; FLUX 0 COS[ANGLE]
   TRACE
   CONSIDER ONLY IMAGE; SELECT ONLY SOURCE -1; PATHS
   SPOTS POS ATT 0 YX PIXEL 1; DISPLAY; $GRAB 'Maximum' PSNIT
   CONSIDER ALL }


DISPLAY ITER
  GRAPH
RETURN

还有一种PST的算法在哪儿呢？有链接吗？

在这个贴子上。
http://www.optzmx.com/thread-12529-1-1.html