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}
(65) 
_{Apr}
(63) 
_{May}
(95) 
_{Jun}
(56) 
_{Jul}
(48) 
_{Aug}
(55) 
_{Sep}
(100) 
_{Oct}
(57) 
_{Nov}
(33) 
_{Dec}
(46) 
2016 
_{Jan}
(76) 
_{Feb}
(53) 
_{Mar}
(88) 
_{Apr}
(79) 
_{May}
(62) 
_{Jun}
(65) 
_{Jul}
(37) 
_{Aug}
(23) 
_{Sep}
(108) 
_{Oct}
(68) 
_{Nov}
(66) 
_{Dec}
(47) 
2017 
_{Jan}
(55) 
_{Feb}
(11) 
_{Mar}
(28) 
_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 





1
(3) 
2
(2) 
3
(6) 
4
(2) 
5

6
(1) 
7
(1) 
8
(2) 
9
(2) 
10

11

12

13

14
(1) 
15
(2) 
16

17

18

19

20
(1) 
21

22

23

24

25

26

27

28

29

30
(3) 
31
(2) 
From: SourceForge.net <noreply@so...>  20111231 20:59:45

Bugs item #3468021, was opened at 20111231 12:59 Message generated for change (Tracker Item Submitted) made by hadronzoo You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3468021&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: Xmaxima or other UI Group: Includes proposed fix Status: Open Resolution: None Priority: 5 Private: No Submitted By: Joshua Griffith (hadronzoo) Assigned to: Nobody/Anonymous (nobody) Summary: imaxima.el uses incorrect tmp directory on OS X Initial Comment: When systemtype is 'darwin, tempdirectory should be "/tmp", not temporaryfiledirectory.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3468021&group_id=4933 
From: SourceForge.net <noreply@so...>  20111231 13:40:23

Bugs item #3467900, was opened at 20111231 05:40 Message generated for change (Tracker Item Submitted) made by stewen You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3467900&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: Roland Stewen (stewen) Assigned to: Nobody/Anonymous (nobody) Summary: fullratsubst  2digits floats behaves unexpected Initial Comment: (%i1) load (lrats)$ (%i2) /* this works fine: */ v1 : [0.5 * t, 1.3 * t]$ (%i3) fullratsubst (v1[1] = x, v1); ... 13 x (%o3) [x, ] 5 (%i4) /* This is unexpected */ v2 : [0.54 * t, 1.31 * t]$ (%i5) fullratsubst (v2[1] = x, v2); ... correct replacements by rat 200 x + 23 t (%o5) [x, ] 100 Can this be a Bug? Maxima version: 5.26.0 Maxima build date: 13:50 12/31/2011 Host type: x86_64unknownlinuxgnu Lisp implementation type: CLISP Lisp implementation version: 2.44.1 (20080223) (built 3487321054) (memory 3534324644) Roland Stewen  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3467900&group_id=4933 
From: SourceForge.net <noreply@so...>  20111230 22:06:01

Bugs item #3467526, was opened at 20111230 10:37 Message generated for change (Comment added) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3467526&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: Ted Woollett (woollett) Assigned to: Nobody/Anonymous (nobody) Summary: expintegral_e: series failed. Initial Comment: float or bfloat of gamma_incomplete fails if second arg is large imag. number.  Maxima 5.25.1 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) fpprintprec:8$ (%i2) block([fpprec:20],bfloat(gamma_incomplete(0,10*%i*log(2)))); (%o2) 1.2240361b1*%i6.9062484b2 (%i3) block([fpprec:20],bfloat(gamma_incomplete(0,100*%i*log(2)))); (%o3) 4.5615238b6*%i5.0797908b5 (%i4) block([fpprec:20],bfloat(gamma_incomplete(0,1000*%i*log(2)))); (%o4) 5.623594b2741.2330921b276*%i (%i5) block([fpprec:20],bfloat(gamma_incomplete(0,10000*%i*log(2)))); expintegral_e: series failed.  an error. To debug this try: debugmode(true); /*I am actually interested in being able to compute the integral */ (%i6) integrate(integrate(cos(x*y*log(2)),y,1,10^4),x,1,10^4); (%o6) (%i*gamma_incomplete(0,100000000*%i*log(2)) 2*%i*gamma_incomplete(0,10000*%i*log(2)) +%i*gamma_incomplete(0,%i*log(2))%i*gamma_incomplete(0,%i*log(2)) +2*%i*gamma_incomplete(0,10000*%i*log(2)) %i*gamma_incomplete(0,100000000*%i*log(2))) /(2*log(2)) /* which involves 10^8*%i*log(2) */ Ted Woollett  >Comment By: Ted Woollett (woollett) Date: 20111230 14:06 Message: Maxima version: 5.25.1 Maxima build date: 10:2 9/6/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 Ted Woollett  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3467526&group_id=4933 
From: SourceForge.net <noreply@so...>  20111230 19:29:40

