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}
(30) 
_{Apr}
(12) 
_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 







1
(1) 
2

3
(1) 
4

5

6

7

8
(1) 
9
(1) 
10
(2) 
11

12

13

14

15

16

17

18

19
(1) 
20
(4) 
21
(6) 
22

23
(1) 
24

25

26

27

28

29
(2) 
30







From: SourceForge.net <noreply@so...>  20120929 21:47:59

Bugs item #3572973, was opened at 20120929 04:06 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3572973&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: Barton Willis (willisbl) Assigned to: Raymond Toy (rtoy) Summary: missing bigfloat::signum for bigfloats Initial Comment: (%i2) y : 5.6b0$ (%i3) :lisp(bigfloat::signum (meval '$y)) Maxima encountered a Lisp error: Error in BIGFLOATIMPL:SIGNUM [or a callee]: No matching method for the genericfunction #<compiledclosure BIGFLOATIMPL:SIGNUM>,  >Comment By: Raymond Toy (rtoy) Date: 20120929 14:47 Message: $y is not a bigfloat object. You need to do something like :lisp (bigfloat:signum (bigfloat:to $y)) > 1.0b0 Did this show up in some other way?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3572973&group_id=4933 
From: SourceForge.net <noreply@so...>  20120929 11:06:36

Bugs item #3572973, was opened at 20120929 04:06 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3572973&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: Barton Willis (willisbl) Assigned to: Raymond Toy (rtoy) Summary: missing bigfloat::signum for bigfloats Initial Comment: (%i2) y : 5.6b0$ (%i3) :lisp(bigfloat::signum (meval '$y)) Maxima encountered a Lisp error: Error in BIGFLOATIMPL:SIGNUM [or a callee]: No matching method for the genericfunction #<compiledclosure BIGFLOATIMPL:SIGNUM>,  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3572973&group_id=4933 
From: SourceForge.net <noreply@so...>  20120923 11:28:48

Bugs item #3570214, was opened at 20120920 12:22 Message generated for change (Comment added) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570214&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  Plotting Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Manuel GonzálezHidalgo (manuelgonha) Assigned to: Nobody/Anonymous (nobody) Summary: polar graphs are wrong Initial Comment:  build_info("5.28.02","20120827 23:16:48","i686pcmingw32","GNU Common Lisp (GCL)","GCL 2.6.8")  If I put in the console wxplot2d(t, [t,0,4*%pi],[plot_format, gnuplot], [gnuplot_preamble, "set polar; set trange [0:4*pi]; set rrange [0:4*pi];"])$ we obtain a incorrect plot of the function as we can see in the attached image. If I put in the gnuplot: set polar; set rrange [0:4*pi]; set trange [0:4*pi]; plot t; I obtain a coorect plot of the function. But in any case, that is for any function in polar coordinates, in Maxima the graph of the function isn't correct. Thank you very much for your attention, and waiting for a solution to this problem, says goodbye with a greeting and best wishes Manuel González Prof. University of the Balearic Islands  Comment By: Aleksas (aleksasd) Date: 20120923 04:28 Message: All correct. Do (%i1) wxplot2d([ph], [ph,0,4*%pi],[y,4*%pi,4*%pi], [gnuplot_preamble, "set polar; set zeroaxis;"], [xlabel, ""],[ ylabel, ""])$ (%t1) << Graphics >> or (%i2) load(draw)$ (%i3) wxdraw2d(user_preamble = "set grid polar", nticks = 200, grid=true, xrange = [12,12], yrange = [12,12], color = blue, line_width = 2, title = "Spiral", polar(theta,theta,0,4*%pi), proportional_axes = xy )$ (%t3) << Graphics >> best Aleksas D  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570214&group_id=4933 
From: SourceForge.net <noreply@so...>  20120921 05:48:00

