On May 16, Raymond Toy wrote:
--------------------------------------
>Just a note that if you update from git again, you'll have to get
>gamma.lisp and numeric.lisp for this to work. I moved the code for the
>Lentz algorithm from gamma.lisp to numeric.lisp.
-------------------------
The new version is behaving nicely:
--------------------------------
(%i1) display2d:false$
(%i2) load("c:/work2/numeric-new.o");
(%o2) "c:/work2/numeric-new.o"
(%i3) load("c:/work2/gamma-new.o");
(%o3) "c:/work2/gamma-new.o"
(%i4) domain:complex;
(%o4) complex
(%i5) ee : integrate(exp(x^5),x,1,2);
(%o5) ((-1)^(4/5)*gamma_incomplete(1/5,-32)-(-1)^(4/5)
*gamma_incomplete(1/5,-1))/5
(%i6) fpprec;
(%o6) 16
(%i7) bfloat(ee);
(%o7)
2.0b-1*(1.0b0*(-1)^(4/5)*(-2.977836145416077b12*%i-4.098639832671666b12)
-1.0b0*(-1)^(4/5)*(-3.599871685328463b0*%i-3.639545917033864b-1))
(%i8) cbfloat(zz):= (expand(bfloat(rectform(zz))))$
(%i9) cbfloat(ee);
(%o9) 1.143595626814653b-5*%i+1.013239489694019b12
(%i10) block([fpprec:20],cbfloat(ee));
(%o10) 1.0132394896940182945b12-6.84948084853770605b-10*%i
---------------------------
Ted