[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: [Axiom-math] special functions

From: yigal
Subject: RE: [Axiom-math] special functions
Date: Sun, 29 Jan 2006 14:18:21 -0800

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 
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.



reply via email to

[Prev in Thread] Current Thread [Next in Thread]