Bugs item #3337674, was opened at 20110627 08:46 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3337674&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: Wont Fix Priority: 5 Private: No Submitted By: Dženan Zukić (dzenanz) Assigned to: Nobody/Anonymous (nobody) Summary: Symmetric matrix yields complex eigenvalues Initial Comment: When using eigenvectors command in wxMaxima, the following symmetric matrix yields complex eigenvalues: matrix([2621.4397,7823.3599,1111.2726],[7823.3599,23347.842,3316.4543],[1111.2726,3316.4543,471.08722]) All eigenvalues of a symmetric matrix should be real: http://en.wikipedia.org/wiki/Symmetric_matrix 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: Raymond Toy (rtoy) Date: 20120818 16:52 Message: In addition, I think algorithms for symmetric matrices should be used, instead of a general eigen solver. I don't consider this a bug in maxima. Marking as pending/wontfix.  Comment By: Barton Willis (willisbl) Date: 20110628 21:00 Message: For a floating point evaluation of eigenvalues, you should use a purely numeric method, not a symbolic method. One (not the only) option is eigens_by_jacobi (symmetric and either binary64 or bigfloat entries).  Comment By: Dženan Zukić (dzenanz) Date: 20110628 05:37 Message: Thanks for suggestions, but I was using Maxima trying to verify some results obtained using numeric library. However after getting this nonsensical result from Maxima I used another numeric library and obtained similar results (difference was after some decimal points). I am not a frequent user of Maxima, and this problem has significantly lowered my faith in it.  Comment By: Barton Willis (willisbl) Date: 20110628 05:23 Message: I think the problem is that the default value of ratepsilon is too small; try this: (also do this same with ratepsilon : 1.0e8) (%i1) load(hypergeometric)$ (%i2) ratepsilon : 1.0e18$ (%i3) m : matrix([2621.4397,7823.3599,1111.2726],[7823.3599,23347.842,3316.4543],[1111.2726,3316.4543,471.08722])$ (%i4) first(eigenvalues(m)), ratprint : false$ (%i5) nfloat(%  conjugate(%),[],100); (%o5) [8.0266455652163197256568351091[46 digits]5913348171925384517960952b197*%i1.3377742608693866209428058515[46 digits]0985558028654230752993492b197,5.3510970434775464837712234061[46 digits]3942232114616923011973968b197*%i1.3377742608693866209428058515[46 digits]0985558028654230752993492b197,1.9934389902195135071021405630[46 digits]0374693317196116973450023b205*%i2.6755485217387732418856117030[46 digits]1971116057308461505986984b197] See also http://en.wikipedia.org/wiki/Casus_irreducibilis  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3337674&group_id=4933 
From: SourceForge.net <noreply@so...>  20120921 05:47:59

Bugs item #3435971, was opened at 20111110 00:02 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3435971&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: Closed Resolution: Invalid Priority: 5 Private: No Submitted By: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: eigenvectors produces wrong results Initial Comment: [vals,vec]:eigenvectors(matrix([0.2273,0.0852],[0.193,0.1794])); rat: replaced 0.0164436 by 1134/68963 = 0.01644360019141 rat: replaced 0.1794 by 897/5000 = 0.1794 rat: replaced 0.2273 by 2273/10000 = 0.2273 rat: replaced 0.01644360019141 by 1091/66348 = 0.01644360040996 rat: replaced 0.1794 by 897/5000 = 0.1794 rat: replaced 0.2273 by 2273/10000 = 0.2273 rat: replaced 1.2057635497678905E12 by 1/829350000000 = 1.2057635497678905E12 rat: replaced 2.1637741760636258E+11 by 216377417606/1 = 2.16377417606E+11 rat: replaced 2.1637741760636258E+11 by 216377417606/1 = 2.16377417606E+11 (%o40) [[[0.07289999842889,0.33380000157111],[1,1]],[[[1,1.812206572769953]],[[1, 1.25]]]] However, the eigenvectors should be [0.62166748, 0.78328126] and [0.46864735,0.88338534].  Comment By: Raymond Toy (rtoy) Date: 20120816 09:06 Message: The eigenvalues computed by maxima are correct. A : matrix([0.2273,0.0852],[0.193,0.1794]); A . vec[1][1]  vals[1][1] * vec[1][1] > matrix([0.0],[0.0]) A . vec[2][1]  vals[1][2] * vec[2][1] > matrix([0.2609],[.4728046948356808]) You were probably expecting the eigenvectors to be normalized to unit length. Eigenvectors are unique only up to a scale factor.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3435971&group_id=4933 
From: SourceForge.net <noreply@so...>  20120921 05:47:58

