[求助]求椭圆逼近程序

sirius

再次感谢ZZXXQQ的帮助 90°鞠躬

PEPSIMAX

好！

guoyan_sw

源码

xyp1964

(defun c:test ()
  (defun ELL (a y b) (* a (sqrt (- 1 (/ (* y y) (* b b) 1.)))))
  (setvar "osmode" 0)
  (setq a  400
 b  200
 y  200
 dy 1
  )
  (command "pline")
  (while (<= (abs y) b)
    (setq x  (ELL a y b)
   pt (list x y)
   y  (- y dy)
    )
    (command pt)
  )
  (setq y -200)
  (while (<= (abs y) b)
    (setq x  (ELL a y b)
   pt (list (* x -1) y)
   y  (+ y dy)
    )
    (command pt)
  )
  (command "c")
  (princ)
)