I made a small mistake but numerically important in my sloppy incomplete gamma code here is a modified version that seems like its working:
x:Float;x:=20.0 a:Integer;a:=7 n : (Integer) -> List Polynomial Float; n(y) == [y-1,[EMAIL PROTECTED] Polynomial Float; numg : (Integer) -> List Polynomial Float; numg(1)== [1,1-a]; numg(i|i>1) == append(numg(i-1),n(i))@List Polynomial Float; d : (Integer) -> List Polynomial Float; d(y) == [x,[EMAIL PROTECTED] Polynomial Float; deng : (Integer) -> List Polynomial Float; deng(1)== [x,1]; deng(i|i>1) == append(deng(i-1),d(i))@List Polynomial Float; num := [numg.i.i for i in 1..]; den := [deng.i.i for i in 1..]; cf := continuedFraction(0,num,den) ccf := convergents cf; gam(i) == exp(-x)*x^a*ccf.i; gamma(n,x) == factorial(n-1)*exp(-x)*reduce(+,[x^i/factorial(i) for i in 0..(n-1)]); where gam.15 = 0.1836881970 165365277 gamma(7,20.)=0.1836881970 165365277 thank you all for your previous help, Yigal Weinstein _______________________________________________ Axiom-math mailing list Axiom-math@nongnu.org http://lists.nongnu.org/mailman/listinfo/axiom-math