[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:37 -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
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