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}
(20) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 




1
(2) 
2
(1) 
3
(1) 
4

5

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

11
(2) 
12
(4) 
13
(4) 
14

15
(3) 
16
(3) 
17

18

19

20

21
(2) 
22

23
(2) 
24
(3) 
25

26

27
(4) 
28
(2) 




From: SourceForge.net <noreply@so...>  20060228 12:12:44

Bugs item #1440286, was opened at 20060228 06:12 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1440286&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: Documentation Group: None Status: Open Resolution: None Priority: 1 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: documentation for 'listofvars' Initial Comment: The user documentation for 'listofvars' doesn't mention the option variable 'listdummyvars.' It ought to: (%i143) sum(i^2,i,1,n); (%o143) sum(i^2,i,1,n) (%i144) listofvars(%o143), listdummyvars : true; (%o144) [i,n] (%i145) listofvars(%o143), listdummyvars : false; (%o145) [n] Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1440286&group_id=4933 
From: SourceForge.net <noreply@so...>  20060228 02:46:31

Bugs item #1440069, was opened at 20060227 20:46 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1440069&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 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: euler numbers & zerobern Initial Comment: The option variable 'zerobern' changes the way euler evaluates. This isn't mentioned in the user documentation. Maybe it's not intended for zerobern to make any difference, or maybe it's a documentation error. I don't know. (%i1) euler(3),zerobern : true; (%o1) 0 (%i2) euler(3),zerobern : false; (%o2) 61 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1440069&group_id=4933 
From: SourceForge.net <noreply@so...>  20060227 18:01:41

Bugs item #541030, was opened at 20020408 10:21 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=541030&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: 9 Submitted By: David Billinghurst (billingd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sqrt(x+1/x2),x,0,1) wrong Initial Comment: The definite integral integrate(sqrt(x+1/x2),x,0,1) => 4/3 with cvs maxima/gcl under windows. The answer should be 4/3 according to Michael Wester, "A Review of CAS Mathematical Capabilities", 15 April 1995. (Problem 84). It definitely should be positive.  >Comment By: Raymond Toy (rtoy) Date: 20060227 13:01 Message: Logged In: YES user_id=28849 Fixed as suggested. The indefinite integral returns the integral, but I think that's ok. Perhaps maxima should ask for the sign of x1?  Comment By: Raymond Toy (rtoy) Date: 20060224 16:43 Message: Logged In: YES user_id=28849 I think the cause of this is in intform in sin.lisp. The very last cond clause sets $radexpand to '$all when doing maximasubstitute. This causes sqrt((x1)^2) to be come (x1), which is not right. If we set $radexpand to '$true, the definite integral is evaluated to 4/3. Hurray! But the indefinite integral returns the integral. Boo!  Comment By: Pedro Fortuny Ayuso (pfortuny) Date: 20020708 13:35 Message: Logged In: YES user_id=519681 FWIW, mupad 2.5 linux gives: int(sqrt(x+1/x2),x); > int(sqrt(x+1/x2),x); (i.e. not evaluated). while int(abs(x1)/sqrt(x),x); > sign(x1)(2*x^(3/2)/32*x^(1/2)) which, when evaluated at 1 gives: 4/3*sign(0). ... ... ...  Comment By: Juan Hierro (buscaideas) Date: 20020611 20:03 Message: Logged In: YES user_id=528795 I am not sure whether this should be a real bug or just a funny consequence of the indetermination in the sign of a square root. Looking at the comment by toy@..., one may realize that the result is correct provided that the sign of the square root in the integrand is shared with that of (x1). This problems appears in much simpler funcions too. For instance, let us say one has f(x):=sqrt(x*x2*x+1); Then, the sin algorithm performs the integration with sqrt always positive, while the risch algorithm works with sqrt sharing the sign of (x1). factor(diff(integrate(f(x),x),x)); ==> sqrt(x*x2*x+1) factor(diff(risch(f(x),x),x)); ==> (x1) Is it wrong this behaviour? Should there be any way to specify which branch to employ when handling multivaluated functions? f(x) in this last situation may be either abs(x1), or (x1), or abs(x1), or (1x). sin seems to work with abs(x1) and risch with (x1). In the same way, sqrt(x+1/x2) is either abs(x1)/sqrt(x) or (x1)/sqrt(x) where both signs are admisible in sqrt(x). In this case, both sin and risch algorithms seem to work with (x1)/sqrt(x).  Comment By: David Billinghurst (billingd) Date: 20020409 23:23 Message: Logged In: YES user_id=365569 toy@... (toy@...) wrote: I guess maxima gets this wrong because it says: (C1) integrate(sqrt(x+1/x2),x); 3/2 2 x  6 SQRT(x) (D1)  3 which is only true if x1 is positive. For some reason it has assumed x1 is positive somewhere during integration. Yet another integration bug. Ray  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=541030&group_id=4933 
From: SourceForge.net <noreply@so...>  20060227 17:58:13