Bugs item #3454370, was opened at 20111208 04:59 Message generated for change (Comment added) made by dgildea You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3454370&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: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: https://www.google.com/accounts () >Assigned to: Dan Gildea (dgildea) Summary: integrate erf fails Initial Comment: I use maxima 5.25.1 on opensuse 11.3. If h(x,a,b):=erf((x+a)/(b*sqrt(2)))erf((xa)/(b*sqrt(2))); C(a,b):=integrate(h(x,a,b),x,inf,inf); then, C(a,1); I got 0, but with maxima 5.24.0 I got 4a.  >Comment By: Dan Gildea (dgildea) Date: 20111230 11:29 Message: Fixed in defint.lisp  avoid inf  inf => 0 in samesheetsubs  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3454370&group_id=4933 
From: SourceForge.net <noreply@so...>  20111230 18:37:49

Bugs item #3467526, was opened at 20111230 10:37 Message generated for change (Tracker Item Submitted) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3467526&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: Ted Woollett (woollett) Assigned to: Nobody/Anonymous (nobody) Summary: expintegral_e: series failed. Initial Comment: float or bfloat of gamma_incomplete fails if second arg is large imag. number.  Maxima 5.25.1 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) fpprintprec:8$ (%i2) block([fpprec:20],bfloat(gamma_incomplete(0,10*%i*log(2)))); (%o2) 1.2240361b1*%i6.9062484b2 (%i3) block([fpprec:20],bfloat(gamma_incomplete(0,100*%i*log(2)))); (%o3) 4.5615238b6*%i5.0797908b5 (%i4) block([fpprec:20],bfloat(gamma_incomplete(0,1000*%i*log(2)))); (%o4) 5.623594b2741.2330921b276*%i (%i5) block([fpprec:20],bfloat(gamma_incomplete(0,10000*%i*log(2)))); expintegral_e: series failed.  an error. To debug this try: debugmode(true); /*I am actually interested in being able to compute the integral */ (%i6) integrate(integrate(cos(x*y*log(2)),y,1,10^4),x,1,10^4); (%o6) (%i*gamma_incomplete(0,100000000*%i*log(2)) 2*%i*gamma_incomplete(0,10000*%i*log(2)) +%i*gamma_incomplete(0,%i*log(2))%i*gamma_incomplete(0,%i*log(2)) +2*%i*gamma_incomplete(0,10000*%i*log(2)) %i*gamma_incomplete(0,100000000*%i*log(2))) /(2*log(2)) /* which involves 10^8*%i*log(2) */ Ted Woollett  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3467526&group_id=4933 
From: SourceForge.net <noreply@so...>  20111220 14:27:01