Bugs item #3479091, was opened at 20120124 10:55 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3479091&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: Invalid Priority: 5 Private: No Submitted By: Ted Woollett (woollett) Assigned to: Nobody/Anonymous (nobody) Summary: realpart(1/e) # 1/realpart(e) case Initial Comment: Example in which e = sqrt ( sin (x) ) using gcl (%i1) fpprintprec:8$ (%i2) r1(x):= realpart(1/sqrt(sin(x)))$ (%i3) map('r1,[2,5,8,11,13]); (%o3) [1/sqrt(sin(2)),0,1/sqrt(sin(8)),0,1/sqrt(sin(13))] (%i4) float(%); (%o4) [1.0486897,0.0,1.0053637,0.0,1.5427268] (%i5) r2(x) := 1/realpart(sqrt(sin(x)))$ (%i6) map('r2,[2,5]); expt: undefined: 0 to a negative exponent. #0: r2(x=5)  an error. To debug this try: debugmode(true); (%i7) build_info()$ Maxima version: 5.26.0 Maxima build date: 22:48 1/15/2012 Host type: i686pcmingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8  Comment By: Raymond Toy (rtoy) Date: 20120816 09:10 Message: What exactly is the issue? I expect 1/realpart(e) to give a divide by zero error for x=5 since sqrt(sin(5)) is purely imaginary.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3479091&group_id=4933 
From: SourceForge.net <noreply@so...>  20120921 05:47:57

Bugs item #3525906, was opened at 20120511 10:25 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3525906&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: ivan antonovich (ognelis) Assigned to: Nobody/Anonymous (nobody) Summary: bug in integrate(x*exp(a*x^2+b*x),x,X_0,inf) Initial Comment: Maxima 5.27.0 http://maxima.sourceforge.net using Lisp CLISP 2.48 (20090728) 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) display2d: false; (%o1) false (%i2) assume(a>0); (%o2) [a > 0] (%i3) expr: exp(a*x^2+b*x)*x; (%o3) x*%e^(b*xa*x^2) (%i4) res1:integrate(expr,x,X_0,inf)$ Is 2*a*X_0b positive, negative, or zero? positive; #============================================ #Let's do the integration once again: #============================================ (%i5) res1new:integrate(expr,x,X_0,inf)$ Is 2*a*X_0b positive, negative, or zero? positive; #============================================ #Let's see the results #============================================ #++++++++first++++++++++++ (%i6) factor(ratsimp(res1)); (%o6) %e^(b^2/(4*a))*(2*gamma_incomplete(1,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a)) *a*abs(2*a*X_0b) +2*gamma_incomplete(1/2, (4*a^2*X_0^24*a*b*X_0+b^2)/(4*a)) *a^(3/2)*b*X_0 gamma_incomplete(1/2,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a)) *sqrt(a)*b^2) /(4*a^2*abs(2*a*X_0b)) #++++++++second++++++++++++ (%i7) factor(ratsimp(res1new)); (%o7) (gamma_incomplete(1/2,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a))*sqrt(a)*b +2*gamma_incomplete(1,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a))*a) *%e^(b^2/(4*a)) /(4*a^2) #============================================ #The results are different!!! #============================================ #============================================ #Let's continue:(now 2 a X_0 b is negative) #============================================ (%i8) res2:integrate(expr,x,X_0,inf)$ Is 2*a*X_0b positive, negative, or zero? negative; (%i9) factor(ratsimp(res2)); (%o9) (gamma_incomplete(1/2,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a))*b 2*sqrt(%pi)*b 2*gamma_incomplete(1,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a))*sqrt(a)) *%e^(b^2/(4*a)) /(4*a^(3/2)) #============================================ #Let'see the difference between expressions with positive and negative 2 a X_0 b #============================================ (%i10) factor(ratsimp(res2res1new)); (%o10) (gamma_incomplete(1/2,(4*a^2*X_0^24*a*b*X_0+b^2)/(4*a))sqrt(%pi)) *b*%e^(b^2/(4*a)) /(2*a^(3/2)) #============================================ #But the results must be the same. #============================================ P.S. integrate(x^n*exp(a*x^2+b*x),x,X_0,inf) can be programmed as diff(integrate(exp(a*x^2+b*x),x,X_0,inf),b,n). To verify this one needs to change the order of integration and differentiation. The result can be obtained in terms of the error function. There is only one conditition  a>0.  Comment By: Raymond Toy (rtoy) Date: 20120815 09:38 Message: In maxima 5.28post, the two integrals are identical for both 2*a*X_0b positive and negative. Marking as pending/worksforme  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3525906&group_id=4933 
From: SourceForge.net <noreply@so...>  20120921 05:47:56

