gamma_incomplete wrong numeric values

gamma_incomplete wrong numeric values

 Example 1
(%i1) integrate(sin(x)/x,x,1,100);
(%o1)
-(%i*gamma_incomplete(0,100*%i))/2+(%i*gamma_incomplete(0,%i))/2-(%i*gamma_incomplete(0,-%i))/2+
     (%i*gamma_incomplete(0,-100*%i))/2
(%i2) float(%), numer$
(%i3) expand(%);
(%o3) 2.102121726749059*10^24

 Example 2
(%i4) integrate(exp(x^4),x,-1,1);
(%o4) -((-1)^(3/4)*gamma_incomplete(1/4,-1)-(-1)^(3/4)*gamma(1/4))/2
(%i5) float(%), numer;
(%o5)
-0.5*(1.0*(-1)^0.75*(.02986697737096389-3.595742930850943*%i)-3.625609908221908*(-1)^0.75)
(%i6) rectform(%)$
(%i7) float(%), numer;
(%o7) 4.440892098500626*10^-16*%i-2.542574209808294

 Example 3
http://www.math.utexas.edu/pipermail/maxima/2011/027242.html
(%i8)
r(n):=expand(float(integrate(integrate(cos(x*y*log(2)),y,1,n),x,1,n)))$
(%i9) makelist(r(n),n,1,31);
(%o9)
[0.0,6.006424694652971*10^-17*%i+.01514411407555013,2.002141564884324*10^-17*%i-1.712607349563953,
2.002141564884324*10^-17*%i-2.037516611552408,-2.058330860423246,-1.7990395639789,1.001070782442162*10^-17*%i-1.282806487348737,
-5.385341733767779,261141.6619068503,5.910683916372842*10^11,3.641253891162185*10^17,2.233034361027964*10^24,
-5.340630456765985*10^31,-1.952422070186322*10^39,-7.020961982238587*10^47,-7.105539688231238*10^57,-7.312154929942063*10^66,
-1.124721377265298*10^78,-8.857175718374467*10^88,-1.31871834594326*10^100,1.06040271806534*10^113,2.63840867785563*10^125,
-6.08147852617669*10^139,-3.32111833808186*10^153,6.589486561463*10^167,8.42086569403558*10^183,1.08007233283605*10^200,
-4.1944716083056*10^214,2.17161643827435*10^233,-2.90547722801773*10^250,9.63249577751765*10^268]
(%i10) r(32);
expintegral_e: series failed.
#0: r(n=32)
 -- an error. To debug this try: debugmode(true);

(%i11) build_info()$
Maxima version: 5.26.0
Maxima build date: 18:42 12/19/2011
Host type: i686-pc-linux-gnu
Lisp implementation type: SBCL
Lisp implementation version: 1.0.50.0.debian

With best wishes for a Happy New Year

Aleksas D