You can subscribe to this list here.
2002 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}
(67) 
_{Jul}
(61) 
_{Aug}
(49) 
_{Sep}
(43) 
_{Oct}
(59) 
_{Nov}
(24) 
_{Dec}
(18) 

2003 
_{Jan}
(34) 
_{Feb}
(35) 
_{Mar}
(72) 
_{Apr}
(42) 
_{May}
(46) 
_{Jun}
(15) 
_{Jul}
(64) 
_{Aug}
(62) 
_{Sep}
(22) 
_{Oct}
(41) 
_{Nov}
(57) 
_{Dec}
(56) 
2004 
_{Jan}
(48) 
_{Feb}
(47) 
_{Mar}
(33) 
_{Apr}
(39) 
_{May}
(6) 
_{Jun}
(17) 
_{Jul}
(19) 
_{Aug}
(10) 
_{Sep}
(14) 
_{Oct}
(74) 
_{Nov}
(80) 
_{Dec}
(22) 
2005 
_{Jan}
(43) 
_{Feb}
(33) 
_{Mar}
(52) 
_{Apr}
(74) 
_{May}
(32) 
_{Jun}
(58) 
_{Jul}
(18) 
_{Aug}
(41) 
_{Sep}
(71) 
_{Oct}
(28) 
_{Nov}
(65) 
_{Dec}
(68) 
2006 
_{Jan}
(54) 
_{Feb}
(37) 
_{Mar}
(82) 
_{Apr}
(211) 
_{May}
(69) 
_{Jun}
(75) 
_{Jul}
(279) 
_{Aug}
(139) 
_{Sep}
(135) 
_{Oct}
(58) 
_{Nov}
(81) 
_{Dec}
(78) 
2007 
_{Jan}
(141) 
_{Feb}
(134) 
_{Mar}
(65) 
_{Apr}
(49) 
_{May}
(61) 
_{Jun}
(90) 
_{Jul}
(72) 
_{Aug}
(53) 
_{Sep}
(86) 
_{Oct}
(61) 
_{Nov}
(62) 
_{Dec}
(101) 
2008 
_{Jan}
(100) 
_{Feb}
(66) 
_{Mar}
(76) 
_{Apr}
(95) 
_{May}
(77) 
_{Jun}
(93) 
_{Jul}
(103) 
_{Aug}
(76) 
_{Sep}
(42) 
_{Oct}
(55) 
_{Nov}
(44) 
_{Dec}
(75) 
2009 
_{Jan}
(103) 
_{Feb}
(105) 
_{Mar}
(121) 
_{Apr}
(59) 
_{May}
(103) 
_{Jun}
(82) 
_{Jul}
(67) 
_{Aug}
(76) 
_{Sep}
(85) 
_{Oct}
(75) 
_{Nov}
(181) 
_{Dec}
(133) 
2010 
_{Jan}
(107) 
_{Feb}
(116) 
_{Mar}
(145) 
_{Apr}
(89) 
_{May}
(138) 
_{Jun}
(85) 
_{Jul}
(82) 
_{Aug}
(111) 
_{Sep}
(70) 
_{Oct}
(83) 
_{Nov}
(60) 
_{Dec}
(16) 
2011 
_{Jan}
(61) 
_{Feb}
(16) 
_{Mar}
(52) 
_{Apr}
(41) 
_{May}
(34) 
_{Jun}
(41) 
_{Jul}
(57) 
_{Aug}
(73) 
_{Sep}
(21) 
_{Oct}
(45) 
_{Nov}
(50) 
_{Dec}
(28) 
2012 
_{Jan}
(70) 
_{Feb}
(36) 
_{Mar}
(71) 
_{Apr}
(29) 
_{May}
(48) 
_{Jun}
(61) 
_{Jul}
(44) 
_{Aug}
(54) 
_{Sep}
(20) 
_{Oct}
(28) 
_{Nov}
(41) 
_{Dec}
(137) 
2013 
_{Jan}
(62) 
_{Feb}
(55) 
_{Mar}
(31) 
_{Apr}
(23) 
_{May}
(54) 
_{Jun}
(54) 
_{Jul}
(90) 
_{Aug}
(46) 
_{Sep}
(38) 
_{Oct}
(60) 
_{Nov}
(92) 
_{Dec}
(17) 
2014 
_{Jan}
(62) 
_{Feb}
(35) 
_{Mar}
(72) 
_{Apr}
(30) 
_{May}
(97) 
_{Jun}
(81) 
_{Jul}
(63) 
_{Aug}
(64) 
_{Sep}
(28) 
_{Oct}
(45) 
_{Nov}
(48) 
_{Dec}
(109) 
2015 
_{Jan}
(106) 
_{Feb}
(36) 
_{Mar}
(56) 
_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 


1
(3) 
2
(2) 
3
(9) 
4
(7) 
5
(4) 
6
(1) 
7
(1) 
8

9
(2) 
10
(2) 
11
(3) 
12
(2) 
13
(1) 
14
(2) 
15
(1) 
16

17
(1) 
18
(1) 
19
(2) 
20
(1) 
21
(2) 
22
(5) 
23
(2) 
24
(7) 
25
(6) 
26
(1) 
27
(3) 
28
(1) 
29

30
(1) 
31




From: SourceForge.net <noreply@so...>  20110822 12:18:16

Bugs item #3377380, was opened at 20110725 14:04 Message generated for change (Comment added) made by dloksnel You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377380&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Solving equations Group: None Status: Open Resolution: Remind Priority: 7 Private: No Submitted By: Art Lenskold (dloksnel) >Assigned to: Barton Willis (willisbl) Summary: 7 nested levels Initial Comment: solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2(b*xt*sqrt(1x^2))*(bsqrt(1x^2))2*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)))^1*(a+(11/(4*((tx)^2+(bsqrt(1x^2))^2)^2)*(4*(b*xt*sqrt(1x^2))^2*(bsqrt(1x^2))^28*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))+4*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))))^(1/2)))^2 ], [x] ); Maxima version: 5.24.0 Maxima build date: 20:39 4/5/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  >Comment By: Art Lenskold (dloksnel) Date: 20110822 08:18 Message: Still fails after my correction below: solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))2*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)))^1*(a+(11/(4*(tx)^2+(bsqrt(1x^2))^2)^2)*(4*(b*xt*sqrt(1x^2))^2*(bsqrt(1x^2))^28*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))+4*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))))^(1/2)))^2], [x]);  Comment By: Art Lenskold (dloksnel) Date: 20110820 15:55 Message: solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))2*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)))^1*(a+(11/(4*(tx)^2+(bsqrt(1x^2))^2)^2)*(4*(b*xt*sqrt(1x^2))^2*(bsqrt(1x^2))^28*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))+4*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))))^(1/2)))^2], [x]);  Comment By: Art Lenskold (dloksnel) Date: 20110802 19:18 Message: Thank you for correcting an error on my part. Did you use a diagnostic tool or was it the result of eyeballing the error ? Attached is a somewhat improved corrected expression.  Comment By: Barton Willis (willisbl) Date: 20110801 20:38 Message: The (2(b*x ... makes this an invalid expression solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2(b*x < Maybe you could post a corrected expression  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377380&group_id=4933 
From: SourceForge.net <noreply@so...>  20110821 22:02:34

