[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: [Axiom-math] special functions
From: |
Vanuxem Grégory |
Subject: |
RE: [Axiom-math] special functions |
Date: |
Mon, 30 Jan 2006 01:31:05 +0100 |
Hi,
> -----Message d'origine-----
> De : address@hidden
> [mailto:address@hidden la part de
> yigal
> Envoyé : dimanche 29 janvier 2006 23:19
> À : address@hidden
> Objet : RE: [Axiom-math] special functions
>
>
> The even parts of the continued fraction do work after specifying that
> all numbers should be converted to floats. So here is my latest rewrite
> of the incomplete gamma function,
>
> )clear all
> j:=120
> nume(a) == cons(1 :: Float,[((a-i)*i) :: Float for i in 1..]);
> dene(a,x) == [(x+2*i+1-a) :: Float for i in 0..];
> cfe(a,x) == continuedFraction(0,nume(a),dene(a,x));
> ccfe(a,x) == convergents cfe(a,x);
> gamcfe(a,x) == exp(-x)*x^a*ccfe(a,x).j;
>
>
> (24) -> gamcfe(0.0,1.0)*exp(1.0)
> (24) ->
> (24) 0.5963473623 2319407232
> Type: Fraction Float
Yes but:
9) -> gamcfe(7.0,0.0)
Compiling function nume with type Float -> Stream Float
Compiling function dene with type (Float,Float) -> Stream Float
Compiling function cfe with type (Float,Float) -> ContinuedFraction Float
Compiling function ccfe with type (Float,Float) -> Stream Fraction Float
>> Error detected within library code:
division by zero
Cheers,
Greg
>
> The form is not so good yet but it is better and a lot less taxing on
> the eyes then previous versions.
>
> Best,
>
> Yigal
>
>
>
> _______________________________________________
> Axiom-math mailing list
> address@hidden
> http://lists.nongnu.org/mailman/listinfo/axiom-math
>