Bugs item #3397562, was opened at 20110824 10:31 Message generated for change (Comment added) made by dgildea You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&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: Closed >Resolution: Fixed Priority: 5 Private: No Submitted By: https://www.google.com/accounts () >Assigned to: Dan Gildea (dgildea) Summary: sin(x^2)/(x^2) improper integral incorrect? Initial Comment: (%i1) integrate(sin(x^2)/(x^2),x,1,inf); defint: integral is divergent. This used to just return the integral as a nounform. Numerical approximation says it's about 0.2862504407259549. Thanks!  >Comment By: Dan Gildea (dgildea) Date: 20111220 06:27 Message: Fixed by adding case for complex infinity to routine simplim%gamma_incomplete in gamma.lisp. (%i2) limit( gamma_incomplete(1/2, %i*x^2), x, inf); (%o2) 0 (%i3) integrate(sin(x^2)/(x^2),x,1,inf); (%o3) %i*gamma_incomplete(1/2,%i)/2^(5/2) gamma_incomplete(1/2,%i)/2^(5/2)%i*gamma_incomplete(1/2,%i)/2^(5/2) gamma_incomplete(1/2,%i)/2^(5/2) (%i4) float(%); (%o4) .1767766952966368*(.2733129188747918*%i.5348723621187728)*%i +.1767766952966368*(.2733129188747918*%i.5348723621187728)*%i .1767766952966368*(.2733129188747918*%i.5348723621187728) .1767766952966368*(.2733129188747918*%i.5348723621187728) (%i5) rectform(%); (%o5) .2857366463228523  Comment By: Dan Gildea (dgildea) Date: 20111204 11:59 Message: Seems to be caused by the following problem with limit: (%i2) limit( gamma_incomplete(1/2, %i*x^2), x, inf); (%o2) gamma_incomplete(1/2,infinity)  Comment By: https://www.google.com/accounts () Date: 20111202 07:56 Message: According to Laurent Fousse, it was caused by commit 59775311e53ef8a8fb5a3ad067a6c1cc153075d2 (see http://trac.sagemath.org/sage_trac/ticket/11737)  Comment By: https://www.google.com/accounts () Date: 20111202 07:45 Message: I was using 5.25.0. And (%i3) quad_qagi(sin(x^2)/(x^2),x,1,inf); ***MESSAGE FROM ROUTINE DQAGI IN LIBRARY SLATEC. ***INFORMATIVE MESSAGE, PROG CONTINUES, TRACEBACK REQUESTED * ABNORMAL RETURN * ERROR NUMBER = 1 * ***END OF MESSAGE (%o3) [.2852680980196833, .007186678748743569, 5985, 1] And Wolfram Alpha gives, for input integrate(sin(x^2)/x^2,x,1,inf) sqrt(pi/2) (12 C(sqrt(2/pi)))+sin(1)~~0.285737 where C is a Fresnel integral.  Comment By: Raymond Toy (rtoy) Date: 20110824 10:59 Message: What version are you using? The current version returns %i*gamma_incomplete(1,2*%i)/2+%i*gamma_incomplete(1,2*%i)/2+1/2 This is approximately .6734567682657728. quad_qagi also says the integral is approximately 0.67335. Marking this as pending/worksforme  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&group_id=4933 
From: SourceForge.net <noreply@so...>  20111215 18:42:48

Bugs item #3460131, was opened at 20111215 10:42 Message generated for change (Tracker Item Submitted) made by lvwarren You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3460131&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: Xmaxima or other UI Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Van (lvwarren) Assigned to: Nobody/Anonymous (nobody) Summary: AVG detects Trojan SHeur4.HKI in xmaxima 5.24.0 Initial Comment: Was doing some system maintenance on my Windows PC this morning and AVG raised an alert: AVG detects Trojan SHeur4.HKI in xmaxima 5.24.0. The file was in C:\Program Files\Maxima5.24.0\bin\xmaxima.exe. There is no entry yet for this trojan in the AVG database. Since xmaxima is used for advanced physics I am wondering if this is a centrifuge attack.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3460131&group_id=4933 
From: SourceForge.net <noreply@so...>  20111215 08:17:40

Bugs item #3459556, was opened at 20111214 09:11 Message generated for change (Comment added) made by riotorto You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3459556&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: Edd Barrett (vext0101) Assigned to: Nobody/Anonymous (nobody) Summary: draw3d of z=0: there is no surface in a draw3d range Initial Comment: Consider this function: convex_env(u_l, u_h, v_l, v_h, z_l, z_h) := wxdraw3d( xrange = [u_l, u_h], yrange = [v_l, v_h], zrange = [z_l, z_h], implicit( z=((v_l*v)  (u_l*u) + (u_l*v_l)), u, u_l, u_h, v, v_l, v_h, z, z_l, z_h ) )$ Most applications of this function yield a surface, however, sometimes the result is a single point and maxima bails out. Eg: convex_env(0, 2, 0, 2, 0, 2); gives: draw3d (implicit): no surface within these ranges. Infact it is the same as plotting: implicit( z=0, u, u_l, u_h, v, v_l, v_h, z, z_l, z_h ) Is z=0 not a surface within z=[0, 2]? I get the same when drawing z=0.1, but not for z=0.2. Is this correct ?  >Comment By: Mario Rodriguez Riotorto (riotorto) Date: 20111215 00:17 Message: Hello, I'll take a look at this. In the meantime, you can consider if the following workaround, or something similar, fits your needs: convex_env(ul, uh, vl, vh, zl, zh) := wxdraw3d( xrange = [ul, uh], yrange = [vl, vh], zrange = [zl, zh], implicit( z=((vl*v)  (ul*u) + (ul*vl)), u, ul0.1, uh+0.1, v ,vl0.1, vh+0.1, z, zl0.1, zh+0.1 ) )$  Mario  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3459556&group_id=4933 
From: SourceForge.net <noreply@so...>  20111214 17:11:45