Bugs item #3395683, was opened at 20110821 15:29 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3395683&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Simplification Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Nelson Luís Dias (nldias) Assigned to: Nobody/Anonymous (nobody) Summary: as of Maxima5.21.1, depends() does not allow a similarity t Initial Comment: Hi: my name is Nelson Dias, and I am an Associate Professor in Environmental Engineering at the Federal University of Paraná, Brazil. Below is a classical problem: it transforms the PDE (dPsi/dy)(d^2Psi/dxdy)  (dPsi/dx)(d^2Psi/dy^2)  nu*(d^3Psi/dy^3) = 0 into the ODE 2f''' + f f'' = 0 through the similarity transform U = U*delta*eta, delta = (nu*x/U)^(1/2), eta = y/delta. This is the classical Blasius solution for a laminar boundarylayer. The code below transforms the PDE into the ODE: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ delta : (nu*x/U)^(1/2)$ eta : y/delta$ depends(f,eta)$ psi : U*delta*f$ eq : diff(psi,y) * diff(psi,x,1,y,1)  diff(psi,x)*diff(psi,y,2)  nu*diff(psi,y,3)$ factor(eq); ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ In Maxima 5.20.1 it gives the right answer: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ (%i1) batch(psibla.max) batching /home/nldias/software/maxima5.20.1/psibla.max nu x 1/2 (%i2) delta : () U y (%i3) eta :  delta (%i4) depends(f, eta) (%i5) psi : U delta f (%i6) eq :  nu diff(psi, y, 3)  diff(psi, x) diff(psi, y, 2) + diff(psi, y) diff(psi, x, 1, y, 1) (%i7) factor(eq) 3 2 d f d f 2 (2  + f ) U y 3 y 2 d() d() nu x nu x sqrt() sqrt() U U (%o7)   2 x ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ but it fails for all the following versions: 5.21.1 5.22.1 5.23.2 5.24.0 5.25.0 giving the following error: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ (%i1) batch(psibla.max) read and interpret file: /home/nldias/software/maxima5.25.0/psibla.max nu x 1/2 (%i2) delta : () U y (%i3) eta :  delta (%i4) depends(f, eta) y depends: argument must be a symbol; found  nu x sqrt() U  an error. To debug this try: debugmode(true); ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ I don't know if this is a bug, but if it is not I would very much like to know what the alternative is for the code I used until Maxima5.20.1. Similarity transforms are a big deal to me, and I would like to be able to continue to use this (and other Maxima examples) in my classes. Sincerely Nelson L. Dias  >Comment By: Dieter Kaiser (crategus) Date: 20110822 00:02 Message: Thank you very much for your bug report. I had a look at the problem. Because of a bug report the routine i$dependencies has been reworked in 14.12.2009. At this time a check has been added to allow only dependencies from a symbol and not an expression. This might be too restrictive, but on the other hand, if we allow dependencies on expressions, we get expressions which have a syntax which is wrong for the function diff. The following is generated with Maxima 5.20: (%i1) depends(f, 2*x); (%o1) [f(2*x)] (%i2) diff(f, x); (%o2) 2*'diff(f,2*x,1) (%i3) ev(%), nouns; diff: second argument must be a variable; found 2*x  an error. To debug this try: debugmode(true); The problem is that the function diff does not allow the differentiation wrt an expression. Perhaps, we can extend the syntax of the function diff to allow the differentiation wrt an expression. I will think again about this problem. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3395683&group_id=4933 
From: SourceForge.net <noreply@so...>  20110821 13:29:12

Bugs item #3395683, was opened at 20110821 11:29 Message generated for change (Tracker Item Submitted) made by nldias You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3395683&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Simplification Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Nelson Luís Dias (nldias) Assigned to: Nobody/Anonymous (nobody) Summary: as of Maxima5.21.1, depends() does not allow a similarity t Initial Comment: Hi: my name is Nelson Dias, and I am an Associate Professor in Environmental Engineering at the Federal University of Paraná, Brazil. Below is a classical problem: it transforms the PDE (dPsi/dy)(d^2Psi/dxdy)  (dPsi/dx)(d^2Psi/dy^2)  nu*(d^3Psi/dy^3) = 0 into the ODE 2f''' + f f'' = 0 through the similarity transform U = U*delta*eta, delta = (nu*x/U)^(1/2), eta = y/delta. This is the classical Blasius solution for a laminar boundarylayer. The code below transforms the PDE into the ODE: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ delta : (nu*x/U)^(1/2)$ eta : y/delta$ depends(f,eta)$ psi : U*delta*f$ eq : diff(psi,y) * diff(psi,x,1,y,1)  diff(psi,x)*diff(psi,y,2)  nu*diff(psi,y,3)$ factor(eq); ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ In Maxima 5.20.1 it gives the right answer: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ (%i1) batch(psibla.max) batching /home/nldias/software/maxima5.20.1/psibla.max nu x 1/2 (%i2) delta : () U y (%i3) eta :  delta (%i4) depends(f, eta) (%i5) psi : U delta f (%i6) eq :  nu diff(psi, y, 3)  diff(psi, x) diff(psi, y, 2) + diff(psi, y) diff(psi, x, 1, y, 1) (%i7) factor(eq) 3 2 d f d f 2 (2  + f ) U y 3 y 2 d() d() nu x nu x sqrt() sqrt() U U (%o7)   2 x ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ but it fails for all the following versions: 5.21.1 5.22.1 5.23.2 5.24.0 5.25.0 giving the following error: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ (%i1) batch(psibla.max) read and interpret file: /home/nldias/software/maxima5.25.0/psibla.max nu x 1/2 (%i2) delta : () U y (%i3) eta :  delta (%i4) depends(f, eta) y depends: argument must be a symbol; found  nu x sqrt() U  an error. To debug this try: debugmode(true); ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ I don't know if this is a bug, but if it is not I would very much like to know what the alternative is for the code I used until Maxima5.20.1. Similarity transforms are a big deal to me, and I would like to be able to continue to use this (and other Maxima examples) in my classes. Sincerely Nelson L. Dias  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3395683&group_id=4933 
From: SourceForge.net <noreply@so...>  20110820 19:55:52