Bugs item #694147, was opened at 20030226 22:42 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=694147&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Wrong integral: integrate(sqrt(x+1/x2),x,0,1) Initial Comment: Maxima gives 4/3 instead of 4/3 for this. An old commercial version from around 1996 gives the right answer.  >Comment By: Raymond Toy (rtoy) Date: 20060227 12:58 Message: Logged In: YES user_id=28849 Duplicate of 541030 integrate(sqrt(x+1/x2),x,0,1) wrong  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=694147&group_id=4933 
From: SourceForge.net <noreply@so...>  20060227 11:38:36

Bugs item #1439566, was opened at 20060227 05:38 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1439566&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: Documentation Group: None Status: Open Resolution: None Priority: 1 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: zerobern & bernpoly Initial Comment: The user documentation for 'bernpoly' doesn't mention that the value of 'zerobern' changes the definition of the Bernoulli polynomial. Actually the documentation refers to "the Bernoulli polynomial," so it's surprising that there are two definitions: (%i7) bernpoly(x,3), zerobern : false; (%o7) x^3(3*x^2)/2+x/21/30 (%i8) bernpoly(x,3), zerobern : true; (%o8) x^3(3*x^2)/2+x/2 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1439566&group_id=4933 
From: SourceForge.net <noreply@so...>  20060227 11:30:43

Bugs item #1439559, was opened at 20060227 05:30 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1439559&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: 1 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: function burn is broken Initial Comment: The function 'burn' is broken: (%i1) load("bffac")$ (%i2) burn(0); argument value `0' to args was not a list (%i3) burn(1); (%o3) 1/2 < OK Until this is fixed, use the function 'bern.' Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1439559&group_id=4933 
From: SourceForge.net <noreply@so...>  20060224 21:43:36