Bugs item #3459556, was opened at 20111214 09:11 Message generated for change (Tracker Item Submitted) made by vext0101 You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3459556&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: Edd Barrett (vext0101) Assigned to: Nobody/Anonymous (nobody) Summary: draw3d of z=0: there is no surface in a draw3d range Initial Comment: Consider this function: convex_env(u_l, u_h, v_l, v_h, z_l, z_h) := wxdraw3d( xrange = [u_l, u_h], yrange = [v_l, v_h], zrange = [z_l, z_h], implicit( z=((v_l*v)  (u_l*u) + (u_l*v_l)), u, u_l, u_h, v, v_l, v_h, z, z_l, z_h ) )$ Most applications of this function yield a surface, however, sometimes the result is a single point and maxima bails out. Eg: convex_env(0, 2, 0, 2, 0, 2); gives: draw3d (implicit): no surface within these ranges. Infact it is the same as plotting: implicit( z=0, u, u_l, u_h, v, v_l, v_h, z, z_l, z_h ) Is z=0 not a surface within z=[0, 2]? I get the same when drawing z=0.1, but not for z=0.2. Is this correct ?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3459556&group_id=4933 
From: SourceForge.net <noreply@so...>  20111209 11:14:56

Bugs item #3455326, was opened at 20111209 03:14 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3455326&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: 4 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: nfloat Initial Comment: Wrong: (should be 0.5) (%i7) nfloat((sqrt(1+x)1)/x,[x=1/10^16],15); (%o7) 0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3455326&group_id=4933 
From: SourceForge.net <noreply@so...>  20111209 05:39:39

Bugs item #3443554, was opened at 20111127 15:31 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3443554&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: Ripley (jensenripley) Assigned to: Nobody/Anonymous (nobody) Summary: Tío different results for an integral Initial Comment: Hello, By integrating sqrt(4*x^2y^2) into the triangle {0<x<1; 0<y<x} in two diffferent order: integrate( integrate( sqrt(4*x^2y^2), x , y, 1), y , 0,1) Integrate( integrate( sqrt(4*x^2y^2), y, 0, x), x, 0,1) WxMaxime give two different results. Where is the mistake?  >Comment By: Raymond Toy (rtoy) Date: 20111208 21:39 Message: You should always check your work. Commercial systems have bugs too. The definite integration routines are in src/defint.lisp. References are included therein along with some description of the individual algorithms. If you're going to debug the code, you'll need to know lisp. If you can even find out that one of the definite or indefinite integrals is wrong would be a great help. Thanks for your help  Comment By: Ripley (jensenripley) Date: 20111208 07:14 Message: I couldn't realize that is simply a mistake... So, should I assume Maxima integrator subroutines aren't trustworthy? I like Maxima but it is not admissible for me to work with a not trusty "scientific" software... Could someone give a hand to have any reference to how Maxima manage integrals internally? I am going to do my best to find the bug...  Comment By: Raymond Toy (rtoy) Date: 20111130 11:05 Message: I have looked very briefly at this. The intermediate indefinite integrals look right, and maxima says the derivative simplifies to the integrand, so I don't really know what's wrong. I think the thing to do is to examine whether the indefinite integrals are actually correct. There could be a problem with choosing the wrong branch since one of the integrals involves inverse trig functions.  Comment By: jcvtrip (jcvtrip) Date: 20111130 09:27 Message: Please, some help on this question? Any comment would be helpful...  Comment By: jcvtrip (jcvtrip) Date: 20111128 10:01 Message: Hello, Here you have a transcript... Maxima version: 5.25.1 Maxima build date: 9:14 9/5/2011 Host type: i686appledarwin10.8.0 Lisp implementation type: SBCL Lisp implementation version: 1.0.47 expres:sqrt(4*x^2y^2); (%o5) sqrt(4*x^2y^2) integrate(integrate(expres, x,y,1),y,0,1); "Is "y1" positive, negative, or zero?"n; "Is "y" positive or negative?"p; (%o6) (3*log(4*sqrt(3)+8)+4*%pi1)/18 ev(%o6,numer); (%o12) 1.093118188249627 integrate(integrate(expres, y, 0,x),x,0,1); "Is "x" positive, negative, or zero?"p; (%o10) (2*%pi+3^(3/2))/18 ev(%o10,numer); (%o13) .6377409849936787  Comment By: Ripley (jensenripley) Date: 20111128 08:28 Message: Thank you rswarbrick but you mistaked the input line ( % i3)... You wrote Integrate(integrate(expr, x,0, x), x,0,1) instead of Integrate(integrate(expr, y, 0, x) ,x ,0,1)  Comment By: Rupert Swarbrick (rswarbrick) Date: 20111128 01:15 Message: Drat. I didn't realise that sourceforge would eat the text formatting. If it isn't clear, the first time I get an answer with logs and sqrts. The second time, I get a more complicated integrate noun form.  Comment By: Rupert Swarbrick (rswarbrick) Date: 20111128 01:13 Message: With a recent version of the sources, I get the following: Maxima 5.25post http://maxima.sourceforge.net using Lisp SBCL 1.0.53.0.debian 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) expr: sqrt(4*x^2y^2); 2 2 (%o1) sqrt(4 x  y ) (%i2) integrate(integrate(expr, x, y, 1), y, 0, 1); Is y  1 positive, negative, or zero? n; Is y positive or negative? pos; 3 log(4 sqrt(3) + 8) + 4 %pi  1 (%o2)  18 (%i3) integrate(integrate(expr, x, 0, x), x, 0, 1); Is x positive, negative, or zero? pos; Is y positive or negative? pos; rat: replaced 1.e8 by 1/100000000 = 1.e8 <snip lots more rat warnings: eugh.> Is y positive or negative? pos; 1 / 2 2 2 2 [ y log(4 sqrt(4 x  y ) + 8 x) y log(4 %i y) (%o3) I (  +  ] 4 4 / 0 2 2 x sqrt(4 x  y ) + ) dx 2 This, while unhelpful isn't wrong. Can you post a transcript of what you actually do and get (use maxima from the command line, maybe). Definitely include the version of Maxima that you're using. Rupert  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3443554&group_id=4933 
From: SourceForge.net <noreply@so...>  20111208 15:14:24