Bugs item #3533747, was opened at 20120608 20:27 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3533747&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: Closed Resolution: Works For Me Priority: 5 Private: No Submitted By: caposar () Assigned to: Nobody/Anonymous (nobody) Summary: plot2d: expression evaluates to nonnumeric value somewhere Initial Comment: Revisando los videos de ejemplo de Javier Arantegui "13. Animaciones (1a parte).mp4" en wxMaxim 12.04.0 me sale este error: plot2d: expression evaluates to nonnumeric value somewhere build_info("5.27.0","20120508 11:27:57","i686pcmingw32","GNU Common Lisp (GCL)","GCL 2.6.8")  Comment By: Raymond Toy (rtoy) Date: 20120815 09:32 Message: I don't know what with_slider does, but plot2d(subst(tau=v, f(t)), [t,0,4]) for various values of tau between .1 and 1 work fine. Marking as pending/worksforme.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3533747&group_id=4933 
From: SourceForge.net <noreply@so...>  20120921 05:47:55

Bugs item #3539220, was opened at 20120630 18:06 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3539220&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: Closed Resolution: Wont Fix Priority: 5 Private: No Submitted By: MukundNaik (mukundnaik) Assigned to: Nobody/Anonymous (nobody) Summary: Romberg Integral yields zero value Initial Comment: build_info("5.27.0","20120508 11:27:57","i686pcmingw32","GNU Common Lisp (GCL)","GCL 2.6.8") m1(x):=(1cos(%pi*x))*(1cos(%pi*x*2/3))*(1cos(%pi*x*2/5))*(1cos(%pi*x*2/7))/16; wxplot2d([m1(x)], [x,32,36], [y,0,0.6], [gnuplot_preamble, "set grid;"])$ romberg(m1(x), x, 32, 36); yields Zero. The plot above shows that it is not zero everywhere in the region 3236. In fact, romberg(m1(x), x, 32, 33) = 0.051305498551289 romberg(m1(x), x, 33, 34)=0.016230140453864 romberg(m1(x), x, 34, 35)=0.0014305334087985 romberg(m1(x), x, 35, 36)=2.3489482722215246*10^4.  Comment By: Raymond Toy (rtoy) Date: 20120815 09:27 Message: 1. Update summary to reflect this is a romberg issue, not symbolic integration. 2. Change category to floatingpoint 3. Changed visibility to nonprivate. Try changing some of the variables that control romberg integration. Changing rombergmin to 1 gives romberg(m1(x),x,32,36) > 0.0692, which is very close to the value of the symbolic integral.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3539220&group_id=4933 
From: SourceForge.net <noreply@so...>  20120920 19:22:17

Bugs item #3570214, was opened at 20120920 12:22 Message generated for change (Tracker Item Submitted) made by manuelgonha You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570214&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  Plotting Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Manuel GonzálezHidalgo (manuelgonha) Assigned to: Nobody/Anonymous (nobody) Summary: polar graphs are wrong Initial Comment:  build_info("5.28.02","20120827 23:16:48","i686pcmingw32","GNU Common Lisp (GCL)","GCL 2.6.8")  If I put in the console wxplot2d(t, [t,0,4*%pi],[plot_format, gnuplot], [gnuplot_preamble, "set polar; set trange [0:4*pi]; set rrange [0:4*pi];"])$ we obtain a incorrect plot of the function as we can see in the attached image. If I put in the gnuplot: set polar; set rrange [0:4*pi]; set trange [0:4*pi]; plot t; I obtain a coorect plot of the function. But in any case, that is for any function in polar coordinates, in Maxima the graph of the function isn't correct. Thank you very much for your attention, and waiting for a solution to this problem, says goodbye with a greeting and best wishes Manuel González Prof. University of the Balearic Islands  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570214&group_id=4933 
From: SourceForge.net <noreply@so...>  20120920 19:05:33