Bugs item #541030, was opened at 20020408 10:21 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=541030&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: 9 Submitted By: David Billinghurst (billingd) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(sqrt(x+1/x2),x,0,1) wrong Initial Comment: The definite integral integrate(sqrt(x+1/x2),x,0,1) => 4/3 with cvs maxima/gcl under windows. The answer should be 4/3 according to Michael Wester, "A Review of CAS Mathematical Capabilities", 15 April 1995. (Problem 84). It definitely should be positive.  >Comment By: Raymond Toy (rtoy) Date: 20060224 16:43 Message: Logged In: YES user_id=28849 I think the cause of this is in intform in sin.lisp. The very last cond clause sets $radexpand to '$all when doing maximasubstitute. This causes sqrt((x1)^2) to be come (x1), which is not right. If we set $radexpand to '$true, the definite integral is evaluated to 4/3. Hurray! But the indefinite integral returns the integral. Boo!  Comment By: Pedro Fortuny Ayuso (pfortuny) Date: 20020708 13:35 Message: Logged In: YES user_id=519681 FWIW, mupad 2.5 linux gives: int(sqrt(x+1/x2),x); > int(sqrt(x+1/x2),x); (i.e. not evaluated). while int(abs(x1)/sqrt(x),x); > sign(x1)(2*x^(3/2)/32*x^(1/2)) which, when evaluated at 1 gives: 4/3*sign(0). ... ... ...  Comment By: Juan Hierro (buscaideas) Date: 20020611 20:03 Message: Logged In: YES user_id=528795 I am not sure whether this should be a real bug or just a funny consequence of the indetermination in the sign of a square root. Looking at the comment by toy@..., one may realize that the result is correct provided that the sign of the square root in the integrand is shared with that of (x1). This problems appears in much simpler funcions too. For instance, let us say one has f(x):=sqrt(x*x2*x+1); Then, the sin algorithm performs the integration with sqrt always positive, while the risch algorithm works with sqrt sharing the sign of (x1). factor(diff(integrate(f(x),x),x)); ==> sqrt(x*x2*x+1) factor(diff(risch(f(x),x),x)); ==> (x1) Is it wrong this behaviour? Should there be any way to specify which branch to employ when handling multivaluated functions? f(x) in this last situation may be either abs(x1), or (x1), or abs(x1), or (1x). sin seems to work with abs(x1) and risch with (x1). In the same way, sqrt(x+1/x2) is either abs(x1)/sqrt(x) or (x1)/sqrt(x) where both signs are admisible in sqrt(x). In this case, both sin and risch algorithms seem to work with (x1)/sqrt(x).  Comment By: David Billinghurst (billingd) Date: 20020409 23:23 Message: Logged In: YES user_id=365569 toy@... (toy@...) wrote: I guess maxima gets this wrong because it says: (C1) integrate(sqrt(x+1/x2),x); 3/2 2 x  6 SQRT(x) (D1)  3 which is only true if x1 is positive. For some reason it has assumed x1 is positive somewhere during integration. Yet another integration bug. Ray  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=541030&group_id=4933 
From: SourceForge.net <noreply@so...>  20060224 16:48:30

Bugs item #928282, was opened at 20040402 09:38 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=928282&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(r^2/(r^2+1)^(3/2) => GCD bug (use spmod instead) Initial Comment: expression: integrate(r^2*(1+r^2)^(03/2),r,0,R); R:=positive Error message: Quotient by a polynomial of higher degree  >Comment By: Raymond Toy (rtoy) Date: 20060224 11:48 Message: Logged In: YES user_id=28849 This doesn't seem to happen with maxima 5.9.2.19cvs. Closing this bug, but perhaps we should add another bug for gcd.  Comment By: Stavros Macrakis (macrakis) Date: 20040408 19:00 Message: Logged In: YES user_id=588346 Sorry for this bug. The problem is in the default GCD routine. Please set gcd: 'spmod$ and try again.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=928282&group_id=4933 
From: SourceForge.net <noreply@so...>  20060224 14:37:52

Bugs item #846112, was opened at 20031120 14:58 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=846112&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: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: scsimp(x,x) / FIX Initial Comment: Consider (C1) scsimp(x+1,x5=0); (D1) 6 Given no rhs, Maxima defaults it to 0 (C2) scsimp(x+1,x5); (D2) 6 The default doesn't work with an atom; this is okay (C3) scsimp(x+1,x=0); (D3) 1 but why an error for this case? (C4) scsimp(x+1,x); Error: $x is not of type LIST. Fast links are on: do (si::usefastlinks nil) for debugging Error signalled by MACSYMATOPLEVEL. Broken at MACSYMATOPLEVEL. Type :H for Help. MAXIMA>> Here is a possible fix (DEFMFUN $SCSIMP N (DO ((I N (f1 I)) (ZRS)) ((= 1 I) (SCS (ARG 1) ZRS)) (setq zrs (cons (meqhk (arg i)) zrs)))) With this definition for scsimp, (C7) scsimp(x+1,x); (D7) 1 (C8) Barton  Comment By: Barton Willis (willisbl) Date: 20060224 08:37 Message: Logged In: YES user_id=895922 The reported bug is not present in the current cvs version of Maxima. Thank you for your report. If you see this bug in a later version of Maxima, please submit a new bug report.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=846112&group_id=4933 
From: SourceForge.net <noreply@so...>  20060223 17:37:49

