Re: [Axiom-mail] Getting numerical values from a closed form integration?

2007-07-12 Thread Francois Maltey
 For example:
 
 integrate(sin(2*x),x=-%pi/3..4*%pi/7)
 
 The answer is in closed form, but I can't convert it to a float.  At least:
 
 integrate(sin(2*x),x=-%pi/3..4*%pi/7)::Float

Try 

numeric integrate(sin(2*x),x=-%pi/3..4*%pi/7)

I get 0.20048...

Don't you find a strong typed CAS is nice for mathematics
even if the final user is a little bored...

Francois


___
Axiom-mail mailing list
Axiom-mail@nongnu.org
http://lists.nongnu.org/mailman/listinfo/axiom-mail


Re: [Axiom-mail] Getting numerical values from a closed form integration?

2007-07-12 Thread Bill Page

On 12 Jul 2007 12:58:49 +0200, Francois Maltey wrote:

 For example:

 integrate(sin(2*x),x=-%pi/3..4*%pi/7)

 The answer is in closed form, but I can't convert it to a float.  At least:

 integrate(sin(2*x),x=-%pi/3..4*%pi/7)::Float

Try

numeric integrate(sin(2*x),x=-%pi/3..4*%pi/7)

I get 0.20048...


I think that is a good solution. But one thing to keep in mind perhaps
for later use is that often when

 integrate(sin(2*x),x=-%pi/3..4*%pi/7)::Float

fails, the following similar coercion will work

 integrate(sin(2*x),x=-%pi/3..4*%pi/7)::Expression Float

I think the issue is a certain weakness in Axiom's system of coercions.


Don't you find a strong typed CAS is nice for mathematics
even if the final user is a little bored...



Francois, could you expand this comment? I am not sure what point you
are making when you say bored. :-)

Regards,
Bill Page.


___
Axiom-mail mailing list
Axiom-mail@nongnu.org
http://lists.nongnu.org/mailman/listinfo/axiom-mail