Bugs item #3570209, was opened at 20120920 12:05 Message generated for change (Tracker Item Submitted) made by manuelgonha You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570209&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  Plotting Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Manuel GonzálezHidalgo (manuelgonha) Assigned to: Nobody/Anonymous (nobody) Summary: polar graphs are wrong Initial Comment:  build_info("5.28.02","20120827 23:16:48","i686pcmingw32","GNU Common Lisp (GCL)","GCL 2.6.8")  If I put in the console wxplot2d(t, [t,0,4*%pi],[plot_format, gnuplot], [gnuplot_preamble, "set polar; set trange [0:4*pi]; set rrange [0:4*pi];"])$ we obtain a incorrect plot of the function as we can see in the attached image. If I put in the gnuplot: set polar; set rrange [0:4*pi]; set trange [0:4*pi]; plot t; I obtain a coorect plot of the function. But in any case, that is for any function in polar coordinates, in Maxima the graph of the function isn't correct. Thank you very much for your attention, and waiting for a solution to this problem, says goodbye with a greeting and best wishes Manuel González Prof. University of the Balearic Islands  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570209&group_id=4933 
From: SourceForge.net <noreply@so...>  20120920 14:36:51

Bugs item #3570112, was opened at 20120920 07:36 Message generated for change (Tracker Item Submitted) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570112&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: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: partfrac might need to gfactor ? Initial Comment: (%i35) partfrac(1/((z+2*%i*w)*(z^2+4*w^2)),z); expt: undefined: 0 to a negative exponent.  an error. To debug this try: debugmode(true); A workaround is to map gfactor on the the expression before computing the pfd: (%i41) partfrac(scanmap('gfactor, 1/((z+2*%i*w)*(z^2+4*w^2))),z); (%o41) 1/(16*w^2*(z+2*%i*w))+%i/(4*w*(z+2*%i*w)^2)1/(16*w^2*(z2*%i*w))  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3570112&group_id=4933 
From: SourceForge.net <noreply@so...>  20120920 07:40:01

Bugs item #3569562, was opened at 20120919 09:22 Message generated for change (Comment added) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&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: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: taylor(elliptic_kc(m), m, 0, 1) fails Initial Comment: Maxima cannot compute the Taylor expansion of elliptic_kc(m) for m = 0. Expansions about other points is ok. I think the problem is that diff(elliptic_kc(m),m) at m = 0 produces a form that is indeterminate at 0. Maxima is unable to evaluate the limit as m approaches 0.  Comment By: Aleksas (aleksasd) Date: 20120920 00:40 Message: From maxima help: Function: elliptic_kc (m) The complete elliptic integral of the first kind, defined as integrate(1/sqrt(1  m*sin(x)^2), x, 0, %pi/2) Example 1 Compute taylor(elliptic_kc(m), m, 0, 1) (%i1) taylor(integrate(1/sqrt(1m*sin(x)^2),x,0,%pi/2),m,0,1); (%o1)/T/ %pi/2+((at(integrate(sin(x)^2,x,0,%pi/2),m=0))*m)/2+... (%i2) ev(%,integrate); (%o2)/R/ (%pi*m+4*%pi)/8 (%i3) taylor(%,m,0,1); (%o3)/T/ %pi/2+(%pi*m)/8+... Example 2 Compute taylor(elliptic_kc(m), m, 0, 3) (%i4) taylor(integrate(1/sqrt(1m*sin(x)^2),x,0,%pi/2),m,0,3)$ (%i5) ev(%,nouns)$ (%i6) taylor(%,m,0,3); (%o6)/T/ %pi/2+(%pi*m)/8+(9*%pi*m^2)/128+(25*%pi*m^3)/512+... best Aleksas D  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&group_id=4933 
From: SourceForge.net <noreply@so...>  20120919 16:22:15

Bugs item #3569562, was opened at 20120919 09:22 Message generated for change (Tracker Item Submitted) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&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: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: taylor(elliptic_kc(m), m, 0, 1) fails Initial Comment: Maxima cannot compute the Taylor expansion of elliptic_kc(m) for m = 0. Expansions about other points is ok. I think the problem is that diff(elliptic_kc(m),m) at m = 0 produces a form that is indeterminate at 0. Maxima is unable to evaluate the limit as m approaches 0.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&group_id=4933 
From: SourceForge.net <noreply@so...>  20120910 15:29:46