Bugs item #3377380, was opened at 20110725 14:04 Message generated for change (Comment added) made by dloksnel You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377380&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Solving equations Group: None Status: Open >Resolution: Remind >Priority: 7 Private: No Submitted By: Art Lenskold (dloksnel) Assigned to: Nobody/Anonymous (nobody) Summary: 7 nested levels Initial Comment: solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2(b*xt*sqrt(1x^2))*(bsqrt(1x^2))2*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)))^1*(a+(11/(4*((tx)^2+(bsqrt(1x^2))^2)^2)*(4*(b*xt*sqrt(1x^2))^2*(bsqrt(1x^2))^28*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))+4*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))))^(1/2)))^2 ], [x] ); Maxima version: 5.24.0 Maxima build date: 20:39 4/5/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  >Comment By: Art Lenskold (dloksnel) Date: 20110820 15:55 Message: solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))2*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)))^1*(a+(11/(4*(tx)^2+(bsqrt(1x^2))^2)^2)*(4*(b*xt*sqrt(1x^2))^2*(bsqrt(1x^2))^28*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))+4*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))))^(1/2)))^2], [x]);  Comment By: Art Lenskold (dloksnel) Date: 20110802 19:18 Message: Thank you for correcting an error on my part. Did you use a diagnostic tool or was it the result of eyeballing the error ? Attached is a somewhat improved corrected expression.  Comment By: Barton Willis (willisbl) Date: 20110801 20:38 Message: The (2(b*x ... makes this an invalid expression solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2(b*x < Maybe you could post a corrected expression  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377380&group_id=4933 
From: SourceForge.net <noreply@so...>  20110819 17:42:10

Bugs item #2974948, was opened at 20100323 00:52 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2974948&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Floating point Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: rectform loses some accuracy Initial Comment: Consider: rectform(1e300/(1e300+3/2*%i)) > 1.0 But float(rectform(10^300/(10^300+3/2*%i))) > 1.0  1.5e300 * %i The inaccuracy is caused by Maxima computing the rectform in two parts. Maxima computes 1/(1e300+3/2*%i) and then multiplies by 1e300. Maxima correctly computes 1/(1e300+3/2*%i) as 1e300. But doing it in two steps causes the very tiny imaginary part to be lost.  >Comment By: Dieter Kaiser (crategus) Date: 20110819 19:42 Message: Sorry, but today I have checked the example again. Now I get the described results with SBCL: Maxima version: 5.25post Maxima build date: 19:8 8/14/2011 Host type: i686pclinuxgnu Lisp implementation type: SBCL Lisp implementation version: 1.0.45 (%i5) rectform(1e300/(1e300+3/2*%i)); (%o5) 1.0 (%i9) float(rectform(10^300/(10^300+3/2*%i))); (%o9) 1.01.5e300*%i Perhaps, I had some own code in my sandbox. I think I have tried to use the bigfloatpackage. Dieter Kaiser  Comment By: Raymond Toy (rtoy) Date: 20110819 06:54 Message: That's odd. I tried it with sbcl, and it had no problems with overflow with the given examples. More investigation needed.  Comment By: Dieter Kaiser (crategus) Date: 20110804 23:05 Message: For the record: With SBCL and Maxima 5.25post I get for the example an overflow error: (%i1) rectform(1e300/(1e300+3/2*%i)); Maxima encountered a Lisp error: arithmetic error FLOATINGPOINTOVERFLOW signalled Automatically continuing. To reenable the Lisp debugger set *debuggerhook* to nil. I have tried the biggest possible exponent: (%i31) float(rectform(1e154/(1e154+3/2*%i))); (%o31) 1.0  1.5e154 %i (%i32) float(rectform(10^154/(10^154+3/2*%i))); (%o32) 1.0  1.5e154 %i Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2974948&group_id=4933 
From: SourceForge.net <noreply@so...>  20110819 04:54:32

Bugs item #2974948, was opened at 20100322 19:52 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2974948&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Floating point Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: rectform loses some accuracy Initial Comment: Consider: rectform(1e300/(1e300+3/2*%i)) > 1.0 But float(rectform(10^300/(10^300+3/2*%i))) > 1.0  1.5e300 * %i The inaccuracy is caused by Maxima computing the rectform in two parts. Maxima computes 1/(1e300+3/2*%i) and then multiplies by 1e300. Maxima correctly computes 1/(1e300+3/2*%i) as 1e300. But doing it in two steps causes the very tiny imaginary part to be lost.  >Comment By: Raymond Toy (rtoy) Date: 20110819 00:54 Message: That's odd. I tried it with sbcl, and it had no problems with overflow with the given examples. More investigation needed.  Comment By: Dieter Kaiser (crategus) Date: 20110804 17:05 Message: For the record: With SBCL and Maxima 5.25post I get for the example an overflow error: (%i1) rectform(1e300/(1e300+3/2*%i)); Maxima encountered a Lisp error: arithmetic error FLOATINGPOINTOVERFLOW signalled Automatically continuing. To reenable the Lisp debugger set *debuggerhook* to nil. I have tried the biggest possible exponent: (%i31) float(rectform(1e154/(1e154+3/2*%i))); (%o31) 1.0  1.5e154 %i (%i32) float(rectform(10^154/(10^154+3/2*%i))); (%o32) 1.0  1.5e154 %i Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2974948&group_id=4933 
From: SourceForge.net <noreply@so...>  20110818 02:49:10

Bugs item #3392909, was opened at 20110817 04:03 Message generated for change (Settings changed) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3392909&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Mikael Samsøe Sørensen (mrverify) >Assigned to: Barton Willis (willisbl) Summary: to_poly_solve missing solution Initial Comment: to_poly_solve(log(x)/log(2)=x1,x); %union([x=1]) But x=2 is also a solution Maxima version 5.25.0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3392909&group_id=4933 
From: SourceForge.net <noreply@so...>  20110817 09:03:39

Bugs item #3392909, was opened at 20110817 11:03 Message generated for change (Tracker Item Submitted) made by mrverify You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3392909&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Share Libraries Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Mikael Samsøe Sørensen (mrverify) Assigned to: Nobody/Anonymous (nobody) Summary: to_poly_solve missing solution Initial Comment: to_poly_solve(log(x)/log(2)=x1,x); %union([x=1]) But x=2 is also a solution Maxima version 5.25.0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3392909&group_id=4933 
From: SourceForge.net <noreply@so...>  20110815 21:04:26

Bugs item #3377380, was opened at 20110725 14:04 Message generated for change (Settings changed) made by dloksnel You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377380&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Solving equations Group: None Status: Open >Resolution: Wont Fix >Priority: 6 Private: No Submitted By: Art Lenskold (dloksnel) Assigned to: Nobody/Anonymous (nobody) Summary: 7 nested levels Initial Comment: solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2(b*xt*sqrt(1x^2))*(bsqrt(1x^2))2*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)))^1*(a+(11/(4*((tx)^2+(bsqrt(1x^2))^2)^2)*(4*(b*xt*sqrt(1x^2))^2*(bsqrt(1x^2))^28*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))^(1/2)*(b*xt*sqrt(1x^2))*(bsqrt(1x^2))+4*(x^2*(b^2t^2)*(x^22+2*t*xt^2)+t^2*(x^2+b^21)(2*t*xx^2)*(1+b^2+t^22*t*x)2*b*sqrt(1x^2)*((tx)^2+t*x*(x^2+2*t*xt^22)))))^(1/2)))^2 ], [x] ); Maxima version: 5.24.0 Maxima build date: 20:39 4/5/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  Comment By: Art Lenskold (dloksnel) Date: 20110802 19:18 Message: Thank you for correcting an error on my part. Did you use a diagnostic tool or was it the result of eyeballing the error ? Attached is a somewhat improved corrected expression.  Comment By: Barton Willis (willisbl) Date: 20110801 20:38 Message: The (2(b*x ... makes this an invalid expression solve([0=x^21+(a(tx)*(t+1/(2*((tx)^2+(bsqrt(1x^2))^2))*(2(b*x < Maybe you could post a corrected expression  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377380&group_id=4933 
From: SourceForge.net <noreply@so...>  20110814 19:48:36

