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
_______________________________________________
Axiom-math mailing list
Axiom-math@nongnu.org
http://lists.nongnu.org/mailman/listinfo/axiom-math