Bugs item #3443554, was opened at 20111127 15:31 Message generated for change (Comment added) made by jensenripley You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3443554&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: Ripley (jensenripley) Assigned to: Nobody/Anonymous (nobody) Summary: Tío different results for an integral Initial Comment: Hello, By integrating sqrt(4*x^2y^2) into the triangle {0<x<1; 0<y<x} in two diffferent order: integrate( integrate( sqrt(4*x^2y^2), x , y, 1), y , 0,1) Integrate( integrate( sqrt(4*x^2y^2), y, 0, x), x, 0,1) WxMaxime give two different results. Where is the mistake?  >Comment By: Ripley (jensenripley) Date: 20111208 07:14 Message: I couldn't realize that is simply a mistake... So, should I assume Maxima integrator subroutines aren't trustworthy? I like Maxima but it is not admissible for me to work with a not trusty "scientific" software... Could someone give a hand to have any reference to how Maxima manage integrals internally? I am going to do my best to find the bug...  Comment By: Raymond Toy (rtoy) Date: 20111130 11:05 Message: I have looked very briefly at this. The intermediate indefinite integrals look right, and maxima says the derivative simplifies to the integrand, so I don't really know what's wrong. I think the thing to do is to examine whether the indefinite integrals are actually correct. There could be a problem with choosing the wrong branch since one of the integrals involves inverse trig functions.  Comment By: jcvtrip (jcvtrip) Date: 20111130 09:27 Message: Please, some help on this question? Any comment would be helpful...  Comment By: jcvtrip (jcvtrip) Date: 20111128 10:01 Message: Hello, Here you have a transcript... Maxima version: 5.25.1 Maxima build date: 9:14 9/5/2011 Host type: i686appledarwin10.8.0 Lisp implementation type: SBCL Lisp implementation version: 1.0.47 expres:sqrt(4*x^2y^2); (%o5) sqrt(4*x^2y^2) integrate(integrate(expres, x,y,1),y,0,1); "Is "y1" positive, negative, or zero?"n; "Is "y" positive or negative?"p; (%o6) (3*log(4*sqrt(3)+8)+4*%pi1)/18 ev(%o6,numer); (%o12) 1.093118188249627 integrate(integrate(expres, y, 0,x),x,0,1); "Is "x" positive, negative, or zero?"p; (%o10) (2*%pi+3^(3/2))/18 ev(%o10,numer); (%o13) .6377409849936787  Comment By: Ripley (jensenripley) Date: 20111128 08:28 Message: Thank you rswarbrick but you mistaked the input line ( % i3)... You wrote Integrate(integrate(expr, x,0, x), x,0,1) instead of Integrate(integrate(expr, y, 0, x) ,x ,0,1)  Comment By: Rupert Swarbrick (rswarbrick) Date: 20111128 01:15 Message: Drat. I didn't realise that sourceforge would eat the text formatting. If it isn't clear, the first time I get an answer with logs and sqrts. The second time, I get a more complicated integrate noun form.  Comment By: Rupert Swarbrick (rswarbrick) Date: 20111128 01:13 Message: With a recent version of the sources, I get the following: Maxima 5.25post http://maxima.sourceforge.net using Lisp SBCL 1.0.53.0.debian 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) expr: sqrt(4*x^2y^2); 2 2 (%o1) sqrt(4 x  y ) (%i2) integrate(integrate(expr, x, y, 1), y, 0, 1); Is y  1 positive, negative, or zero? n; Is y positive or negative? pos; 3 log(4 sqrt(3) + 8) + 4 %pi  1 (%o2)  18 (%i3) integrate(integrate(expr, x, 0, x), x, 0, 1); Is x positive, negative, or zero? pos; Is y positive or negative? pos; rat: replaced 1.e8 by 1/100000000 = 1.e8 <snip lots more rat warnings: eugh.> Is y positive or negative? pos; 1 / 2 2 2 2 [ y log(4 sqrt(4 x  y ) + 8 x) y log(4 %i y) (%o3) I (  +  ] 4 4 / 0 2 2 x sqrt(4 x  y ) + ) dx 2 This, while unhelpful isn't wrong. Can you post a transcript of what you actually do and get (use maxima from the command line, maybe). Definitely include the version of Maxima that you're using. Rupert  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3443554&group_id=4933 
From: SourceForge.net <noreply@so...>  20111208 12:59:33

