请各位高手帮忙改一下生成混凝土搅拌车拌筒和螺旋线的GRIP程序的错误

DUHUIDUHUI

编译和链接都没有错误，只是在载入ug中时报错说数组超出范围，错误在第48行，已用红色标记
5 `- O  g; N0 j' w: V/ P0 lNUMBER/r(3),h(3),bt(3),i,st,stm(3),al,a3,p1,p3
) e' ^+ H2 w2 a                                           $$1-定义变量
0 m: X& a' `9 e: @* p. Y( GENTITY/ax,ps1(1001),ps2(1001),m(3),g(3),sl1,$
" k* k5 d  [% }) P  l2 e; b9 e      sl2,pt(4),pt0(2)
4 _3 I7 {. K/ Hpds1:                                      $$2-生成输入参数的对话框
    PARAM/'请输入罐体和螺旋线参数:',$
) h4 Z1 j% q1 Q! l& L      '前锥小端半径:',r(1),$
      '前锥高度:',h(1),$
4 @# y1 ~  B( M- q3 x, e1 g      '中柱半径:',r(2),$
      '中柱高度:',h(2),$
" M5 z9 m4 t9 [, f& ]; ^      '后锥小端半径:',r(3),$
      '后锥高度:',h(3),$
      '前锥螺旋角:',bt(1),$
      '中锥螺旋角:',bt(2),$
( I: ], k( N, w& n  r0 e6 G& b! _) t" {      '后锥螺旋角:',bt(3),rsp
- j1 F/ I/ }- l& \   JUMP/pds1:,stop1:,rsp
* j, N( S# N# S) }   ax=LINE/(pt0(1)=POINT/0,0,0),(pt0(2)=POINT/0,0,600)
3 H0 ~' D) w# v# C% Z* A# z# C8 _4 I                                           $$3-生成罐体
0 X/ G' X! k; N" U5 o* W   g(1)=LINE/(pt(1)=POINT/r(1),0,0),(pt(2)=POINT/r(2),0,h(1))
+ J* {) P5 l( l, w& |7 T   g(2)=LINE/pt(2),(pt(3)=POINT/r(2),0,(h(1)+h(2)))
0 Q: F0 \2 C+ S5 i- I  O% D   g(3)=LINE/PT(3),(pt(4)=POINT/r(3),0,(h(1)+h(2)+h(3)))
   m(1)=REVSRF/g(1),AXIS,ax,0,360
   m(2)=REVSRF/g(2),AXIS,ax,0,360
! I% w2 J- f4 R; _   m(3)=REVSRF/g(3),AXIS,ax,0,360
   a1=ATANF((r(2)-r(1)/h(1)))          $$4-计算基本参数
( O! s) O+ \: J; y. O) N   a3=ATANF((r(3)-r(2)/h(3)))
5 }* ~: w  k0 i   p1=(r(2)-r(1))/h(1)*COSF(a1)*COSF(bt(1))/SINF(bt(1))
   p3=(r(3)-r(2))/h(3)*COSF(a3)*COSF(bt(3))/SINF(bt(3))
   stm(1)=LOGF((r(2)-r(1))/r(1)+1)/p1
   stm(2)=stm(1)+h(2)/r(2)*SINF(bt(2))/cosf(bt(2))
* c0 T9 v- a  Z   stm(3)=stm(2)+LOGF((r(3)-r(2))/r(2)+1)/p3
   i=1
  u6 ^' U7 d/ [0 U# N+ V   STR1:                              $$5-计算曲线上点的坐标
   st=(i-1)*stm(3)/1000
6 {! x" k) U- Z; Y; O  IFTHEN/st<=stm(1)
' z0 s4 k& c, N8 I5 f! F   xc=r(1)*EXPF(p1*st)*COSF(st*180/3.14159)
   yc=r(1)*EXPF(p1*st)*SINF(st*180/3.14159)
   zc=r(1)*h(1)/(r(2)-r(1))*(EXPF(pl*st)-1)
; D7 j% F6 N2 l7 Y  ELSEIF/stm(1)<st AND st<=stm(2)
/ H5 {# d3 n' S& v+ }   xc=r(2)*COSF(st*180/3.14159)
8 f+ d) K( t" g" z9 F) ~! w: m   yc=r(2)*SINF(st*180/3.14159)
" w& G$ _- L( W6 l! [' `: R" K   zc=r(2)*COSF(bt(2))/SINF(bt(2))*(st-stm(1))+h(1)
  ELSEIF/stm(2)<st AND st<=stm(3)
8 I* f' o6 e  `1 m   xc=r(2)*EXPF(p3*(st-stm(2)))*COSF(st*180/3.14159)
   yc=r(2)*EXPF(p3*(st-stm(2)))*SINF(st*180/3.14159)
   zc=r(2)*h(3)/(r(3)-r(2))*(EXPF(p3*(st-stm(2)))-1)+h(1)+h(2)
  ENDIF
2 `4 N/ j2 [( e( s+ b! |; Z( d9 p5 C   ps1(i)=POINT/xc,yc,zc            $$6-生成螺旋线
   ps2(i)=POINT/-xc,-yc,zc
   i=i+1
  IF/st<=stm(3),JUMP/STR1:
  sl1=SPLINE/ps1(1..(i-1))
  sl2=SPLINE/ps2(1..(i-1))
BLANK/ax,g,pt,pt0,ps1(1..(i-1)),ps1(1..(i-1))
                                       $$结束处理
: w' X, n  u0 e& J. c5 E3 \% K6 [  stop1:
   HALT

27699885

亲自测试过这个程序，对数螺旋线grip吧，也曾出现过这问题qq137303928