Bugs item #3377114, was opened at 20110725 14:20 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377114&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None Status: Open Resolution: None Priority: 2 Private: No Submitted By: Dieter Kaiser (crategus) Assigned to: Nobody/Anonymous (nobody) Summary: Linear display of long symbols and strings Initial Comment: In linear display expressions with long symbols or strings are not formated as expected. This is an example for 2ddisplay and for linear display. The option variable linel is set to 10: (%i1) display2d:true$ (%i2) a+b+c+ddddddd+eeeeeee+ffffffffffffff; (%o2) fff\ fffffffff\ ff + eeee\ eee + ddd\ dddd + c + b + a Again for linear display. The spaces in the output are replaced by points to show the format in this posting: (%i3) display2d:false$ (%i4) a+b+c+ddddddd+eeeeeee+ffffffffffffff; (%o4) ffffffffffffff .......+eeeeeee .......+ddddddd .......+c .......+b .......+a In linear display there are two problems. First, long symbols and strings are not broken up. This is a problem in the function mprint. Second, in linear display the available space is not used fully. This problem is caused in the function lineardisplay. This function prints out the output label directly. Therefore the format of the output label is not under control of the function mprint. The second problem is solved by providing a function msizelabel, which does the formating of the output label under control of the functions msize and mprint. This is an example with a new msizelabel function: (%i1) a+b+c+ddddddd+eeeeeee+ffffffffffffff; (%o1) ffffffffffffff +eeeeeee +ddddddd +c+b+a Still, the first problem is present. Long symbols are not broken up. Dieter Kaiser  >Comment By: Dieter Kaiser (crategus) Date: 20110814 21:48 Message: In linear display we have no functions for formating mlable and mtextexpressions. The formating is done in the function lineardispla. Therefore, the output of such expressions in lineardispla is not under the control of the function mprint which is responsible for the output of a formated expression. In 2ddisplay this problem does not exist. The following two functions implement the formating of mlable and mtextexpression for linear display: ;;;  ;; Formating a mlable expression (defprop mlable msizemlabel grind) (defun msizemlabel (x l r) (declare (special *displaylabelsp*)) (cond (*displaylabelsp* (setq l (msize (cadr x) (list #\( ) (list #\) #\ ) nil nil) r (msize (caddr x) nil r 'mparen 'mparen)) (cons (+ (car l) (car r)) (cons l (cons r nil)))) (t (msize (caddr x) l r 'mparen 'mparen)))) ;;;  ;; Formating a mtext expression (defprop mtext msizemtext grind) (defun msizemtext (x l r) (setq x (cdr x)) (if (null x) (msz nil l r) (do ((nl) (w 0)) ((null (cdr x)) (setq nl (cons (if (atom (car x)) (msz (makestring (car x)) l r) (msize (car x) l r lop rop)) nl)) (cons (+ w (caar nl)) (nreverse nl))) (setq nl (cons (if (atom (car x)) (msz (makestring (car x)) l r) (msize (car x) l nil lop rop)) nl) w (+ w (caar nl)) x (cdr x) l nil)))) ;;;  With this definitions the function lineardispla reduces to the following code: (defun lineardispla (form) (freshline *standardoutput*) (mgrind form *standardoutput*) (terpri)) Furthermore, the global variable fortranp is no longer necessary and can be cut out completely. It is set in lineardispla to signal the function msizeatom that a mtextexpression has to be formated. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3377114&group_id=4933 
From: SourceForge.net <noreply@so...>  20110814 18:01:17

Bugs item #3390551, was opened at 20110812 14:31 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3390551&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None >Status: Closed >Resolution: Fixed Priority: 3 Private: No Submitted By: Dieter Kaiser (crategus) Assigned to: Nobody/Anonymous (nobody) Summary: Error break in rtest2 with display_all=true Initial Comment: In 2Ddisplay and with display_all=true we have an error break in the test file rtest2.mac for problem 55. (%i1) run_testsuite(display_all=true, tests=[rtest2]); Running tests in rtest2: [...] ********************** Problem 55 *************** Input: Caused an error break: rtest2 Error summary: Error found in rtest2, problem: (error break) 0 tests failed out of 0 total tests. This is the example 55 which fails. The error break occurs only if the flag display_all is set to true. /* should trigger an error */ errcatch (matrix ([1], [1, 2])); []; The reason is, that in a test mode an unsimplified expression is displayed. But in this case it is a bad formed unsimplified expression. The matrix is not well formed. The wrong matrix causes an error break in the function dimension$matrix. The problem has been introduced in displa.lisp with revision 15.03.2011. The check for an unsimplified matrix has been cut out by me. This has been done, because the linear display of an unsimplified matrix has been considered as a bug in ID: 635708 Bad display of unsimplified expr. A possible correction is to check more careful, if we have a valid matrix in dim$matrix. The following code allows the display of bad formed matrices without an error break. The code has been extended with a check for a equal number of columns. (defun dim$matrix (form result) (prog (dmstr rstr cstr consp cols) (setq cols (if ($listp (cadr form)) (length (cadr form)) 0)) (if (or (null (cdr form)) (memalike '((mlist simp)) (cdr form)) ;; Check if the matrix has lists as rows with a equal number of ;; columns. (dolist (row (cdr form)) (if (or (not ($listp row)) (not (eql cols (length row)))) (return t)))) ;; The matrix is not well formed. Display the matrix in linear mode. (return (dimensionfunction form result))) [...] Dieter Kaiser  >Comment By: Dieter Kaiser (crategus) Date: 20110814 20:01 Message: Fixed in displa.lisp revision 14.08.2011. Closing this bug report as fixed. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3390551&group_id=4933 
From: SourceForge.net <noreply@so...>  20110813 08:50:10

Bugs item #3133029, was opened at 20101209 15:48 Message generated for change (Comment added) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3133029&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Complex Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: Incorrect result of residue Initial Comment: From https://bugs.launchpad.net/ubuntu/+source/maxima/+bug/323221: kubuntu intrepid, maxima 5.13.0 Binary package hint: maxima (%i1) residue(%e**(1/z),z,0); (%o1) 0 but according to the laurent series: %e**(1/z) = 1 + 1/z + 1/2z^2 + ... the correct result is 1.  I can confirm this bug in Maxima 5.22.1. Laurent serie with `powerseries' is correct: (%i2) powerseries(%e^(1/z), z, 0); inf ==== \ 1 (%o2) >  / i2 ==== i2! z i2 = 0 but with `taylor': (%i3) taylor(%e^(1/z), z, 0, 8); 1 (%o3)/T/  + . . .  1/z %e  Comment By: https://www.google.com/accounts () Date: 20110813 10:50 Message: This bug is still present in Maxima 5.25.0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3133029&group_id=4933 
From: SourceForge.net <noreply@so...>  20110812 14:42:43