Bugs item #1030837, was opened at 20040919 16:25 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1030837&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: wrong results for integrate(sqrt(cos(x)+1)); Initial Comment: integrate(sqrt(cos(x)+1),x); ==> 2 SQRT(2) SIN(x)  2 SIN (x) (COS(x) + 1) SQRT( + 1) 2 (COS(x) + 1) This is wrong!  >Comment By: Raymond Toy (rtoy) Date: 20060223 12:37 Message: Logged In: YES user_id=28849 Closing this report. The derivative matches the integrand, after some manipulations.  Comment By: Raymond Toy (rtoy) Date: 20041101 17:32 Message: Logged In: YES user_id=28849 Why is this wrong? If I differentiate the result, I get: trigsimp(diff(%,x)); (%o35) (sqrt(2)*sin(x)^22*sqrt(2)*cos(x)2*sqrt(2)) /((cos(x)+1)*sqrt(2*cos(x)+2)) (%i36) factor(%); (%o36) (sin(x)^22*cos(x)2)/(cos(x)+1)^(3/2) (%i37) ev(%,sin(x)^2=1cos(x)^2); (%o37) (cos(x)^22*cos(x)1)/(cos(x)+1)^(3/2) (%i38) factor(%); (%o38) sqrt(cos(x)+1)  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1030837&group_id=4933 
From: SourceForge.net <noreply@so...>  20060223 11:00:08

Bugs item #1435602, was opened at 20060220 20:27 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1435602&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: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: string length function Initial Comment: There is no string_length function in Maxima 5.9.1: In Macsyma we count the length of a string by: string_length("aaaaaaaaaa"); 10 This is awkward when one wishes to manipulate strings.  >Comment By: Barton Willis (willisbl) Date: 20060223 05:00 Message: Logged In: YES user_id=895922 In November 2005, Volker van Nek contributed a package for working with strings in Maxima. You can download his code from the Maxima CVS. Look for the directory /share/contrib/stringproc. Or you can wait for the next release of Maxima. (%i4) load("stringproc")$ (%i5) slength("aaaaaaaaaa"); (%o5) 10 Barton  Comment By: Barton Willis (willisbl) Date: 20060223 05:00 Message: Logged In: YES user_id=895922 The reported bug is not present in the current cvs version of Maxima. Thank you for your report. If you see this bug in a later version of Maxima, please submit a new bug report.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1435602&group_id=4933 
From: SourceForge.net <noreply@so...>  20060221 02:27:19

Bugs item #1435602, was opened at 20060220 18:27 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1435602&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: string length function Initial Comment: There is no string_length function in Maxima 5.9.1: In Macsyma we count the length of a string by: string_length("aaaaaaaaaa"); 10 This is awkward when one wishes to manipulate strings.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1435602&group_id=4933 
From: SourceForge.net <noreply@so...>  20060221 02:24:41

Bugs item #1435600, was opened at 20060220 18:24 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1435600&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: equal function Initial Comment: is(equal(true,true)) will return true is(equal(false,false)) will return false but is(equal(true,false) or is(equal(false,true)) will give a bug ... syntax error. Why not provide a k_delta function like in Macsymas? e.g. k_delta(true,false) returns false k_delta(true,true) returns true k_delta(1,1) returns 1 k_delta(1,0) returns 0 HuenYK v 5.9.1 Maxima Comment: there are some good features in Maxima especially on Random(n). Macsymas only give a ceiling of n = 10^8.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1435600&group_id=4933 
From: SourceForge.net <noreply@so...>  20060216 15:27:00

