[求助]我绘制的日照曲线，参数有问题！！

MRHDF

(defun c: rz ()
 
  (setq　pt1 (getpoint "\请输入基点："))
  　
  (setq　an1　 (dtr ( getreal　 "地理纬度")))
  　
  (setq　an2　 (dtr (getreal "大寒日或冬至日的赤纬角")))
  (setq　an3　0.0)
  　
  (setq　ta　 (/　 (*　5　pi) 　180))
  (command　 "ucs" 　 "o" 　pt1)
  (command　 "pline" 　 (list　r　0.0　0.0))
  (repeat　16
    (setq r
    (sqrt
      (-
        1000000.0
        (*
   (/
     1000.0
     ((+ (* (sin an1) (sin an2)) (* (cos an1) (cos an2) (cos an3)))
     )
   )
   (/
     1000.0
     ((+ (* (sin an1) (sin an2)) (* (cos an1) (cos an2) (cos an3)))
     )
   )
        )
      )
    )
    )
    (setq　y　 (*　r　 (cos　an3)))
    (setq　x　 (*　r　 (sin　an3)))
    (setq　pt2　 (list　x　y))
    (command　pt2)
    (setq　an3　 (+　ta　an3))
  )
)
;;;弧度转换为角度的设置数
(defun dtr (a) (* a (/ pi 180.0)))

可以用 an1=32    an2=20。  合肥大寒日。

问题： 说我参数太多。可是每句话单独执行还都没有问题。

yshf

程序中有空格作分隔的地方全部输入为中文状态（全角）的空格，将其改为英文状态的空格即可。另外：

“(/
     1000.0
     ((+ (* (sin an1) (sin an2)) (* (cos an1) (cos an2) (cos an3)))
     )
   )
   (/
     1000.0
     ((+ (* (sin an1) (sin an2)) (* (cos an1) (cos an2) (cos an3)))
     )
   )”

中“+”号多了“(”，后面又多了“)”，应改为：

“(/ 1000.0 (+ (* (sin an1) (sin an2)) (* (cos an1) (cos an2) (cos an3))))
              (/ 1000.0 (+ (* (sin an1) (sin an2)) (* (cos an1) (cos an2) (cos an3))))”

但开方数出现了负数，可能计算公式输入有误，请仔细查正。改过的程序一同附上。

MRHDF

这这么快就回了  太感谢你了。负数是我搞错了：减法的被减数和减数反了这个。

曲线我们搞小区的时候常用，不懂得地方要住店输入很麻烦。

MRHDF

曲线还是不像样子   这回应该是我公式的问题了。