RE: [Axiom-math] special functions

[yigal]

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

The form is not so good yet but it is better and a lot less taxing on
the eyes then previous versions.

Best,

Yigal