Bugs item #2998214, was opened at 20100507 17:18 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2998214&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None >Status: Closed Resolution: Works For Me Priority: 5 Private: No Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Cell can't end with a comment Initial Comment: When the last line in a cell is comment, maxima comlains about premature termination of input. To reproduce: Just enter a comment only in a input cell. /*Comment*/ and evaluate. Why do you need that? If you want to do somehting like this: a:99 /*assign a the magic value*/  >Comment By: Dieter Kaiser (crategus) Date: 20110812 16:42 Message: Setting the status to closed. Dieter Kaiser  Comment By: Jose SanchezMarin (jsm) Date: 20110809 19:55 Message: Hi It's not so elegant, but I always ends a comment at the end of a cell with 1$: /* comment*/1$ Of course, the 1 is optative: whatever acceptable character is good. It is the simplest way of not ending the cell with a comment. Your example should look like this: a:99 /*assign a the magic value*/1$  Comment By: Dieter Kaiser (crategus) Date: 20110804 22:50 Message: As written in one of the posts, it is save to use the syntax a:99 /* Comment */ ; in wxMaxima. The following syntax a:99; /* Comment */ might give an error in wxMaxima, because wxMaxima automatically appends an ;char and a comment followed by a ;char is not a valid command. Furthermore, I do not see a problem with a comment at the end of a file with SBCL. I have tried the Maxima versions 5.9, 5.10, 5.11, 5.12, 5.13, 5.14, 5.15, 5.16, 5.17, 5.18, 5.19, 5.20, and 5.25post and had no problems loading the test file of the last posting. I have tried several variations of the file and had a look at the code to search for a potential problem. Perhaps, we can close this bug report and open a new bug report, if we have a report with a known Lisp and Maxima version. Setting the status to pending and the resolution to "works for me". Dieter Kaiser  Comment By: Robert Dodier (robert_dodier) Date: 20100603 17:54 Message: This is actually a bug in the Maxima parser if I am not mistaken; it is not specific to wxMaxima so far as I know. At one time (maybe it is still true) the parser would barf on a batch script if it ended with a comment. I just tried that with Maxima 5.20.0 + GCL + Windows but the script was loaded successfully. Maybe the bug has been fixed? Need to test on other platforms and Lisp implementations. To the original poster: what does build_info(); report? Here is my test script: x:1; /* comment */  Comment By: Nobody/Anonymous (nobody) Date: 20100527 08:55 Message: Got it. A comment can't be the only text in a field. /*comment*/a:99 and a:99/*comment*/ both work ok. wxMaxima appears to append a semicolon to the line if it is absent. /*comment*/ gets a semicolon appended, then an error message "incorrect syntax: Premature termination of input at ;. /*comment*/; ^" With the caret indicating the closing slash of the comment. Using a $ terminator gets a similar message, no semicolon is appended, and the caret indicates the closing slash as before.  Comment By: Peter Cusack (pogo1) Date: 20100527 03:53 Message: I note that wmMaxima version 0.8.5 helpfully adds a semicolon and that fixes the problem sometimes. I'm still trying to make it reproducibly fail. If the line contains a $ terminator before the comment ( like this B0:i*Mu0/(2*a)/*Field at centre of single turn*/$) then adding the semicolon at the end makes Maxima hang.  Comment By: l_butler () Date: 20100520 10:16 Message: Hi, I think the problem is that Maxima input must be terminated with a semicolon: a:99; /* comment */ a:99 /* comment */; Both are valid inputs that assign 99 to a.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2998214&group_id=4933 
From: SourceForge.net <noreply@so...>  20110812 12:31:23

Bugs item #3390551, was opened at 20110812 14:31 Message generated for change (Tracker Item Submitted) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3390551&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None Status: Open Resolution: None Priority: 3 Private: No Submitted By: Dieter Kaiser (crategus) Assigned to: Nobody/Anonymous (nobody) Summary: Error break in rtest2 with display_all=true Initial Comment: In 2Ddisplay and with display_all=true we have an error break in the test file rtest2.mac for problem 55. (%i1) run_testsuite(display_all=true, tests=[rtest2]); Running tests in rtest2: [...] ********************** Problem 55 *************** Input: Caused an error break: rtest2 Error summary: Error found in rtest2, problem: (error break) 0 tests failed out of 0 total tests. This is the example 55 which fails. The error break occurs only if the flag display_all is set to true. /* should trigger an error */ errcatch (matrix ([1], [1, 2])); []; The reason is, that in a test mode an unsimplified expression is displayed. But in this case it is a bad formed unsimplified expression. The matrix is not well formed. The wrong matrix causes an error break in the function dimension$matrix. The problem has been introduced in displa.lisp with revision 15.03.2011. The check for an unsimplified matrix has been cut out by me. This has been done, because the linear display of an unsimplified matrix has been considered as a bug in ID: 635708 Bad display of unsimplified expr. A possible correction is to check more careful, if we have a valid matrix in dim$matrix. The following code allows the display of bad formed matrices without an error break. The code has been extended with a check for a equal number of columns. (defun dim$matrix (form result) (prog (dmstr rstr cstr consp cols) (setq cols (if ($listp (cadr form)) (length (cadr form)) 0)) (if (or (null (cdr form)) (memalike '((mlist simp)) (cdr form)) ;; Check if the matrix has lists as rows with a equal number of ;; columns. (dolist (row (cdr form)) (if (or (not ($listp row)) (not (eql cols (length row)))) (return t)))) ;; The matrix is not well formed. Display the matrix in linear mode. (return (dimensionfunction form result))) [...] Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3390551&group_id=4933 
From: SourceForge.net <noreply@so...>  20110811 23:51:50

Bugs item #3387042, was opened at 20110806 00:04 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None >Status: Pending >Resolution: Works For Me Priority: 5 Private: No Submitted By: Petr Kartsev (kpf) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x exp(x^2)) : strange jump Initial Comment: (%i1) integrate( t*exp((ta)^2), t, inf, x); (%o1) (gamma_incomplete(1/2,x^22*a*x+a^2)*a*x)/(2*abs(xa))+(gamma_incomplete(1/2,x^22*a*x+a^2)*a^2)/(2*abs(xa))gamma_incomplete(1,x^22*a*x+a^2)/2 (1st part on screenshot) why is it so complex? when plotting this integral, it has strange jump at point x=1 (and even becomes negative) while the integrand is smooth and mostly always positive) (%i2) plot2d([x*exp((x1)^2), integrate(t*exp((t1)^2), t, inf,x)], [x,5,5])$ (2nd part on screenshot) However, if I simply add and subtract unity, the actual result is recovered: (%i3) integrate( (t1)*exp((t1)^2), t, inf, x)+ integrate( exp((t1)^2), t, inf, x); (%o3) (sqrt(%pi)*erf(x1))/2%e^(x^2+2*x1)/2+sqrt(%pi)/2 (3rd part on screenshot) Looks like additional signum or abs takes place somewhere? Sorry for not giving solution, I am not that skillful in lisp... Thanks for good system, anyway! :)  Maxima version: 5.24.0 Maxima build date: 10:35 8/1/2011 Host type: x86_64pclinuxgnu Lisp implementation type: SBCL Lisp implementation version: 1.0.19gentoo  PS. Just for information, if this can help, run_testsuite() gives 1 error. May be just my installation is wrong? Running tests in rtest16: ********************** Problem 386 *************** Input: closeto(zeta(%i + 3)  (1.10721440843141  .1482908671781754 %i), 1.e15) Result: 3.3157171357748244e9 This differed from the expected result: true   >Comment By: Dieter Kaiser (crategus) Date: 20110812 01:51 Message: One way to avoid questions is to give Maxima facts which are stored in the assume database, e.g. (%i1) assume(x>1); (%o1) [x > 1] (%i2) res:integrate( t*exp((t1)^2), t,minf,x); (%o2) gamma_incomplete(1/2,x^22*x+1)*x/(2*(1x)) gamma_incomplete(1/2,x^22*x+1)/(2*(1x)) gamma_incomplete(1,x^22*x+1)/2+sqrt(%pi) The result can be transformed to a representation in terms of the erfc function with the option variable gamma_expand: (%i3) res, gamma_expand:true; (%o3) sqrt(%pi)*x*erfc(sqrt(x^22*x+1))/(2*(1x)) sqrt(%pi)*erfc(sqrt(x^22*x+1))/(2*(1x))%e^(x^2+2*x1)/2+sqrt(%pi) In addition it is possible to transform to the erf function when setting the option variable erf_representation: (%i4) res,gamma_expand:true, erf_representation:erf; (%o4) sqrt(%pi)*x*(1erf(sqrt(x^22*x+1)))/(2*(1x)) sqrt(%pi)*(1erf(sqrt(x^22*x+1)))/(2*(1x))%e^(x^2+2*x1)/2 +sqrt(%pi) At this point I think we do not have a bug. I suggest to close this bug report as "works for me". Setting the status to pending. Dieter Kaiser  Comment By: Petr Kartsev (kpf) Date: 20110812 00:30 Message: willisbl, thanks for answer. I upgraded to version 5.25 (appeared in Gentoo Linux repository ) and now the behavior is exactly as you posted! However, my question transforms to the following, may be I should file it as dedicated bug report: Can you remove the question "positive/negative/zero" in this situation, as we can see the result is the same (already known erf+exp) for all cases ? This can be possible, in my opinion, since this question is already removed by simply shifting the variable: (%i6) e1(x): integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?" pos; (%o6) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2+sqrt(%pi) but (%i8) e1: integrate( (t+1)*exp(t^2), t,minf,x1); (%o8) (%e^(x^21)*(sqrt(%pi)*%e^(x^2+1)*erf(x1)%e^(2*x)))/2+sqrt(%pi)/2 (My reason is that in my work I differentiate this function and then do some numerical actions for programmed solution with external parameters. Such automatic analysis would get in trouble if it needs to answer positive/negative) Maybe gamma_incomplete is too universal to let simple erf appear in the formula, but erf is sometimes better, so can some special case be implemented? Sorry if I see something wrong. Thanks for attention and for this good tool! Sincerely, Petr F. Kartsev  Comment By: Barton Willis (willisbl) Date: 20110806 03:25 Message: First, you should check your initialization filemaybe you have set some option variables to non default values ; second try this: (%i23) e1 : integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?"neg; (%o23) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))+gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2 (%i24) e2 : integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?"pos; (%o24) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2+sqrt(%pi) (%i25) plot2d(if x < 1 then ''e1 else ''e2,[x,10,10]); The graph appears to be continuous.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 
From: SourceForge.net <noreply@so...>  20110811 22:30:46

Bugs item #3387042, was opened at 20110806 02:04 Message generated for change (Comment added) made by kpf You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Petr Kartsev (kpf) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x exp(x^2)) : strange jump Initial Comment: (%i1) integrate( t*exp((ta)^2), t, inf, x); (%o1) (gamma_incomplete(1/2,x^22*a*x+a^2)*a*x)/(2*abs(xa))+(gamma_incomplete(1/2,x^22*a*x+a^2)*a^2)/(2*abs(xa))gamma_incomplete(1,x^22*a*x+a^2)/2 (1st part on screenshot) why is it so complex? when plotting this integral, it has strange jump at point x=1 (and even becomes negative) while the integrand is smooth and mostly always positive) (%i2) plot2d([x*exp((x1)^2), integrate(t*exp((t1)^2), t, inf,x)], [x,5,5])$ (2nd part on screenshot) However, if I simply add and subtract unity, the actual result is recovered: (%i3) integrate( (t1)*exp((t1)^2), t, inf, x)+ integrate( exp((t1)^2), t, inf, x); (%o3) (sqrt(%pi)*erf(x1))/2%e^(x^2+2*x1)/2+sqrt(%pi)/2 (3rd part on screenshot) Looks like additional signum or abs takes place somewhere? Sorry for not giving solution, I am not that skillful in lisp... Thanks for good system, anyway! :)  Maxima version: 5.24.0 Maxima build date: 10:35 8/1/2011 Host type: x86_64pclinuxgnu Lisp implementation type: SBCL Lisp implementation version: 1.0.19gentoo  PS. Just for information, if this can help, run_testsuite() gives 1 error. May be just my installation is wrong? Running tests in rtest16: ********************** Problem 386 *************** Input: closeto(zeta(%i + 3)  (1.10721440843141  .1482908671781754 %i), 1.e15) Result: 3.3157171357748244e9 This differed from the expected result: true   >Comment By: Petr Kartsev (kpf) Date: 20110812 02:30 Message: willisbl, thanks for answer. I upgraded to version 5.25 (appeared in Gentoo Linux repository ) and now the behavior is exactly as you posted! However, my question transforms to the following, may be I should file it as dedicated bug report: Can you remove the question "positive/negative/zero" in this situation, as we can see the result is the same (already known erf+exp) for all cases ? This can be possible, in my opinion, since this question is already removed by simply shifting the variable: (%i6) e1(x): integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?" pos; (%o6) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2+sqrt(%pi) but (%i8) e1: integrate( (t+1)*exp(t^2), t,minf,x1); (%o8) (%e^(x^21)*(sqrt(%pi)*%e^(x^2+1)*erf(x1)%e^(2*x)))/2+sqrt(%pi)/2 (My reason is that in my work I differentiate this function and then do some numerical actions for programmed solution with external parameters. Such automatic analysis would get in trouble if it needs to answer positive/negative) Maybe gamma_incomplete is too universal to let simple erf appear in the formula, but erf is sometimes better, so can some special case be implemented? Sorry if I see something wrong. Thanks for attention and for this good tool! Sincerely, Petr F. Kartsev  Comment By: Barton Willis (willisbl) Date: 20110806 05:25 Message: First, you should check your initialization filemaybe you have set some option variables to non default values ; second try this: (%i23) e1 : integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?"neg; (%o23) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))+gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2 (%i24) e2 : integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?"pos; (%o24) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2+sqrt(%pi) (%i25) plot2d(if x < 1 then ''e1 else ''e2,[x,10,10]); The graph appears to be continuous.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 
From: SourceForge.net <noreply@so...>  20110811 19:29:38

