RE: [Axiom-math] special functions

[Vanuxem Grégory]

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
>