login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

Submitted by : (unknown) at: 2007-11-17T22:28:24-08:00 (16 years ago)
Name :
Axiom Version :
Category : Severity : Status :
Optional subject :  
Optional comment :

The first real calculus is right, the second fails.

axiom
asin 1.0

\label{eq1}1.5707963267948966192(1)
Type: Float
axiom
asin (1.0+0*%i)
>> Error detected within library code: catdef: division by zero

The problem is perhaps line 162 of trigcat.spad. There is no test for the ... / sqrt (1-x^2) in asin definition.

asin x == atan(x/sqrt(1-x**2))

patch possible. --Francois Maltey, Fri, 05 Jan 2007 07:57:48 -0600 reply
I propose to remplace in asin x == atan (x/sqrt(1-x**2)) in trigcat.spad line 162 by

asin x == x=1 => 2atan 1 x=-1 => -2atan 1 atan(x/sqrt(1-x**2))

The result seems right...

Category: Axiom Mathematics => Axiom Library Status: open => duplicate




  Subject:   Be Bold !!
  ( 14 subscribers )  
Please rate this page: