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}
(54) 
_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 







1
(1) 
2

3

4

5

6

7
(3) 
8

9

10
(1) 
11
(1) 
12

13

14

15

16
(1) 
17

18
(1) 
19
(1) 
20

21
(4) 
22

23
(3) 
24
(4) 
25
(4) 
26
(3) 
27
(2) 
28
(3) 
29
(1) 
30
(6) 
31
(4) 





From: SourceForge.net <noreply@so...>  20050124 18:31:45

Bugs item #1097982, was opened at 20050107 11:30 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1097982&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: limit(x/(x^(log(x)))); returns wrong answer Initial Comment: In the following, the first two output are correct, but the third one is wrong. I do not know why.  (C1) limit(x/(1+x^log(x)),x,inf); (D1) 0 (C2) limit(x/(1+x^log(x)),x,inf); (D2) 0 (C3) limit(x/(x^log(x)),x,inf); (D3) INF  >Comment By: Stavros Macrakis (macrakis) Date: 20050124 13:31 Message: Logged In: YES user_id=588346 More closeby examples which it gets right: limit(x^(alog(x)),x,inf) (a>0) => 0 OK limit(x^(alog(x)),x,inf) (a<0) => 0 OK limit(x^( log(x)),x,inf) => 0 OK limit(inf^a) (a>0) => inf OK limit(inf^a) (a<0) => 0 OK limit(x^minf) (abs(x)1>0, x>0) => 0 OK limit(x^(alog(x)),x,inf) (a=0) => nounform ???  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1097982&group_id=4933 
From: SourceForge.net <noreply@so...>  20050124 02:20:39

Bugs item #1108066, was opened at 20050123 18:20 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=1108066&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: How to turn off browser when startup? Initial Comment: I don't know why browser is shown every time I start Maxima. It seems useless for me. Can I ture it off automatically when start Maxima?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1108066&group_id=4933 
From: SourceForge.net <noreply@so...>  20050123 20:50:52

Bugs item #1107943, was opened at 20050123 12:50 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=1107943&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: compile_file doesn't escape spaces in paths Initial Comment: Using Maxima 5.9.1 with GCL under MS Windows, when trying to compile the file, SYSTEM call for gcc didn't escape spaces in paths to the file and to the includes directory (yet the default instalation path for Maxima is in "Program Files") which caused gcc report a (misleading) error. Lukasz Stafiniak <lukstafi@...>  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1107943&group_id=4933 
From: SourceForge.net <noreply@so...>  20050123 16:14:55

Bugs item #1107784, was opened at 20050123 15:40 Message generated for change (Comment added) made by franco68tn You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1107784&group_id=4933 Category: Lisp Core Group: Fix for 5.9.0 Status: Open Resolution: None Priority: 5 Submitted By: Gian Paolo Bronzetti (a_toy_lab) Assigned to: Nobody/Anonymous (nobody) Summary: Plot2/3D: If a Limit is an Integer Fraction Initial Comment: (%i14) plot2d(2*x,[x,5,5])$ is OK, while, if 5 > 5/2 (%i15) plot2d(2*x,[x,5,5/2])$ Bad Range [x,5,5/2] must be of the form [variable,min,max]  an error. Quitting. To debug this try DEBUGMODE(TRUE); The same happens using plot3d(). If 5/2 > 5.0/2 then it works OK (Amundson) Paolo  Comment By: FrancoB (franco68tn) Date: 20050123 17:14 Message: Logged In: YES user_id=1074823 plot2d(2*x,[x,5,float(5/2)])$ is also OK. Franco  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1107784&group_id=4933 
From: SourceForge.net <noreply@so...>  20050123 14:40:42

Bugs item #1107784, was opened at 20050123 15:40 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=1107784&group_id=4933 Category: Lisp Core Group: Fix for 5.9.0 Status: Open Resolution: None Priority: 5 Submitted By: Gian Paolo Bronzetti (a_toy_lab) Assigned to: Nobody/Anonymous (nobody) Summary: Plot2/3D: If a Limit is an Integer Fraction Initial Comment: (%i14) plot2d(2*x,[x,5,5])$ is OK, while, if 5 > 5/2 (%i15) plot2d(2*x,[x,5,5/2])$ Bad Range [x,5,5/2] must be of the form [variable,min,max]  an error. Quitting. To debug this try DEBUGMODE(TRUE); The same happens using plot3d(). If 5/2 > 5.0/2 then it works OK (Amundson) Paolo  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1107784&group_id=4933 
From: SourceForge.net <noreply@so...>  20050121 19:07:06

Bugs item #1106912, was opened at 20050121 14:07 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=1106912&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(x/sin(x)^2@... Initial Comment: limit(x/sin(x)^2,x,inf) => UND actually = inf  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1106912&group_id=4933 
From: SourceForge.net <noreply@so...>  20050121 19:04:29

Bugs item #593351, was opened at 20020809 22:48 Message generated for change (Settings changed) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593351&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) >Summary: limit/sin(inf)etc. should give 0, not IND Initial Comment: limit(cos(1/x)*sin(x)sin(x),x,inf) should give 0, not IND.  Comment By: Stavros Macrakis (macrakis) Date: 20020918 23:55 Message: Logged In: YES user_id=588346 Strangely enough, limit even gets this wrong if you feed it the factored form: limit ( (cos(1/x)1) * sin(x), x, inf) even though it correctly gets limit(cos(1/x)1,x,inf) => 0 and limit(sin(x),x,inf) => ind and 0*ind is always 0. On the other hand, it does get it right if you factor the exponentialized form: limit(factor(ev(...,exponentialize)),x,inf) => 0  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593351&group_id=4933 
From: SourceForge.net <noreply@so...>  20050121 19:01:59

Bugs item #593353, was opened at 20020809 22:54 Message generated for change (Settings changed) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593353&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) >Summary: Limit/log/+/^ returns 2 instead of 1 Initial Comment: expr: (log(x)^a+log(x)^b)/(log(x)^c+log(x)^d) limit(expr,x,inf) Is a positive, negative, or zero? pos; Is c positive, negative, or zero? pos; Is d positive, negative, or zero? pos; Is a an integer? yes; Is d  c positive, negative, or zero? zero; Is c  a positive, negative, or zero? zero; Is c an integer? yes; Is b positive, negative, or zero? pos; Is b an integer? yes; Is c  b positive, negative, or zero? zero; (That is, all the exponents are equal and positive  and integral, as if that mattered.) Returns the answer ** 2 **. Should be 1.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=593353&group_id=4933 
From: SourceForge.net <noreply@so...>  20050121 19:00:17

Bugs item #1102908, was opened at 20050115 08:08 Message generated for change (Settings changed) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1102908&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Daniel van Leeuwen (laeusen) Assigned to: Nobody/Anonymous (nobody) >Summary: limit/atan/exp returns complex expr with wrong principal val Initial Comment: Maxima version: 5.9.1 Maxima build date: 14:17 9/30/2004 host type: i686redhatlinuxgnu lispimplementationtype: CLISP lispimplementationversion: 2.33.2 (20040602) (built on sting.unl.edu [129.93.181.39]) (C1) limit(atan(x)/(1/exp(1)exp((1+x)^2)),x,inf,plus); %E %I LOG( 1) (D1)  2 /* Since LOG(1) equals %I*%PI, this should be the same as */ (C2) subst(%I*%PI,log(1),%th(1)); %E %PI (D2)   2 /* But doing the calculation stepbystep */ (C3) limit(atan(x),x,inf); %PI (D3)  2 (C4) limit((x+1)^2,x,inf); (D4) MINF (C5) limit(exp((x+1)^2),x,inf); (D5) 0 /* shows that the solution should be: * * exp(1)*%pi/2 * */  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1102908&group_id=4933 
From: SourceForge.net <noreply@so...>  20050119 16:47:37

Bugs item #1105366, was opened at 20050119 08:47 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=1105366&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: SOLVE gives a wrong answer Initial Comment: When I try this: kill(all)$ x(t):=exp(t)*cos(t)$ y(t):=exp(t)*sin(t)$ z(t):=exp(t)*sqrt(3)$ g(t):=[x(t),y(t),z(t)]$ d(t):=sqrt(g(t).g(t))$ d(t)=2; solve(%,t); I get the following output: [t=log(2/sqrt(3)),t=log(2/sqrt(3))] but the correct answer should be t= 0. Using TRIGSIMP or TRIGREDUCE I get te correct answer: trigreduce(d(t)=2); solve(%,t); Now I get: [t = 0] Is this a bug? Franco Buratti (Italy) bufranz@...  Maxima version: 5.9.1 Maxima build date: 7:34 9/24/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.5   You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1105366&group_id=4933 
From: SourceForge.net <noreply@so...>  20050118 21:17:12

Bugs item #1073338, was opened at 20041125 10:21 Message generated for change (Comment added) made by nobody You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1073338&group_id=4933 Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Submitted By: Robert Dodier (robert_dodier) Assigned to: Nobody/Anonymous (nobody) Summary: integrate yields incorrect result on rational function Initial Comment: "integrate" yields incorrect results on some rational functions. "Division by 0" is strange. The definite integral below is certainly greater than 0 as the integrand is positive over [0, 1]. "integrate (1/((x3)^4+1/2), x)" returns the noun form, so maybe (maybe) what happens is that the noun form is evaluated at the limits of integration and it's the same, hence 0 is the result. (Just guessing there.) Note that the difference between the two integrands is that one is 1/(something + 1), while the other is 1/(same something + 1/2).  (%i1) integrate (1/((x3)^4+1), x, 0, 1); Division by 0  an error. Quitting. To debug this try DEBUGMODE(TRUE); (%i2) integrate (1/((x3)^4+1/2), x, 0, 1); (%o2) 0 (%i3) build_info (); Maxima version: 5.9.1 Maxima build date: 21:24 9/23/2004 host type: i686pclinuxgnu lispimplementationtype: CMU Common Lisp lispimplementationversion: 19a  Same behavior observed in CVS build of 2004/11/24.  Comment By: Nobody/Anonymous (nobody) Date: 20050118 13:17 Message: Logged In: NO Two remarks: The method that Maxima uses to solve such integrals is integration by residues. Trace (residue) shows that this function is called with the correct arguments but answers zeroes for integrate (1/((x3)^4+1/2), x, 0, 1); and runs into a division by zero for the other integral. The indefinite integrals can be solved with changevar: 'Integrate(1/((x3)^4+1/2), x); changevar (%, x  3  y ,y ,x); ev (%, Integrate);  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1073338&group_id=4933 
From: SourceForge.net <noreply@so...>  20050116 21:45:30

Bugs item #1103515, was opened at 20050116 16:45 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=1103515&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: limit(atan2(x,1),x,0) wrong Initial Comment: limit(atan2(x,1),x,0,plus) => %pi OK limit(atan2(x,1),x,0,minus) => %pi OK limit(atan2(x,1),x,0) => %pi !!! should be IND  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1103515&group_id=4933 
From: SourceForge.net <noreply@so...>  20050111 13:05:18

Bugs item #1100129, was opened at 20050111 05:05 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=1100129&group_id=4933 Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: limit finds not a limit from defint, but finds it from topl Initial Comment: I found a problem Maxima hangs at, could someone take a look? I tried to solve a definite integration problem, and then I tried to track down the problem in maxima internals. I found that $limit fails on an expression which is easily solved from toplevel. Any ideas on what's happening? One more question: unless I explicilty do (trace sratsimp behaviorbydiff), :bt reporst a stop in $limit. It take quite a time to type in most of limit.lisp function names to get to the actual timeeater. Please, show me a technique to dig an actual call sequence?  Andrei Zorine Maxima 5.9.1 http://maxima.sourceforge.net Using Lisp Kyoto Common Lisp GCL 2.6.5 (aka GCL) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. This is a development version of Maxima. The function bug_report() provides bug reporting information. (%i1) (sqrt(sqrt(x^2+1)+x)sqrt(sqrt(x^2+1)x))/x/(x^2+1); 2 2 SQRT(SQRT(x + 1) + x)  SQRT(SQRT(x + 1)  x) (%o1)  2 x (x + 1) (%i2) limit(x*%,x,0); (%o2) 0 (%i3) defint(%o1,x,0,inf); Maxima encountered a Lisp error: Console interrupt. Automatically continuing. To reenable the Lisp debugger set *debuggerhook* to nil. (%i4) :lisp(trace sratsimp) (SRATSIMP) (%i4) :lisp(trace behaviorbydiff) (BEHAVIORBYDIFF) (%i5) :lisp(setq *debuggerhook* nil) NIL (%i5) defint(%o1,x,0,inf); console interrupt. Fast links are on: do (usefastlinks nil) for debugging Broken at SYSTEM::CLCSTERMINALINTERRUPT. Type :H for Help. 1 (Continue) Continues execution. 2 (Abort) Return to top level. dbl:MAXIMA>>:bt #0 CLCSTERMINALINTERRUPT {loc0=t} [ihs=17] #1 TRACECALL {tempname=nil,args=nil,cond=nil,entrycond=((system::arglist (#))),entry=nil,exi...} [ihs=16] #2 SRATSIMP {(#0=(mplus . #1=(simp)) (#2=(mtimes . #1#) (#3=# 1 2) (#4=# #5=# 1) ...) (#2#...} [ihs=15] #3 TRACECALL {tempname=nil,args=(#0=(mplus . #1=(simp)) (#2=(mtimes . #1#) (#3=# 1 2) (#4=#...} [ihs=14] #4 BEHAVIORBYDIFF {(#0=(mplus . #1=(simp)) (#2=(mtimes . #1#) 1 (#3=# # 1) ...) (#2# (#3# # 1) ...} [ihs=13] #5 $LIMIT {loc0=((mtimes . #0=(simp)) (#1=(mexpt . #0#) (#2=# 1 #) 1) (#2# (# 1 #) (#1# ...} [ihs=12] #6 LIMITNOERR {loc0=((mtimes . #0=(simp)) (#1=(mexpt . #0#) (#2=# 1 #) 1) (#2# (# 1 #) (#1# ...} [ihs=11] #7 MACSYMATOPLEVEL {inputstream=:internal,batchflag=$inf,loc2=nil,loc3=stringchar,loc4="0",loc...} [ihs=10] #8 RUN {} [ihs=6] #9 TOPLEVEL {loc0=nil,loc1=nil,loc2=nil,loc3=(lambdablock run nil ...)} [ihs=5] #10 FUNCALL {loc0=#<compiledfunction system:toplevel>,loc1=nil,loc2=0,loc3=0,loc4=nil,loc5...} [ihs=4] NIL  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1100129&group_id=4933 
From: SourceForge.net <noreply@so...>  20050110 09:41:11

Bugs item #1050304, was opened at 20041019 20:41 Message generated for change (Comment added) made by hggraebe You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1050304&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Frank Thieme (lefloyd) Assigned to: Nobody/Anonymous (nobody) Summary: primep() sometimes takes for ever Initial Comment: we had to test 3^n2 for being prime from n=1 up to 100. In MuPAD the results just need seconds, but in maxima5.9.1 (with clisp or cmucl) it takes really long for some n. For example n=37.  Comment By: HG Graebe (hggraebe) Date: 20050110 09:41 Message: Logged In: YES user_id=1193894 Frank is right, and I'm wondering about that. Here is as a first hack an implementation of the RabinMillerTest. May be it is better to have it in the lisp part, but I did not (yet) realize how and where. div(a,b):=first(divide(a,b)); modpower(a,n,m):= if n<0 then error("modpower only for n>0") else if n=0 then 1 else if n=1 then a else if oddp(n) then mod(a*modpower(mod(a*a,m),div(n,2),m),m) else mod(modpower(mod(a*a,m),div(n,2),m),m); RabinMillerTest(m):=block([q,t,j,isComposite], if m<10 then return(primep(m)), /* avoids some trivialities, a good prime test makes anyway a smallPrimesTest for a list of precompiled small primes */ q:m1, t:0, isComposite:FALSE, while evenp(q) do (q:q/2, t:t+1), /* print(q,t, m1=2^t*q), */ thru 20 do ( block([a,b], /* print("start a try"), */ a:random(m), b:modpower(a,q,m), /* print(a,b), */ if b=1 or b=1 then return(0), for j:1 thru (t1) do ( b:mod(b*b,m), /* print("loop with j=",j,", b=",b), */ if b=1 then (isComposite:TRUE,return(0)), if b=1 then j:t ), /* print("loop finished with b=",b,isComposite), */ if isComposite then return(0), if b#1 then (isComposite:TRUE,return(0)) /* , print(a,"no information in the loop") */ ), if isComposite then return(0) ), not(isComposite)); Programming first time a Maxima piece was a hard task. May be I did not yet understand it well but my opinion was that Maxima has an ugly old fashioned syntax, uncommon organization of control flow etc. Also I missed a description of the elementary arithmetic functions in the Reference Manual and the DESCRIBE part (even a remark how to include a comment), so I defined my own. And now test it: RabinMillerTest(79); CarmichaelZahlen:[1105,1729,294409, 56052361, 118901521, 172947529, 216821881, 228842209, 1299963601]; map(RabinMillerTest, CarmichaelZahlen); SmallTestExamples:makelist(3^n2,n,1,10); map(RabinMillerTest, SmallTestExamples); map(primep, SmallTestExamples); BigTestExamples:makelist(3^n2,n,40,50); map(RabinMillerTest, BigTestExamples); map(primep, BigTestExamples); Compare with MuPAD's result: BigTestExamples:=[3^n2 $ n=40..50]; map(BigTestExamples, isprime); [FALSE, TRUE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE, FALSE]  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1050304&group_id=4933 
From: SourceForge.net <noreply@so...>  20050107 16:30:53

Bugs item #1097982, was opened at 20050107 08: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=1097982&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: limit(x/(x^(log(x)))); returns wrong answer Initial Comment: In the following, the first two output are correct, but the third one is wrong. I do not know why.  (C1) limit(x/(1+x^log(x)),x,inf); (D1) 0 (C2) limit(x/(1+x^log(x)),x,inf); (D2) 0 (C3) limit(x/(x^log(x)),x,inf); (D3) INF  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1097982&group_id=4933 
From: SourceForge.net <noreply@so...>  20050107 15:21:59

Bugs item #1097915, was opened at 20050107 10:21 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=1097915&group_id=4933 Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Submitted By: Raymond Toy (rtoy) Assigned to: Nobody/Anonymous (nobody) Summary: hgfred([a,a],[1/2],z) vs hgfred([a,a],[1/2],z) Initial Comment: A&S 15.1.11 gives the result for hgfred([a,a],[1/2],z^2). maxima gets this right. But hgfred([a,a],[1/2],z^2) = hgfred([a,a],[1/2],z^2) but maxima gets the wrong answer for hgfred([a,a],[1/2],z^2).  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1097915&group_id=4933 
From: SourceForge.net <noreply@so...>  20050107 12:32:17

Bugs item #1089988, was opened at 20041222 13:41 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1089988&group_id=4933 Category: Lisp Core Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: tex does not handle bigfloats well Initial Comment: (%i1) [%pi, %e]; (%o1) [%pi,%e] (%i2) bfloat(%); (%o2) [3.14159 2653589793B0,2.718281828459045B0] (%i3) tex(%); $$\left[ 3.141592653589793B02.718281828459045B0$$ (%o3) false (%i4) x+%e; (%o4) x+%e (%i5) bfloat(%); (%o5) x+2.718281828459045B0 (%i6) tex(%); $$x2.718281828459045B0$$ (%o6) false (%i7) build_info()$ Maxima version: 5.9.1 Maxima build date: 11:36 12/13/2004 host type: i686pclinuxgnu lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.5 (%i8)  >Comment By: Barton Willis (willisbl) Date: 20050107 06:32 Message: Logged In: YES user_id=895922 I think the patch is okay. I suggest that we apply it and that we close the bug. Barton  Comment By: Robert Dodier (robert_dodier) Date: 20041231 13:58 Message: Logged In: YES user_id=501686 Patch by Barton Willis seems to work as intended  rerunning the examples above I get reasonable outputs. Maxima version: 5.9.1.1cvs Maxima build date: 0:4 12/31/2004 host type: i686redhatlinuxgnu lispimplementationtype: CLISP lispimplementationversion: 2.31 (released 20030901)  Comment By: Barton Willis (willisbl) Date: 20041226 19:32 Message: Logged In: YES user_id=895922 Please try this fix (defun texbigfloat (x l r) (append l (fpformat x) r)) Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1089988&group_id=4933 
From: SourceForge.net <noreply@so...>  20050101 17:04:11

Bugs item #1094108, was opened at 20050101 09:04 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=1094108&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Maxima crashed Dell 8400 Initial Comment: I installed Maxima on two dells for trial runs: Dell 8400 (P4, 3.2GHz) and Dell Latitude C610 (P3, 1.2 GHz). All installed tests cases run and my simple tests also work. Maxima is granted with full network access in firewall setting. And Maxima is brought up in the startup as installed. Dell 8400 has multiple user setup. When the machine is idled, sometime the machine will freeze. When the user switched, the machine always froze. I uninstalled Maxima on Dell 8400. Everything else function properly.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1094108&group_id=4933 