Bugs item #3454370, was opened at 20111208 04:59 Message generated for change (Tracker Item Submitted) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3454370&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 erf fails Initial Comment: I use maxima 5.25.1 on opensuse 11.3. If h(x,a,b):=erf((x+a)/(b*sqrt(2)))erf((xa)/(b*sqrt(2))); C(a,b):=integrate(h(x,a,b),x,inf,inf); then, C(a,1); I got 0, but with maxima 5.24.0 I got 4a.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3454370&group_id=4933 
From: SourceForge.net <noreply@so...>  20111207 13:16:54

Bugs item #3453553, was opened at 20111207 05:16 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3453553&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: 2 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: subst allows nested lists Initial Comment: The user documentation for subst doesn't mention that the first argument of subst can be a nested list; for example (%i25) subst([[a=2]],a+b); (%o25) b+2 For lists L1 and L2, subst(L1, subst(L2,e)) = subst(append(L1,L2), e) isn't an identity; example: (%i31) subst([[b=1]], subst([[a=1]], a+b)); (%o31) 2 (%i32) subst(append([[b=1]],[[a=1]]),a+b); subst: improper argument: [[b=1],[a=1]] (%i33) subst([b=1], subst([a=1], a+b)); (%o33) 2 (%i34) subst(append([b=1],[a=1]),a+b); (%o34) 2  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3453553&group_id=4933 
From: SourceForge.net <noreply@so...>  20111206 14:26:19

Bugs item #3452434, was opened at 20111206 06:26 Message generated for change (Tracker Item Submitted) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3452434&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: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: Contour integral apparently depends on center of circle Initial Comment: See the following sagesupport thread: http://groups.google.com/group/sagesupport/browse_thread/thread/3aee0a551bfac1d9/a0143d963486e38f Apparently this integral evaluation depends on what circle containing the (only) pole one uses. I haven't checked the math by hand personally, just reporting. Maxima 5.25.1 http://maxima.sourceforge.net using Lisp SBCL 1.0.24 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) f(z):=(z%i)*(z1)^2/(z(1/2%i/3)); 2 (z  %i) (z  1) (%o1) f(z) :=   1 %i z  (  ) 2 3 (%i2) display2d:false; (%o2) false (%i3) integrate(f(1/2%i/3+exp(%i*t))*%i*exp(%i*t),t,0,2*%pi); (%o3) %i*((4*%i13)/3((724*%i57)*%pi+144*%i468)/108) (%i4) expand(ratsimp(%i*((4*%i13)/3((724*%i57)*%pi+144*%i468)/108))); (%o4) 19*%i*%pi/36+181*%pi/27 (%i5) integrate(f(2*exp(%i*t))*2*%i*exp(%i*t),t,0,2*%pi); (%o5) 0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3452434&group_id=4933 
From: SourceForge.net <noreply@so...>  20111204 19:59:50

