Hi, > -----Message d'origine----- > De : [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] la part de > yigal > Envoyé : dimanche 29 janvier 2006 23:19 > À : [email protected] > Objet : RE: [Axiom-math] special functions > > > 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
Yes but: 9) -> gamcfe(7.0,0.0) Compiling function nume with type Float -> Stream Float Compiling function dene with type (Float,Float) -> Stream Float Compiling function cfe with type (Float,Float) -> ContinuedFraction Float Compiling function ccfe with type (Float,Float) -> Stream Fraction Float >> Error detected within library code: division by zero Cheers, Greg > > 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 > [email protected] > http://lists.nongnu.org/mailman/listinfo/axiom-math > _______________________________________________ Axiom-math mailing list [email protected] http://lists.nongnu.org/mailman/listinfo/axiom-math