Bugs item #3389830, was opened at 20110811 00:19 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3389830&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None >Status: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: Dieter Kaiser (crategus) Assigned to: Nobody/Anonymous (nobody) Summary: Error break in rtest15 with linear display Initial Comment: I have run the testfile rtest15.mac in linear display and with display_all=true and I have got an error break. (%i1) display2d:false; (%o1) false (%i2) run_testsuite(display_all=true, tests=[rtest15]); Running tests in rtest15: [...] ********************** Problem 203 *************** Input: Caused an error break: rtest15 Error summary: Error found in rtest15, problem: (error break) 0 tests failed out of 0 total tests. This is a short example to show the problem: (%i1) display2d:false; (%o1) false (%i2) postfix("f"); (%o2) "f" (%i3) "f"(x) := x; (%o3) Maxima encountered a Lisp error: The value #\x is not of type LIST. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. Maxima can not display the definition in (%i3) for a postfixoperator in linear display. It works in 2ddisplay. It is an old bug. I have observed it in all Maxima versions since 5.9. Dieter Kaiser  >Comment By: Dieter Kaiser (crategus) Date: 20110811 21:29 Message: Fixed in grind.lisp revision 11.08.2011. The function mszmdef has been extended to handle prefix and postfix operators. Closing this bug report as fixed. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3389830&group_id=4933 
From: SourceForge.net <noreply@so...>  20110810 22:19:02

Bugs item #3389830, was opened at 20110811 00:19 Message generated for change (Tracker Item Submitted) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3389830&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Dieter Kaiser (crategus) Assigned to: Nobody/Anonymous (nobody) Summary: Error break in rtest15 with linear display Initial Comment: I have run the testfile rtest15.mac in linear display and with display_all=true and I have got an error break. (%i1) display2d:false; (%o1) false (%i2) run_testsuite(display_all=true, tests=[rtest15]); Running tests in rtest15: [...] ********************** Problem 203 *************** Input: Caused an error break: rtest15 Error summary: Error found in rtest15, problem: (error break) 0 tests failed out of 0 total tests. This is a short example to show the problem: (%i1) display2d:false; (%o1) false (%i2) postfix("f"); (%o2) "f" (%i3) "f"(x) := x; (%o3) Maxima encountered a Lisp error: The value #\x is not of type LIST. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. Maxima can not display the definition in (%i3) for a postfixoperator in linear display. It works in 2ddisplay. It is an old bug. I have observed it in all Maxima versions since 5.9. Dieter Kaiser  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3389830&group_id=4933 
From: Mike Welland <mike@mi...>  20110810 14:54:31

Hello, When I run maxima from the terminal, I just get 'Segmentation fault' and that's all. About one month ago it was working fine. I am running Archlinux 64bit. Thanks, Mike 
From: SourceForge.net <noreply@so...>  20110809 17:55:43

Bugs item #2998214, was opened at 20100507 17:18 Message generated for change (Comment added) made by jsm You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2998214&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None Status: Pending Resolution: Works For Me Priority: 5 Private: No Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Cell can't end with a comment Initial Comment: When the last line in a cell is comment, maxima comlains about premature termination of input. To reproduce: Just enter a comment only in a input cell. /*Comment*/ and evaluate. Why do you need that? If you want to do somehting like this: a:99 /*assign a the magic value*/  Comment By: Jose SanchezMarin (jsm) Date: 20110809 19:55 Message: Hi It's not so elegant, but I always ends a comment at the end of a cell with 1$: /* comment*/1$ Of course, the 1 is optative: whatever acceptable character is good. It is the simplest way of not ending the cell with a comment. Your example should look like this: a:99 /*assign a the magic value*/1$  Comment By: Dieter Kaiser (crategus) Date: 20110804 22:50 Message: As written in one of the posts, it is save to use the syntax a:99 /* Comment */ ; in wxMaxima. The following syntax a:99; /* Comment */ might give an error in wxMaxima, because wxMaxima automatically appends an ;char and a comment followed by a ;char is not a valid command. Furthermore, I do not see a problem with a comment at the end of a file with SBCL. I have tried the Maxima versions 5.9, 5.10, 5.11, 5.12, 5.13, 5.14, 5.15, 5.16, 5.17, 5.18, 5.19, 5.20, and 5.25post and had no problems loading the test file of the last posting. I have tried several variations of the file and had a look at the code to search for a potential problem. Perhaps, we can close this bug report and open a new bug report, if we have a report with a known Lisp and Maxima version. Setting the status to pending and the resolution to "works for me". Dieter Kaiser  Comment By: Robert Dodier (robert_dodier) Date: 20100603 17:54 Message: This is actually a bug in the Maxima parser if I am not mistaken; it is not specific to wxMaxima so far as I know. At one time (maybe it is still true) the parser would barf on a batch script if it ended with a comment. I just tried that with Maxima 5.20.0 + GCL + Windows but the script was loaded successfully. Maybe the bug has been fixed? Need to test on other platforms and Lisp implementations. To the original poster: what does build_info(); report? Here is my test script: x:1; /* comment */  Comment By: Nobody/Anonymous (nobody) Date: 20100527 08:55 Message: Got it. A comment can't be the only text in a field. /*comment*/a:99 and a:99/*comment*/ both work ok. wxMaxima appears to append a semicolon to the line if it is absent. /*comment*/ gets a semicolon appended, then an error message "incorrect syntax: Premature termination of input at ;. /*comment*/; ^" With the caret indicating the closing slash of the comment. Using a $ terminator gets a similar message, no semicolon is appended, and the caret indicates the closing slash as before.  Comment By: Peter Cusack (pogo1) Date: 20100527 03:53 Message: I note that wmMaxima version 0.8.5 helpfully adds a semicolon and that fixes the problem sometimes. I'm still trying to make it reproducibly fail. If the line contains a $ terminator before the comment ( like this B0:i*Mu0/(2*a)/*Field at centre of single turn*/$) then adding the semicolon at the end makes Maxima hang.  Comment By: l_butler () Date: 20100520 10:16 Message: Hi, I think the problem is that Maxima input must be terminated with a semicolon: a:99; /* comment */ a:99 /* comment */; Both are valid inputs that assign 99 to a.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2998214&group_id=4933 
From: SourceForge.net <noreply@so...>  20110809 03:04:10