Bugs item #3397562, was opened at 20110824 10:31 Message generated for change (Comment added) made by dgildea You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&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: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: sin(x^2)/(x^2) improper integral incorrect? Initial Comment: (%i1) integrate(sin(x^2)/(x^2),x,1,inf); defint: integral is divergent. This used to just return the integral as a nounform. Numerical approximation says it's about 0.2862504407259549. Thanks!  >Comment By: Dan Gildea (dgildea) Date: 20111204 11:59 Message: Seems to be caused by the following problem with limit: (%i2) limit( gamma_incomplete(1/2, %i*x^2), x, inf); (%o2) gamma_incomplete(1/2,infinity)  Comment By: https://www.google.com/accounts () Date: 20111202 07:56 Message: According to Laurent Fousse, it was caused by commit 59775311e53ef8a8fb5a3ad067a6c1cc153075d2 (see http://trac.sagemath.org/sage_trac/ticket/11737)  Comment By: https://www.google.com/accounts () Date: 20111202 07:45 Message: I was using 5.25.0. And (%i3) quad_qagi(sin(x^2)/(x^2),x,1,inf); ***MESSAGE FROM ROUTINE DQAGI IN LIBRARY SLATEC. ***INFORMATIVE MESSAGE, PROG CONTINUES, TRACEBACK REQUESTED * ABNORMAL RETURN * ERROR NUMBER = 1 * ***END OF MESSAGE (%o3) [.2852680980196833, .007186678748743569, 5985, 1] And Wolfram Alpha gives, for input integrate(sin(x^2)/x^2,x,1,inf) sqrt(pi/2) (12 C(sqrt(2/pi)))+sin(1)~~0.285737 where C is a Fresnel integral.  Comment By: Raymond Toy (rtoy) Date: 20110824 10:59 Message: What version are you using? The current version returns %i*gamma_incomplete(1,2*%i)/2+%i*gamma_incomplete(1,2*%i)/2+1/2 This is approximately .6734567682657728. quad_qagi also says the integral is approximately 0.67335. Marking this as pending/worksforme  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&group_id=4933 
From: SourceForge.net <noreply@so...>  20111204 12:03:18

Bugs item #3450204, was opened at 20111204 04:03 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3450204&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: 3 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: subst allows more than 3 arguments Initial Comment: Generally, Maxima gives an error when the number of arguments is wrong. But not subst: (%i1) subst(a,b,x+b,p,q,w); (%o1) x+a  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3450204&group_id=4933 
From: SourceForge.net <noreply@so...>  20111203 18:54:48

Bugs item #3442168, was opened at 20111125 15:37 Message generated for change (Comment added) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3442168&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: Ted Woollett (woollett) Assigned to: Nobody/Anonymous (nobody) Summary: real and imaginary parts asech(x) Initial Comment: Despite asech(x) having both real and imaginary parts, (%i1) realpart(asech(x)); (%o1) asech(x) (%i2) imagpart(asech(x)); (%o2) 0  >Comment By: Ted Woollett (woollett) Date: 20111203 10:54 Message: system info: Maxima version: 5.25.1 Maxima build date: 10:2 9/6/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3442168&group_id=4933 
From: SourceForge.net <noreply@so...>  20111203 18:53:54

Bugs item #3442167, was opened at 20111125 15:35 Message generated for change (Comment added) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3442167&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: Ted Woollett (woollett) Assigned to: Nobody/Anonymous (nobody) Summary: real and imaginary parts acosh Initial Comment: despite acosh(x) for real x having both real and imaginary parts, (%i1) realpart(acosh(x)); (%o1) acosh(x) (%i2) imagpart(acosh(x)); (%o2) 0  >Comment By: Ted Woollett (woollett) Date: 20111203 10:53 Message: system info: Maxima version: 5.25.1 Maxima build date: 10:2 9/6/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3442167&group_id=4933 
From: SourceForge.net <noreply@so...>  20111203 18:52:22

Bugs item #3442165, was opened at 20111125 15:32 Message generated for change (Comment added) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3442165&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: Ted Woollett (woollett) Assigned to: Nobody/Anonymous (nobody) Summary: inappropriate messages from integrate Initial Comment: (%i1) integrate(acos(%i*x),x,2,3); Is W35620 positive or negative? pos; Is W37394 positive or negative? pos; (%o1) %i*(3*log(sqrt(10)+3)2*log(sqrt(5)2)+sqrt(10)sqrt(5))+5*%pi/2  >Comment By: Ted Woollett (woollett) Date: 20111203 10:52 Message: system info: Maxima version: 5.25.1 Maxima build date: 10:2 9/6/2011 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3442165&group_id=4933 
From: SourceForge.net <noreply@so...>  20111203 16:21:39

Bugs item #3449468, was opened at 20111203 08:20 Message generated for change (Settings changed) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3449468&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: Duplicate Priority: 5 Private: No Submitted By: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: nfloat Initial Comment: Wrong: (%i32) nfloat((sqrt(1+x)1)/x,[x=1.0e16],15); (%o32) 0 OK: (%i33) nfloat((sqrt(1+x)1)/x,[x=1.0e16],16); (%o33) 5.0b1  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3449468&group_id=4933 
From: SourceForge.net <noreply@so...>  20111203 16:20:11

