From: SourceForge.net <no...@so...> - 2012-01-14 20:09:34
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Tracker Item Submitted) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-01-14 20:22:08
|
Bugs item #3473860, was opened at 2012-01-14 12:22 Message generated for change (Tracker Item Submitted) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473860&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473860&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-01-14 21:04:39
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Comment added) made by dgildea You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- >Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-01-15 17:59:07
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Comment added) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 09:59 Message: With the latest version of the Maxima again wrong: (%i1) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 (%i2) domain:complex; (%o2) complex (%i3) integrate(x^(1/3)*sin(x),x,0,1); (%o3) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i4) float(%),expand; (%o4) -0.40863415733948 (%i5) quad_qags(x^(1/3)*sin(x),x,0,1); (%o5) [0.39140213472074,2.2152992977098117*10^-9,231,0] What I am doing wrong? Aleksas D ---------------------------------------------------------------------- Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-01-15 19:47:24
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Comment added) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 11:47 Message: Using 5.26.0 (see mailing list) and just adding the latest versions of defint.lisp and limit.lisp (no rebuild) , together with domain:complex, results in the correct answer here: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o1) 1.091879396853239 (%i2) integrate(x^(1/3)*sin(x),x,1,2); (%o2) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -%i*gamma_incomplete(1/3,%i)/12 +%i*gamma_incomplete(1/3,-%i)/12 -%i*gamma_incomplete(1/3,-2*%i)/12 (%i3) expand(float(rectform(%))); (%o3) 0.96415974107533 (%i4) load("limit-new.lisp"); (%o4) "c:/work2/limit-new.lisp" (%i5) load("defint-new.lisp"); (%o5) "c:/work2/defint-new.lisp" (%i6) domain:false$ (%i7) integrate(x^(1/3)*sin(x),x,1,2); (%o7) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -gamma_incomplete(1/3,2*%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,%i)/12 +gamma_incomplete(1/3,%i)/(4*sqrt(3)) +%i*gamma_incomplete(1/3,-%i)/12 +gamma_incomplete(1/3,-%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-2*%i)/12 -gamma_incomplete(1/3,-2*%i)/(4*sqrt(3)) (%i8) expand(float(rectform(%))); (%o8) 1.09187939685324 (%i9) build_info(); Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 09:59 Message: With the latest version of the Maxima again wrong: (%i1) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 (%i2) domain:complex; (%o2) complex (%i3) integrate(x^(1/3)*sin(x),x,0,1); (%o3) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i4) float(%),expand; (%o4) -0.40863415733948 (%i5) quad_qags(x^(1/3)*sin(x),x,0,1); (%o5) [0.39140213472074,2.2152992977098117*10^-9,231,0] What I am doing wrong? Aleksas D ---------------------------------------------------------------------- Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-01-15 20:51:43
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Comment added) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 12:51 Message: However, some other integrals apparently need current git, rather than just 5.26.0 plus updated limit.lisp and defint.lisp. Two examples of history stack overflow errors are shown below: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) load("limit-new.lisp"); (%o1) "c:/work2/limit-new.lisp" (%i2) load("defint-new.lisp"); (%o2) "c:/work2/defint-new.lisp" (%i3) cfloat(z):= expand(float(rectform(z)))$ (%i4) domain:complex$ (%i5) first(quad_qags(x^3,x,0,1)); (%o5) 0.25 (%i6) integrate(x^3,x,0,1); (%o6) 1/4 (%i7) first(quad_qags(x^5,x,0,1)); (%o7) 0.16666666666667 (%i8) integrate(x^5,x,0,1); (%o8) 1/6 (%i9) float(%); (%o9) 0.16666666666667 (%i10) first(quad_qags(exp(x^3),x,0,1)); (%o10) 1.34190441797742 (%i11) integrate(exp(x^3),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i12) integrate(exp(sqrt(x^3)),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i13) first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o13) 0.39140213472074 (%i14) integrate(x^(1/3)*sin(x),x,0,1); (%o14) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i15) cfloat(%); (%o15) 0.3914021347234 Ted Woollett ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 11:47 Message: Using 5.26.0 (see mailing list) and just adding the latest versions of defint.lisp and limit.lisp (no rebuild) , together with domain:complex, results in the correct answer here: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o1) 1.091879396853239 (%i2) integrate(x^(1/3)*sin(x),x,1,2); (%o2) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -%i*gamma_incomplete(1/3,%i)/12 +%i*gamma_incomplete(1/3,-%i)/12 -%i*gamma_incomplete(1/3,-2*%i)/12 (%i3) expand(float(rectform(%))); (%o3) 0.96415974107533 (%i4) load("limit-new.lisp"); (%o4) "c:/work2/limit-new.lisp" (%i5) load("defint-new.lisp"); (%o5) "c:/work2/defint-new.lisp" (%i6) domain:false$ (%i7) integrate(x^(1/3)*sin(x),x,1,2); (%o7) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -gamma_incomplete(1/3,2*%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,%i)/12 +gamma_incomplete(1/3,%i)/(4*sqrt(3)) +%i*gamma_incomplete(1/3,-%i)/12 +gamma_incomplete(1/3,-%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-2*%i)/12 -gamma_incomplete(1/3,-2*%i)/(4*sqrt(3)) (%i8) expand(float(rectform(%))); (%o8) 1.09187939685324 (%i9) build_info(); Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 09:59 Message: With the latest version of the Maxima again wrong: (%i1) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 (%i2) domain:complex; (%o2) complex (%i3) integrate(x^(1/3)*sin(x),x,0,1); (%o3) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i4) float(%),expand; (%o4) -0.40863415733948 (%i5) quad_qags(x^(1/3)*sin(x),x,0,1); (%o5) [0.39140213472074,2.2152992977098117*10^-9,231,0] What I am doing wrong? Aleksas D ---------------------------------------------------------------------- Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-01-15 21:03:05
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Comment added) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- >Comment By: Aleksas (aleksasd) Date: 2012-01-15 13:03 Message: For integrate(x^(1/3)*sin(x),x,1,2) I got a moore nice solution (compare with Mathematica or Maple solutions): (%i1) %emode:false$ (%i2) (%e^(-(%i*%pi)/6)*gamma_incomplete(4/3, %i ))/2 +(%e^( (%i*%pi)/6)*gamma_incomplete(4/3, -%i ))/2 -(%e^( (%i*%pi)/6)*gamma_incomplete(4/3,-2*%i))/2 -(%e^(-(%i*%pi)/6)*gamma_incomplete(4/3, 2*%i))/2$ (%i3) rectform(%)$ (%i4) float(%), numer; (%o4) 1.091879396853239 (%i5) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o5) 1.091879396853239 Aleksas D ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 12:51 Message: However, some other integrals apparently need current git, rather than just 5.26.0 plus updated limit.lisp and defint.lisp. Two examples of history stack overflow errors are shown below: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) load("limit-new.lisp"); (%o1) "c:/work2/limit-new.lisp" (%i2) load("defint-new.lisp"); (%o2) "c:/work2/defint-new.lisp" (%i3) cfloat(z):= expand(float(rectform(z)))$ (%i4) domain:complex$ (%i5) first(quad_qags(x^3,x,0,1)); (%o5) 0.25 (%i6) integrate(x^3,x,0,1); (%o6) 1/4 (%i7) first(quad_qags(x^5,x,0,1)); (%o7) 0.16666666666667 (%i8) integrate(x^5,x,0,1); (%o8) 1/6 (%i9) float(%); (%o9) 0.16666666666667 (%i10) first(quad_qags(exp(x^3),x,0,1)); (%o10) 1.34190441797742 (%i11) integrate(exp(x^3),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i12) integrate(exp(sqrt(x^3)),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i13) first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o13) 0.39140213472074 (%i14) integrate(x^(1/3)*sin(x),x,0,1); (%o14) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i15) cfloat(%); (%o15) 0.3914021347234 Ted Woollett ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 11:47 Message: Using 5.26.0 (see mailing list) and just adding the latest versions of defint.lisp and limit.lisp (no rebuild) , together with domain:complex, results in the correct answer here: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o1) 1.091879396853239 (%i2) integrate(x^(1/3)*sin(x),x,1,2); (%o2) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -%i*gamma_incomplete(1/3,%i)/12 +%i*gamma_incomplete(1/3,-%i)/12 -%i*gamma_incomplete(1/3,-2*%i)/12 (%i3) expand(float(rectform(%))); (%o3) 0.96415974107533 (%i4) load("limit-new.lisp"); (%o4) "c:/work2/limit-new.lisp" (%i5) load("defint-new.lisp"); (%o5) "c:/work2/defint-new.lisp" (%i6) domain:false$ (%i7) integrate(x^(1/3)*sin(x),x,1,2); (%o7) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -gamma_incomplete(1/3,2*%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,%i)/12 +gamma_incomplete(1/3,%i)/(4*sqrt(3)) +%i*gamma_incomplete(1/3,-%i)/12 +gamma_incomplete(1/3,-%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-2*%i)/12 -gamma_incomplete(1/3,-2*%i)/(4*sqrt(3)) (%i8) expand(float(rectform(%))); (%o8) 1.09187939685324 (%i9) build_info(); Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 09:59 Message: With the latest version of the Maxima again wrong: (%i1) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 (%i2) domain:complex; (%o2) complex (%i3) integrate(x^(1/3)*sin(x),x,0,1); (%o3) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i4) float(%),expand; (%o4) -0.40863415733948 (%i5) quad_qags(x^(1/3)*sin(x),x,0,1); (%o5) [0.39140213472074,2.2152992977098117*10^-9,231,0] What I am doing wrong? Aleksas D ---------------------------------------------------------------------- Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-03-16 15:56:53
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Settings changed) made by dgildea You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 13:03 Message: For integrate(x^(1/3)*sin(x),x,1,2) I got a moore nice solution (compare with Mathematica or Maple solutions): (%i1) %emode:false$ (%i2) (%e^(-(%i*%pi)/6)*gamma_incomplete(4/3, %i ))/2 +(%e^( (%i*%pi)/6)*gamma_incomplete(4/3, -%i ))/2 -(%e^( (%i*%pi)/6)*gamma_incomplete(4/3,-2*%i))/2 -(%e^(-(%i*%pi)/6)*gamma_incomplete(4/3, 2*%i))/2$ (%i3) rectform(%)$ (%i4) float(%), numer; (%o4) 1.091879396853239 (%i5) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o5) 1.091879396853239 Aleksas D ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 12:51 Message: However, some other integrals apparently need current git, rather than just 5.26.0 plus updated limit.lisp and defint.lisp. Two examples of history stack overflow errors are shown below: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) load("limit-new.lisp"); (%o1) "c:/work2/limit-new.lisp" (%i2) load("defint-new.lisp"); (%o2) "c:/work2/defint-new.lisp" (%i3) cfloat(z):= expand(float(rectform(z)))$ (%i4) domain:complex$ (%i5) first(quad_qags(x^3,x,0,1)); (%o5) 0.25 (%i6) integrate(x^3,x,0,1); (%o6) 1/4 (%i7) first(quad_qags(x^5,x,0,1)); (%o7) 0.16666666666667 (%i8) integrate(x^5,x,0,1); (%o8) 1/6 (%i9) float(%); (%o9) 0.16666666666667 (%i10) first(quad_qags(exp(x^3),x,0,1)); (%o10) 1.34190441797742 (%i11) integrate(exp(x^3),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i12) integrate(exp(sqrt(x^3)),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i13) first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o13) 0.39140213472074 (%i14) integrate(x^(1/3)*sin(x),x,0,1); (%o14) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i15) cfloat(%); (%o15) 0.3914021347234 Ted Woollett ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 11:47 Message: Using 5.26.0 (see mailing list) and just adding the latest versions of defint.lisp and limit.lisp (no rebuild) , together with domain:complex, results in the correct answer here: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o1) 1.091879396853239 (%i2) integrate(x^(1/3)*sin(x),x,1,2); (%o2) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -%i*gamma_incomplete(1/3,%i)/12 +%i*gamma_incomplete(1/3,-%i)/12 -%i*gamma_incomplete(1/3,-2*%i)/12 (%i3) expand(float(rectform(%))); (%o3) 0.96415974107533 (%i4) load("limit-new.lisp"); (%o4) "c:/work2/limit-new.lisp" (%i5) load("defint-new.lisp"); (%o5) "c:/work2/defint-new.lisp" (%i6) domain:false$ (%i7) integrate(x^(1/3)*sin(x),x,1,2); (%o7) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -gamma_incomplete(1/3,2*%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,%i)/12 +gamma_incomplete(1/3,%i)/(4*sqrt(3)) +%i*gamma_incomplete(1/3,-%i)/12 +gamma_incomplete(1/3,-%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-2*%i)/12 -gamma_incomplete(1/3,-2*%i)/(4*sqrt(3)) (%i8) expand(float(rectform(%))); (%o8) 1.09187939685324 (%i9) build_info(); Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 09:59 Message: With the latest version of the Maxima again wrong: (%i1) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 (%i2) domain:complex; (%o2) complex (%i3) integrate(x^(1/3)*sin(x),x,0,1); (%o3) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i4) float(%),expand; (%o4) -0.40863415733948 (%i5) quad_qags(x^(1/3)*sin(x),x,0,1); (%o5) [0.39140213472074,2.2152992977098117*10^-9,231,0] What I am doing wrong? Aleksas D ---------------------------------------------------------------------- Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-03-16 15:57:26
|
Bugs item #3473860, was opened at 2012-01-14 12:22 Message generated for change (Settings changed) made by dgildea You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473860&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: Duplicate Priority: 5 Private: No Submitted By: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473860&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-08-15 15:46:41
|
Bugs item #3473858, was opened at 2012-01-14 12:09 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&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: Works For Me Priority: 5 Private: No Submitted By: Aleksas (aleksasd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x^(1/3)*sin(x),x,0,1) Initial Comment: Wrong (%i1) integrate(x^(1/3)*sin(x),x,0,1); (%o1) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i2) float(%),expand; (%o2) -0.408634157339482 (%i3) /* test */ first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o3) 0.391402134720738 Correct exact answer(compare with Mathematica or Maple solutions): (%i4) ((%i-sqrt(3))*gamma_incomplete(4/3,%i))/4-((%i+sqrt(3))*gamma_incomplete(4/3,-%i))/4+gamma(1/3)/(2*sqrt(3))$ (%i5) float(%),expand; (%o5) 0.391402134723405 (%i6) build_info()$ Maxima version: 5.25.1 Maxima build date: 9:23 9/6/2011 Host type: i686-pc-mingw32 Lisp implementation type: Clozure Common Lisp Lisp implementation version: Version 1.7-r14925M (WindowsX8632) Aleksas D. ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 13:03 Message: For integrate(x^(1/3)*sin(x),x,1,2) I got a moore nice solution (compare with Mathematica or Maple solutions): (%i1) %emode:false$ (%i2) (%e^(-(%i*%pi)/6)*gamma_incomplete(4/3, %i ))/2 +(%e^( (%i*%pi)/6)*gamma_incomplete(4/3, -%i ))/2 -(%e^( (%i*%pi)/6)*gamma_incomplete(4/3,-2*%i))/2 -(%e^(-(%i*%pi)/6)*gamma_incomplete(4/3, 2*%i))/2$ (%i3) rectform(%)$ (%i4) float(%), numer; (%o4) 1.091879396853239 (%i5) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o5) 1.091879396853239 Aleksas D ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 12:51 Message: However, some other integrals apparently need current git, rather than just 5.26.0 plus updated limit.lisp and defint.lisp. Two examples of history stack overflow errors are shown below: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) load("limit-new.lisp"); (%o1) "c:/work2/limit-new.lisp" (%i2) load("defint-new.lisp"); (%o2) "c:/work2/defint-new.lisp" (%i3) cfloat(z):= expand(float(rectform(z)))$ (%i4) domain:complex$ (%i5) first(quad_qags(x^3,x,0,1)); (%o5) 0.25 (%i6) integrate(x^3,x,0,1); (%o6) 1/4 (%i7) first(quad_qags(x^5,x,0,1)); (%o7) 0.16666666666667 (%i8) integrate(x^5,x,0,1); (%o8) 1/6 (%i9) float(%); (%o9) 0.16666666666667 (%i10) first(quad_qags(exp(x^3),x,0,1)); (%o10) 1.34190441797742 (%i11) integrate(exp(x^3),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i12) integrate(exp(sqrt(x^3)),x,0,1); Maxima encountered a Lisp error: Error in SETQ [or a callee]: Invocation history stack overflow. Automatically continuing. To enable the Lisp debugger set *debugger-hook* to nil. (%i13) first(quad_qags(x^(1/3)*sin(x),x,0,1)); (%o13) 0.39140213472074 (%i14) integrate(x^(1/3)*sin(x),x,0,1); (%o14) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i15) cfloat(%); (%o15) 0.3914021347234 Ted Woollett ---------------------------------------------------------------------- Comment By: Ted Woollett (woollett) Date: 2012-01-15 11:47 Message: Using 5.26.0 (see mailing list) and just adding the latest versions of defint.lisp and limit.lisp (no rebuild) , together with domain:complex, results in the correct answer here: Maxima 5.26.0 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) (%i1) first(quad_qags(x^(1/3)*sin(x),x,1,2)); (%o1) 1.091879396853239 (%i2) integrate(x^(1/3)*sin(x),x,1,2); (%o2) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -%i*gamma_incomplete(1/3,%i)/12 +%i*gamma_incomplete(1/3,-%i)/12 -%i*gamma_incomplete(1/3,-2*%i)/12 (%i3) expand(float(rectform(%))); (%o3) 0.96415974107533 (%i4) load("limit-new.lisp"); (%o4) "c:/work2/limit-new.lisp" (%i5) load("defint-new.lisp"); (%o5) "c:/work2/defint-new.lisp" (%i6) domain:false$ (%i7) integrate(x^(1/3)*sin(x),x,1,2); (%o7) -2^(1/3)*cos(2)+cos(1)+%i*gamma_incomplete(1/3,2*%i)/12 -gamma_incomplete(1/3,2*%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,%i)/12 +gamma_incomplete(1/3,%i)/(4*sqrt(3)) +%i*gamma_incomplete(1/3,-%i)/12 +gamma_incomplete(1/3,-%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-2*%i)/12 -gamma_incomplete(1/3,-2*%i)/(4*sqrt(3)) (%i8) expand(float(rectform(%))); (%o8) 1.09187939685324 (%i9) build_info(); Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-01-15 09:59 Message: With the latest version of the Maxima again wrong: (%i1) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8 (%i2) domain:complex; (%o2) complex (%i3) integrate(x^(1/3)*sin(x),x,0,1); (%o3) -cos(1)+(%i*gamma_incomplete(1/3,%i))/12-(%i*gamma_incomplete(1/3,-%i))/12 (%i4) float(%),expand; (%o4) -0.40863415733948 (%i5) quad_qags(x^(1/3)*sin(x),x,0,1); (%o5) [0.39140213472074,2.2152992977098117*10^-9,231,0] What I am doing wrong? Aleksas D ---------------------------------------------------------------------- Comment By: Dan Gildea (dgildea) Date: 2012-01-14 13:04 Message: Seems OK in current git with domain:complex. (%i4) domain:complex; (%o4) complex (%i5) integrate(x^(1/3)*sin(x),x,0,1); (%o5) -cos(1)+%i*gamma_incomplete(1/3,%i)/12 -gamma_incomplete(1/3,%i)/(4*sqrt(3)) -%i*gamma_incomplete(1/3,-%i)/12 -gamma_incomplete(1/3,-%i)/(4*sqrt(3))+gamma(1/3)/(2*sqrt(3)) (%i6) float(%),expand; (%o6) .3914021347234045 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3473858&group_id=4933 |