[教学]：请各位大虾帮我看一哈摆线轮的这个程序

pearce

老师最近要我用AutoLISP编写一个摆线方程的程序，进行摆线轮的绘制：下面的程序尚未编辑，请各位仁兄帮我检查，修改，先行谢过，
；角度与弧度转换函数
(defun C:dtr(a)                                                      
        (* pi (/ a 180.0))
         )
；主函数
(defun C:main()
；系统准备                                                   
        (setq sblip ( getvar "blipmode" ) )
        (setq scmde ( getvar "cmdecho" ) )
        (setvar "blipmode" 0 )
        (setvar "cmdecho" 0 )                                            
；参数输入
        (setq r2 ( getreal "\n输入针轮半径："))
        (setq r3 ( getreal "\n输入摆线轮半径："))
        (setq k ( getreal "\n输入短幅系数："))
        (setq zb ( getreal "\n输入针轮齿数："))
        (setq zg ( getreal "\n输入摆线轮齿数："))
        (setq d ( getreal "\n输入计算步长："))  
；参数转化
        (setq s 0.0)
        (setq s ( dtr s ))
        (setq d (dtr d ))                                                
；寻找摆线上东方零度起点
        (setq  l (* r3 (- (sin s ) (* (/ k zb )(sin (* zb s))))))
        (setq  m (* r3 (- (cos s ) (* (/ k zb )(cos (* zb s))))))
        (setq  w (sqrt ( - (+ (expt k 2) 1) (* k 2 (cos(* zg s ))))))
        (setq  o (/ (- (* k (sin (* zb s))) (sin s)) w))
        (setq  p (/ (- (* -1 k (cos(* zb s ))) (sin s )) w))
        (setq  x (+ l (* r2 o)))
        (setq  y (- m (* r2 p)))
        (setq  u (list x y ))                                             
；迭代循环绘制摆线轮廓上的点
        (command "pline" u)
        (while (< s (+ pi 0.001))
                (setq  l (* r3 (- (sin s ) (* (/ k zb )(sin (* zb s))))))
                (setq  m (* r3 (- (cos s ) (* (/ k zb )(cos (* zb s))))))
                (setq  w (sqrt ( - (+ (expt k 2) 1) (* k 2 (cos(* zg s ))))))
                (setq  o (/ (- (* k (sin (* zb s))) (sin s)) w))
                (setq  p (/ (- (* -1 k (cos(* zb s ))) (sin s )) w))
                (setq  x (+ l (* r2 o)))
                (setq  y (- m (* r2 p)))
                (setq  u (list x y ))
                (setq  s (+ s d))
                (command u)
        )                                                               
；两维面内的曲线光滑拟合
        (command " " "pedit" "last" "f" " ")                             
； 恢复系统原来状态                  
        (setvar "blipmode" sblip)
        (setvar "cmdecho" scmde)                                       
        (princ)
  )                                                                     
;程序结束
   
；加载操作：
command:(load "dtr(a))
command:(load "main")
command:main

zxq0220

函数：
(defun C:dtr(a)
(* pi (/ a 180.0))
)
应定义成：
(defun dtr(a)
(* pi (/ a 180.0))
)

pearce

谢过二楼楼主，我还有一个问题，也是这个程序的关键就是，寻找摆线轮上的起点的问题，不知道我的编写是否有误，关于这条曲线的方程我上传了一下但是，传不上来，下面是我的联系方式：song.cs2001@163.com  QQ:276048553