Bugs item #1051437, was opened at 20041021 08:21 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1051437&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 Submitted By: David Billinghurst (billingd) Assigned to: Nobody/Anonymous (nobody) Summary: Trig integral error Initial Comment: The integral of 2*COT(x)^2*COS(2*x)/(CSC(2*x)+COT(2*x)); is wrong for maxima5.9.1 (%i1) display2d:false; (%o1) FALSE (%i2) h: 2*COT(x)^2*COS(2*x)/(CSC(2*x)+COT(2*x)); (%o2) 2*COT(x)^2*COS(2*x)/(CSC(2*x)+COT(2*x)) (%i3) ih:integrate(h,x); (%o3) (2*LOG(SIN(x)^2+COS(x)^2+2*COS(x)+1) +2*LOG(SIN(x)^2+COS(x)^22*COS(x)+1) +COS(2*x)) /2 (%i4) ev(ih,x=1.0,numer)ev(ih,x=0.5,numer); (%o4) .6469013090248041 (%i5) quad_qags(h,x,0.5,1); (%o5) [.1686767378171631,3.37999776996994E 15,21,0] (%i6) h2:trigsimp(trigexpand(h)); (%o6) (4*COS(x)^32*COS(x))/SIN(x) (%i7) ih2:integrate(h2,x); (%o7) LOG(COS(x)+1)+LOG(COS(x)1)+2*COS(x)^2 (%i8) ev(ih2,x=1.0,numer)ev(ih2,x=0.5,numer); (%o8) .1686767378171636 The integral over 0.5 < x < 1.0 at %o4 differs from the numerical integral %o5 and the analytic integral of an equivalent expression %o8.  >Comment By: Raymond Toy (rtoy) Date: 20060216 10:26 Message: Logged In: YES user_id=28849 This seems to be a bug in the Risch integrator. If you trace(?rischint), you can see h is converted to exponential form and the Risch integrator returns log(exp(%i*x)+1)+log(exp(%i*x)1)+(exp(2*%i*x)4*%i*x)/4. However integrate(trigsimp(exponentialize(h)),x); (which uses the Risch integrator too!) returns 2*log(%e^(%i*x)+1)+2*log(%e^(%i*x)1)+%e^(2*%i*x)/2+%e^(2*%i*x)/22*%i*x I think this simplifies to log(2*cos(x)+2)+log(22*cos(x))+cos(2*x) Differentiating this produces something equal to h.  Comment By: David Billinghurst (billingd) Date: 20041021 08:29 Message: Logged In: YES user_id=365569 Once this is fixed, activate equation (22) in share/contrib/diffequations/tests/rtestode_murphy1.mac  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1051437&group_id=4933 
From: SourceForge.net <noreply@so...>  20060216 14:35:36

Bugs item #1290363, was opened at 20050913 14:39 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1290363&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 Submitted By: Vadim V. Zhytnikov (vvzhy) Assigned to: Nobody/Anonymous (nobody) Summary: integrate((tan(x)^2+1)/tan(x),x,%pi/6,%pi/3)  error Initial Comment: (%i1) integrate((tan(x)^2+1)/tan(x),x,%pi/6,%pi/3); `sign' called on an imaginary argument: 1/4 ( 1)  an error. Quitting. To debug this try debugmode(true); Right result is log(3).  >Comment By: Raymond Toy (rtoy) Date: 20060216 09:35 Message: Logged In: YES user_id=28849 This integral is converted to integrate((1+tan(x+%pi/6)^2)/tan(x+%pi/6),x,0,%pi/6) Eventually, rischint is called and the error comes from somewhere in rischint.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1290363&group_id=4933 
From: SourceForge.net <noreply@so...>  20060216 01:26:42

Bugs item #635357, was opened at 20021108 07:57 Message generated for change (Comment added) made by andrejv You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=635357&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 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: factor infinite loop Initial Comment: expr: (%i*a1)/(a^2+1)^2 factor(expr) runs forever but factor((%i*a1)/(a^2+1)) works fine, as does gfactor(expr)  >Comment By: Andrej Vodopivec (andrejv) Date: 20060216 02:26 Message: Logged In: YES user_id=1179910 This is a bug in the simplifier: ratsimp(expr), ratfac=true; goes to infinite loop. The infinite loop is in function lgcd1 in lesfac.lisp: aloop (cond ((setq t1 (testdivide ai c)) (setq ai t1 d1 (f1+ d1)) (go aloop))) ai=a*%i1 and c=%i. Since division by %i is multiplication by %i, maxima keeps dividing ai with c. It is not clear if maxima should never get to this loop or if testdivide should not perform the division by %i. Andrej  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=635357&group_id=4933 
From: SourceForge.net <noreply@so...>  20060215 20:23:39

