[Axiom-developer] cot and asec

[Juergen Weiss]

Here the lisp version of the C versions in csl
One problem: the formula may work for real arguments
only.


;; stolen from csl
#+(or :cmu :akcl)
(defun cot (a)
  (if (or (> a 1000.0) (< a -1000.0))
      (/ (cos a) (sin a))
    (/ 1.0 (tan a))))

;; stolen from csl
#+(or :cmu :akcl)
(defun asec (a)
  (if (and (a > -1.0) (a < 1.0))
      0.0
    (acos (/ 1.0 a))))

-- 
Juergen Weiss     | Universitaet Mainz, Zentrum fuer Datenverarbeitung,
address@hidden| 55099 Mainz, Tel: +49(6131)39-26361, FAX: +49(6131)39-26407