Bugs item #3566334, was opened at 20120910 08:29 Message generated for change (Settings changed) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3566334&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: gcfactor fail Initial Comment: Apparently gcfactor requires an *explicit* Gaussian integer argument of the form AAA + BBB * %i, not an expression equivalent to that. This should not be necessary. But at the very least, it should be documented. foo: 74*%i11*(%i+1)^2*(%i+2)*(2*%i+1) $ <<< simplified but unexpanded number gcfactor(foo) => 74*%i+110 <<< expanded but not factored gcfactor(gcfactor(foo)) == gcfactor(expand(foo)) == gcfactor(ratsimp(foo)) => %i*(1+%i)^3*(2+3*%i)^3 <<< correct factorization  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3566334&group_id=4933 
From: SourceForge.net <noreply@so...>  20120910 15:29:11

Bugs item #3566334, was opened at 20120910 08:29 Message generated for change (Tracker Item Submitted) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3566334&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: gcfactor fail Initial Comment: Apparently gcfactor requires an *explicit* Gaussian integer argument of the form AAA + BBB * %i, not an expression equivalent to that. This should not be necessary. But at the very least, it should be documented. foo: 74*%i11*(%i+1)^2*(%i+2)*(2*%i+1) $ <<< simplified but unexpanded number gcfactor(foo) => 74*%i+110 <<< expanded but not factored gcfactor(gcfactor(foo)) == gcfactor(expand(foo)) == gcfactor(ratsimp(foo)) => %i*(1+%i)^3*(2+3*%i)^3 <<< correct factorization  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3566334&group_id=4933 
From: SourceForge.net <noreply@so...>  20120909 19:48:11

Bugs item #3566123, was opened at 20120909 12:48 Message generated for change (Tracker Item Submitted) made by woollett You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3566123&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: domain=complex integrate error Initial Comment: domain = real: no problem: (%i2) domain; (%o2) real (%i3) foo :((%pi)/2+acos(y/l)+acos(x/l))/(2*%pi)$ (%i4) assume(l>0,x>0, x<l); (%o4) [l > 0,x > 0,l > x] (%i5) integrate(integrate(foo,y,0,sqrt(l^2  x^2)),x,0,l); (%o5) l^2/(4*%pi) domain = complex: integrate error return: (%i6) domain:complex; (%o6) complex (%i7) integrate(integrate(foo,y,0,sqrt(l^2  x^2)),x,0,l); expt: undefined: 0 to a negative exponent.  an error. To debug this try: debugmode(true); (%i8) build_info(); (%o8) ?%build_info("5.28.02","20120827 23:16:48","i686pcmingw32", "GNU Common Lisp (GCL)","GCL 2.6.8") (%i9) display2d:true; (%o9) true (%i10) build_info(); (%o10) Maxima version: "5.28.02" Maxima build date: "20120827 23:16:48" 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=3566123&group_id=4933 
From: SourceForge.net <noreply@so...>  20120908 04:50:21

Bugs item #3565710, was opened at 20120907 21:50 Message generated for change (Tracker Item Submitted) made by antonvoropaev You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3565710&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: Anton Voropaev (antonvoropaev) Assigned to: Nobody/Anonymous (nobody) Summary: Maxima aborts (radcan, simp=false) Initial Comment: INPUT build_info(); radcan((6^(log(12)/log(6))+1)^(1/2)), simp=false; INPUT AND OUTPUT (%i1) build_info(); (%o1) Maxima version: "5.28.02" Maxima build date: "20120827 23:16:48" Host type: "i686pcmingw32" Lisp implementation type: "GNU Common Lisp (GCL)" Lisp implementation version: "GCL 2.6.8" (%i2) radcan((6^(log(12)/log(6))+1)^(1/2)), simp=false; <ABORT>  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3565710&group_id=4933 
From: SourceForge.net <noreply@so...>  20120903 21:17:10