Bugs item #1309432, was opened at 20050930 08:19 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1309432&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(1/cosh(a*x)^2,x,inf,inf); Initial Comment:  Maxima version: 5.9.1 Maxima build date: 16:35 2/10/2005 host type: i686pclinuxgnu lispimplementationtype: GNU Common Lisp (GCL) lispimplementationversion: GCL 2.6.6  integrate(1/cosh(a*x)^2,x,inf,inf); Is a positive, negative, or zero? p; (%o3) 0 Correct answer is 2/a.  >Comment By: Raymond Toy (rtoy) Date: 20060215 15:23 Message: Logged In: YES user_id=28849 I think this integral fails because polelist is unable to find the roots of z^(4*a)+2*z^(2*a)+1.  Comment By: Stavros Macrakis (macrakis) Date: 20051109 14:39 Message: Logged In: YES user_id=588346 Interestingly, if you exponentialize the expression first, it gets the right answer. But still asks, unnecessarily, the sign of 'a'.  Comment By: Stavros Macrakis (macrakis) Date: 20051109 14:34 Message: Logged In: YES user_id=588346 Interestingly, if you exponentialize the expression first, it gets the right answer. But still asks, unnecessarily, the sign of 'a'.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1309432&group_id=4933 
From: SourceForge.net <noreply@so...>  20060215 19:55:09

Bugs item #1405931, was opened at 20060114 11:12 Message generated for change (Settings changed) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1405931&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(log(x)/(x^2+1)^2,x,0,inf)  wrong result Initial Comment: (%i1) integrate(log(x)/(x^2+1)^2,x,0,inf); (%o1) (3*%i*%pi^24*%pi)/32 Correct result should be %pi/4.  Maxima version: 5.9.2.15cvs Maxima build date: 17:26 1/9/2006 host type: powerpcappledarwin8.3.0 lispimplementationtype: SBCL lispimplementationversion: 0.9.8   >Comment By: Raymond Toy (rtoy) Date: 20060215 14:55 Message: Logged In: YES user_id=28849 Fixed in CVS. The problem was that we were not computing the residues of all the poles.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1405931&group_id=4933 
From: SourceForge.net <noreply@so...>  20060215 12:18:22

Bugs item #1432121, was opened at 20060215 06:18 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1432121&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 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: limit(infinity^inf) Initial Comment: limit(infinity^inf) throws Maxima into an infinite loop: .2 Enter limit [log(foo),foo,infinity] .2 Exit limit infinity .2 Enter limit [foo,foo,inf] .2 Exit limit inf .2 Enter limit [log(foo),foo,infinity] .2 Exit limit infinity .2 Enter limit [foo,foo,inf] .2 Exit limit inf Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1432121&group_id=4933 
From: SourceForge.net <noreply@so...>  20060213 18:10:57

Bugs item #1430843, was opened at 20060213 13:10 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430843&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 Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: integrate(x/(1x^2),x) and integrate(x/(x^21),x) Initial Comment: integrate(x/(1x^2),x) and integrate(x/(x^21),x) don't ask if x < 1 or x > 1. Because of this we get two different results, when they should be equivalent, depending on the range of integration.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430843&group_id=4933 
From: SourceForge.net <noreply@so...>  20060213 18:07:03