Bugs item #3388801, was opened at 20110809 11:04 Message generated for change (Tracker Item Submitted) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3388801&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: None Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: integrate error Initial Comment: I use maxima on windows xp, Maxima 5.25.0 http://maxima.sourceforge.net using Lisp Clozure Common Lisp Version 1.7r14925M (WindowsX8632) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) f1:sqrt(107669.2969(x474.9481828977571)^2) + 621.87$ (%i2) f2:sqrt(107669.2969(x474.9481828977571)^2) $ (%i3) A1 : float (integrate(2*%pi * f1 * sqrt(1+diff(f1,x)^2) , x ,181.4663110563288,474.9481828977571))$ (%i4) A2 : float (integrate(2*%pi * f2 * sqrt(1+diff(f2,x)^2) , x ,181.4663110563288,474.9481828977571))$ (%i5) [A1,A2]; (%o5)[605072.0511807054, 605072.0511807053] The A1 and A2 should not be the same 。 And I use maxima 5.13 on Debian 5.0 Maxima 5.13.0 http://maxima.sourceforge.net Using Lisp GNU Common Lisp (GCL) GCL 2.6.7 (aka GCL) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. This is a development version of Maxima. The function bug_report() provides bug reporting information. (%i1) f1:sqrt(107669.2969(x474.9481828977571)^2) + 621.87$ (%i2) f2:sqrt(107669.2969(x474.9481828977571)^2) $ (%i3) A1 : float (integrate(2*%pi * f1 * sqrt(1+diff(f1,x)^2) , x ,181.4663110563288,474.9481828977571))$ (%i4) A2 : float (integrate(2*%pi * f2 * sqrt(1+diff(f2,x)^2) , x ,181.4663110563288,474.9481828977571))$ (%i5) [A1,A2]; (%o5)[2024501.89874211, 605072.0511807054] I check A1 and A2 are right the Maxima 5.21.0 by Debian 6.0 ， the calculation was error. So, why 5.25.0's or 5.21.0's calculation is error?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3388801&group_id=4933 
From: SourceForge.net <noreply@so...>  20110807 09:56:59

Bugs item #3387673, was opened at 20110807 04:56 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387673&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: None Group: None Status: Open Resolution: None Priority: 4 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: newdet limitation Initial Comment: The function newdet sometimes gives a Lisp error: (%i1) m : genmatrix(lambda([a,b], random(10)),30,30)$ (%i2) newdet(m); Maxima encountered a Lisp error: value 310235040 is not of the expected type (UNSIGNEDBYTE 24).  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387673&group_id=4933 
From: SourceForge.net <noreply@so...>  20110806 01:25:06

Bugs item #3387042, was opened at 20110805 17:04 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Petr Kartsev (kpf) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x exp(x^2)) : strange jump Initial Comment: (%i1) integrate( t*exp((ta)^2), t, inf, x); (%o1) (gamma_incomplete(1/2,x^22*a*x+a^2)*a*x)/(2*abs(xa))+(gamma_incomplete(1/2,x^22*a*x+a^2)*a^2)/(2*abs(xa))gamma_incomplete(1,x^22*a*x+a^2)/2 (1st part on screenshot) why is it so complex? when plotting this integral, it has strange jump at point x=1 (and even becomes negative) while the integrand is smooth and mostly always positive) (%i2) plot2d([x*exp((x1)^2), integrate(t*exp((t1)^2), t, inf,x)], [x,5,5])$ (2nd part on screenshot) However, if I simply add and subtract unity, the actual result is recovered: (%i3) integrate( (t1)*exp((t1)^2), t, inf, x)+ integrate( exp((t1)^2), t, inf, x); (%o3) (sqrt(%pi)*erf(x1))/2%e^(x^2+2*x1)/2+sqrt(%pi)/2 (3rd part on screenshot) Looks like additional signum or abs takes place somewhere? Sorry for not giving solution, I am not that skillful in lisp... Thanks for good system, anyway! :)  Maxima version: 5.24.0 Maxima build date: 10:35 8/1/2011 Host type: x86_64pclinuxgnu Lisp implementation type: SBCL Lisp implementation version: 1.0.19gentoo  PS. Just for information, if this can help, run_testsuite() gives 1 error. May be just my installation is wrong? Running tests in rtest16: ********************** Problem 386 *************** Input: closeto(zeta(%i + 3)  (1.10721440843141  .1482908671781754 %i), 1.e15) Result: 3.3157171357748244e9 This differed from the expected result: true   >Comment By: Barton Willis (willisbl) Date: 20110805 20:25 Message: First, you should check your initialization filemaybe you have set some option variables to non default values ; second try this: (%i23) e1 : integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?"neg; (%o23) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))+gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2 (%i24) e2 : integrate( t*exp((t1)^2), t,minf,x); "Is "x1" positive, negative, or zero?"pos; (%o24) (gamma_incomplete(1/2,x^22*x+1)*x)/(2*(1x))gamma_incomplete(1/2,x^22*x+1)/(2*(1x))gamma_incomplete(1,x^22*x+1)/2+sqrt(%pi) (%i25) plot2d(if x < 1 then ''e1 else ''e2,[x,10,10]); The graph appears to be continuous.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 
From: SourceForge.net <noreply@so...>  20110805 23:59:58

Bugs item #3387081, was opened at 20110805 18:59 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387081&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Limit Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: question from limit involving a float Initial Comment: The limit code has far too much code that fills the fact data base with comparisons to 1.0*10^8 or 1.0*10^8. One of these silly facts leaks out: (%i16) integrate( t*exp((ta)^2), t)$ Silly question: (%i17) limit(%,t,minf); "Is "t+a+1.0*10^8" positive or negative?"neg; (%o17) 0 (%i18) build_info(); Maxima version: 5.25.0 Maxima build date: 12:0 8/2/2011 Host type: i686pcmingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7r14925M (WindowsX8632) (%o18)  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387081&group_id=4933 
From: SourceForge.net <noreply@so...>  20110805 22:04:04

Bugs item #3387042, was opened at 20110806 02:04 Message generated for change (Tracker Item Submitted) made by kpf You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core  Integration Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Petr Kartsev (kpf) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x exp(x^2)) : strange jump Initial Comment: (%i1) integrate( t*exp((ta)^2), t, inf, x); (%o1) (gamma_incomplete(1/2,x^22*a*x+a^2)*a*x)/(2*abs(xa))+(gamma_incomplete(1/2,x^22*a*x+a^2)*a^2)/(2*abs(xa))gamma_incomplete(1,x^22*a*x+a^2)/2 (1st part on screenshot) why is it so complex? when plotting this integral, it has strange jump at point x=1 (and even becomes negative) while the integrand is smooth and mostly always positive) (%i2) plot2d([x*exp((x1)^2), integrate(t*exp((t1)^2), t, inf,x)], [x,5,5])$ (2nd part on screenshot) However, if I simply add and subtract unity, the actual result is recovered: (%i3) integrate( (t1)*exp((t1)^2), t, inf, x)+ integrate( exp((t1)^2), t, inf, x); (%o3) (sqrt(%pi)*erf(x1))/2%e^(x^2+2*x1)/2+sqrt(%pi)/2 (3rd part on screenshot) Looks like additional signum or abs takes place somewhere? Sorry for not giving solution, I am not that skillful in lisp... Thanks for good system, anyway! :)  Maxima version: 5.24.0 Maxima build date: 10:35 8/1/2011 Host type: x86_64pclinuxgnu Lisp implementation type: SBCL Lisp implementation version: 1.0.19gentoo  PS. Just for information, if this can help, run_testsuite() gives 1 error. May be just my installation is wrong? Running tests in rtest16: ********************** Problem 386 *************** Input: closeto(zeta(%i + 3)  (1.10721440843141  .1482908671781754 %i), 1.e15) Result: 3.3157171357748244e9 This differed from the expected result: true   You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3387042&group_id=4933 