Bugs item #3564492, was opened at 20120903 14:17 Message generated for change (Tracker Item Submitted) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3564492&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: carg(1),numer => %pi Initial Comment: carg(1),numer => %pi Should be float 3.14. Similarly for polarform(1),numer  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3564492&group_id=4933 
From: SourceForge.net <noreply@so...>  20120901 22:24:26

Bugs item #3558096, was opened at 20120815 17:46 Message generated for change (Comment added) made by jyoberle You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3558096&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: None Priority: 5 Private: No Submitted By: JeanYves (jyoberle) Assigned to: Nobody/Anonymous (nobody) Summary: to_poly_solve gives a wrong solution for cos(x)=sin(3x) Initial Comment: Hi, When doing: load(to_poly_solve); algexact:true; to_poly_solve(cos(x)sin(3*x),x); I get: to_poly_solve: to_poly_solver.mac is obsolete; I'm loading to_poly_solve.mac instead. %union([x=(4*%pi*%z0+%pi)/4],[x=(4*%pi*%z1+%pi)/8]) But I think that the first solution should be (based on hand solving): (4*%pi*%z0+%pi)/4 (no minus sign). For example: if we consider %z0 = 0 in the to_poly_solve solution, we get x=%pi/4 which is not a solution of the equation cos(x)sin(3*x). On the other hand, if we set %z0 = 0 in the hand found solution, we get x=%pi/4 which is a solution. The build_info is: build_info("5.27.0","20120508 11:27:57","i686pcmingw32","GNU Common Lisp (GCL)","GCL 2.6.8"). Best regards, JeanYves  >Comment By: JeanYves (jyoberle) Date: 20120901 15:24 Message: Aleksasd, I agree with you regarding your calculations. However, when I apply the function "trigsolve" to the equation, I get: trigsolve(cos(x)sin(3*x),%pi,%pi); {(7*%pi)/8,(3*%pi)/8,%pi/8,(5*%pi)/8} Several solutions are missing (e.g. %pi/4) because of the bug I highlighted. Best regards, JeanYves  Comment By: Aleksas (aleksasd) Date: 20120831 08:03 Message: To finding all solutions of trigonometric equation eq from interval [a, b] we define function "trigsolve": (%i1) trigsolve(eq,a,b):=block([s,i,ats,algebraic], algebraic:true, to_poly_solve([eq], [x],'simpfuncs = ['rootscontract,'expand,'radcan,'nicedummies]), s:makelist(rhs(part(%%,k)[1]),k,1,length(%%)), ats:[], for i:1 thru length(s) do (makelist(ev(s[i],%z0=k),k,10,10), ats:append(ats,%%)), sublist(ats,lambda([e],e>=a and e<=b and float(ev(abs(lhs(eq)rhs(eq)),x=e))<ratepsilon)), sort(%%), setify(%%) )$ Example: solve cos(x)sin(3*x)=0 (%i2) eq:cos(x)sin(3*x)=0$ (%i3) cos(x)cos(y)=2*sin(1/2*x+1/2*y)*sin(1/2*x1/2*y)$ (%i4) subst(y=3*x%pi/2,%),expand; (%o4) cos(x)sin(3*x)=2*sin(x%pi/4)*sin(2*x%pi/4) (%i5) eq1:sin(x%pi/4)=0$ (%i6) eq2:sin(2*x%pi/4)=0$ (%i7) S1:trigsolve(eq1,%pi,%pi); to_poly_solve: to_poly_solver.mac is obsolete; I'm loading to_poly_solve.mac instead. Loading maximagrobner $Revision: 1.6 $ $Date: 20090602 07:49:49 $ (%o7) {(3*%pi)/4,%pi/4} (%i8) S2:trigsolve(eq2,%pi,%pi); (%o8) {(7*%pi)/8,(3*%pi)/8,%pi/8,(5*%pi)/8} (%i9) S:union(S1,S2); (%o9) {(7*%pi)/8,(3*%pi)/4,(3*%pi)/8,%pi/8,%pi/4,(5*%pi)/8} (%i10) float(%), numer; (%o10) {2.748893571891069,2.356194490192345,1.178097245096172,0.39269908169872,0.78539816339745,1.963495408493621} Answer: x=a+2*%pi*k, where a  any from S, k  any integer (%i11) plot2d([cos(x)sin(3*x)], [x,%pi,%pi])$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3558096&group_id=4933 