;;由 sin 值返回角度(弧度)，返回值位于 -pi/2 ~ pi/2
(defun csin (num / n tan)
  (if (and (<= num 1.)
	   (>= num -1.)
      )
    (progn
      (setq n	(sqrt (- 1. (* num num)))
	    tan	(/ n num)
      )
      (atan tan)
    )
  )
)