Bugs item #3449468, was opened at 20111203 08:20 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3449468&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: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: nfloat Initial Comment: Wrong: (%i32) nfloat((sqrt(1+x)1)/x,[x=1.0e16],15); (%o32) 0 OK: (%i33) nfloat((sqrt(1+x)1)/x,[x=1.0e16],16); (%o33) 5.0b1  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3449468&group_id=4933 
From: SourceForge.net <noreply@so...>  20111203 13:44:41

Bugs item #3449369, was opened at 20111203 05:44 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3449369&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: Barton Willis (willisbl) Assigned to: Barton Willis (willisbl) Summary: nfloat Initial Comment: Wrong: (%i32) nfloat((sqrt(1+x)1)/x,[x=1.0e16],15); (%o32) 0 OK: (%i33) nfloat((sqrt(1+x)1)/x,[x=1.0e16],16); (%o33) 5.0b1  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3449369&group_id=4933 
From: SourceForge.net <noreply@so...>  20111202 15:56:52

Bugs item #3397562, was opened at 20110824 10:31 Message generated for change (Comment added) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&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: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: sin(x^2)/(x^2) improper integral incorrect? Initial Comment: (%i1) integrate(sin(x^2)/(x^2),x,1,inf); defint: integral is divergent. This used to just return the integral as a nounform. Numerical approximation says it's about 0.2862504407259549. Thanks!  >Comment By: https://www.google.com/accounts () Date: 20111202 07:56 Message: According to Laurent Fousse, it was caused by commit 59775311e53ef8a8fb5a3ad067a6c1cc153075d2 (see http://trac.sagemath.org/sage_trac/ticket/11737)  Comment By: https://www.google.com/accounts () Date: 20111202 07:45 Message: I was using 5.25.0. And (%i3) quad_qagi(sin(x^2)/(x^2),x,1,inf); ***MESSAGE FROM ROUTINE DQAGI IN LIBRARY SLATEC. ***INFORMATIVE MESSAGE, PROG CONTINUES, TRACEBACK REQUESTED * ABNORMAL RETURN * ERROR NUMBER = 1 * ***END OF MESSAGE (%o3) [.2852680980196833, .007186678748743569, 5985, 1] And Wolfram Alpha gives, for input integrate(sin(x^2)/x^2,x,1,inf) sqrt(pi/2) (12 C(sqrt(2/pi)))+sin(1)~~0.285737 where C is a Fresnel integral.  Comment By: Raymond Toy (rtoy) Date: 20110824 10:59 Message: What version are you using? The current version returns %i*gamma_incomplete(1,2*%i)/2+%i*gamma_incomplete(1,2*%i)/2+1/2 This is approximately .6734567682657728. quad_qagi also says the integral is approximately 0.67335. Marking this as pending/worksforme  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&group_id=4933 
From: SourceForge.net <noreply@so...>  20111202 15:45:53

Bugs item #3397562, was opened at 20110824 10:31 Message generated for change (Comment added) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&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: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: sin(x^2)/(x^2) improper integral incorrect? Initial Comment: (%i1) integrate(sin(x^2)/(x^2),x,1,inf); defint: integral is divergent. This used to just return the integral as a nounform. Numerical approximation says it's about 0.2862504407259549. Thanks!  >Comment By: https://www.google.com/accounts () Date: 20111202 07:45 Message: I was using 5.25.0. And (%i3) quad_qagi(sin(x^2)/(x^2),x,1,inf); ***MESSAGE FROM ROUTINE DQAGI IN LIBRARY SLATEC. ***INFORMATIVE MESSAGE, PROG CONTINUES, TRACEBACK REQUESTED * ABNORMAL RETURN * ERROR NUMBER = 1 * ***END OF MESSAGE (%o3) [.2852680980196833, .007186678748743569, 5985, 1] And Wolfram Alpha gives, for input integrate(sin(x^2)/x^2,x,1,inf) sqrt(pi/2) (12 C(sqrt(2/pi)))+sin(1)~~0.285737 where C is a Fresnel integral.  Comment By: Raymond Toy (rtoy) Date: 20110824 10:59 Message: What version are you using? The current version returns %i*gamma_incomplete(1,2*%i)/2+%i*gamma_incomplete(1,2*%i)/2+1/2 This is approximately .6734567682657728. quad_qagi also says the integral is approximately 0.67335. Marking this as pending/worksforme  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3397562&group_id=4933 