Bugs item #1374700, was opened at 20051206 13:41 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1374700&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 Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: integrate((1+tan(x)^2)/tan(x),x); Initial Comment: Nonreal result  >Comment By: Raymond Toy (rtoy) Date: 20060213 13:07 Message: Logged In: YES user_id=28849 This integral is transformed to cos(x)/sin(x)*(sin(x)^2/cos(x)^2+1). Then maxima uses the substitution y=sin(x) to get 1/y*(y^2/(1y^2)+1. However: integrate(1/y*(y^2/(1y^2)+1),y) > log(y)log(y^21)/2. But integrate(expand(1/y*(y^2/(1y^2)+1)),y) > log(y)log(1y^2)/2. The former is wrong for our integration problem; the latter would produce the desired answer.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1374700&group_id=4933 
From: SourceForge.net <noreply@so...>  20060213 11:09:54

Bugs item #1430379, was opened at 20060212 19:10 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430379&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 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: algsys & algebraic == true Initial Comment: (%o7) [b*c+a^2a,b*d+a*bb,c*d+a*cc,d^2d+b*c] (%i8) algsys(%,[a,b,c,d]); Maxima encountered a Lisp error: (%i9) algsys(%o7,[a,b,c,d]),algebraic : true; (%o9) [[a=(sqrt(14*%r7*%r8)1)/2, ...]] Would it be OK to have algsys (or maybe just resultant) set algebraic to true? Barton  >Comment By: Barton Willis (willisbl) Date: 20060213 05:09 Message: Logged In: YES user_id=895922 Proposed fix: (defmfun $resultant (a b mainvar) (prog (varlist formflag $ratfac res ans genvar $keepfloat $algebraic) (setq varlist (list mainvar) $ratfac t ans 1 $algebraic t) It seems to fix the problem: (%o9) [b*c+a^2a,b*d+a*bb,c*d+a*cc,d^2d+b*c] (%i10) algsys(%,[a,b,c,d]); (%o10) [[a=(sqrt(14*%r7*%r8)1)/2,b=%r7 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430379&group_id=4933 
From: SourceForge.net <noreply@so...>  20060213 01:10:35

Bugs item #1430379, was opened at 20060212 19:10 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430379&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 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: algsys & algebraic == true Initial Comment: (%o7) [b*c+a^2a,b*d+a*bb,c*d+a*cc,d^2d+b*c] (%i8) algsys(%,[a,b,c,d]); Maxima encountered a Lisp error: (%i9) algsys(%o7,[a,b,c,d]),algebraic : true; (%o9) [[a=(sqrt(14*%r7*%r8)1)/2, ...]] Would it be OK to have algsys (or maybe just resultant) set algebraic to true? Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430379&group_id=4933 
From: SourceForge.net <noreply@so...>  20060212 19:24:21

Bugs item #1430245, was opened at 20060212 11:24 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430245&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 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: error using ode2 Initial Comment: Maxima 5.9.2 for Windows crashes when using this statement ode2('diff(y,x,2)'diff(y,x)+y=cos(x),y,x); Paolo Ferraresi.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1430245&group_id=4933 
From: SourceForge.net <noreply@so...>  20060212 19:16:08

Bugs item #1376860, was opened at 20051208 22:53 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1376860&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: specint(gammaincomplete(v,a*t)*exp(p*t),t) seems wrong Initial Comment: (%i12) assume(v>0,p>0); (%o12) [v>0,p>0] (%i13) specint(gammaincomplete(v,a*t)*exp(p*t),t); SIMP2F1WILLCONTINUEIN (%o13) a^v*(p+a)^(v1)*gamma(v+1)*%f[2,1]([v+1,3/2],[2],p/(p+a)) Compare this with (%i14) specint(gammagreek(v,a*t)*exp(p*t),t); (%o14) a^v*p^(v1)*gamma(v+1)/((a/p+1)^v*v) This matches formula 34, p 179 in Tables of Transforms. Considering gammaincomplete = gamma(n)gammagreek, the expression for gammaincomplete seems wrong. It might still be right if the hypergeometric function simplifies, but maxima can't, and I can't think of any way to simplify it either.  >Comment By: Raymond Toy (rtoy) Date: 20060212 14:15 Message: Logged In: YES user_id=28849 The fix for transforming %w causes this return the result in terms of the associated Legendre function Q. Somewhat better, but I do not know if this is equivalent.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1376860&group_id=4933 