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

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 







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

5
(5) 
6
(1) 
7
(3) 
8
(1) 
9

10
(3) 
11
(5) 
12
(3) 
13

14

15
(5) 
16
(1) 
17
(1) 
18
(7) 
19
(1) 
20
(3) 
21

22

23

24

25

26

27
(3) 
28
(3) 
29
(1) 
30
(5) 






From: SourceForge.net <noreply@so...>  20031115 20:10:13

Bugs item #842863, was opened at 20031115 14: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=842863&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: qq.lisp is broken (numerical integration) Initial Comment: (C1) load("qq.lisp"); Load failed for C:/maxima/Maxima/share/maxima/5.9.0/share/numeric/qq .lisp  an error. Quitting. To debug this try DEBUGMODE (TRUE);) (C2) build_info(); Maxima version: 5.9.0 Maxima build date: 19:10 2/9/2003 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL25.0 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=842863&group_id=4933 
From: SourceForge.net <noreply@so...>  20031115 20:03:11

Bugs item #834417, was opened at 20031102 00:09 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=834417&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Luc Maisonobe (maisonobe) Assigned to: Nobody/Anonymous (nobody) Summary: wrong TeX output for indexed variables with exponent Initial Comment: This simple session shows the error: (C1) tex(expand((1+alpha[1])^2)); $$\alpha^2\left(1\right)+2\,\alpha_{1}+1$$ (D1) FALSE Instead of \alpha^2\left(1\right), I would have expected \alpha_{1}^2  >Comment By: Barton Willis (willisbl) Date: 20031115 14:03 Message: Logged In: YES user_id=895922 The file texmexpt.lisp has a possible fix for this bug. (C1) load("l:/texmexpt.lisp")$ (C2) tex(expand((1+alpha[1])^2)); $$\alpha_{1}^2+2\,\alpha_{1}+1$$ (D2) FALSE Here is another problem: Powers of subscripted functions don't TeX with the exponent immediately following the function (C3) f[5](x)^9; (D3) f[5](x)^9 (C4) tex(%); $$f_{5}(x)^9$$ Barton  Comment By: Barton Willis (willisbl) Date: 20031115 14:02 Message: Logged In: YES user_id=895922 The file texmexpt.lisp has a possible fix for this bug. (C1) load("l:/texmexpt.lisp")$ (C2) tex(expand((1+alpha[1])^2)); $$\alpha_{1}^2+2\,\alpha_{1}+1$$ (D2) FALSE Here is another problem: Powers of subscripted functions don't TeX with the exponent immediately following the function (C3) f[5](x)^9; (D3) f[5](x)^9 (C4) tex(%); $$f_{5}(x)^9$$ Barton  Comment By: Barton Willis (willisbl) Date: 20031115 14:01 Message: Logged In: YES user_id=895922 The file texmexpt.lisp has a possible fix for this bug. (C1) load("l:/texmexpt.lisp")$ (C2) tex(expand((1+alpha[1])^2)); $$\alpha_{1}^2+2\,\alpha_{1}+1$$ (D2) FALSE Here is another problem: Powers of subscripted functions don't TeX with the exponent immediately following the function (C3) f[5](x)^9; (D3) f[5](x)^9 (C4) tex(%); $$f_{5}(x)^9$$ Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=834417&group_id=4933 
From: SourceForge.net <noreply@so...>  20031115 20:02:26

Bugs item #834417, was opened at 20031102 00:09 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=834417&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Luc Maisonobe (maisonobe) Assigned to: Nobody/Anonymous (nobody) Summary: wrong TeX output for indexed variables with exponent Initial Comment: This simple session shows the error: (C1) tex(expand((1+alpha[1])^2)); $$\alpha^2\left(1\right)+2\,\alpha_{1}+1$$ (D1) FALSE Instead of \alpha^2\left(1\right), I would have expected \alpha_{1}^2  >Comment By: Barton Willis (willisbl) Date: 20031115 14:02 Message: Logged In: YES user_id=895922 The file texmexpt.lisp has a possible fix for this bug. (C1) load("l:/texmexpt.lisp")$ (C2) tex(expand((1+alpha[1])^2)); $$\alpha_{1}^2+2\,\alpha_{1}+1$$ (D2) FALSE Here is another problem: Powers of subscripted functions don't TeX with the exponent immediately following the function (C3) f[5](x)^9; (D3) f[5](x)^9 (C4) tex(%); $$f_{5}(x)^9$$ Barton  Comment By: Barton Willis (willisbl) Date: 20031115 14:01 Message: Logged In: YES user_id=895922 The file texmexpt.lisp has a possible fix for this bug. (C1) load("l:/texmexpt.lisp")$ (C2) tex(expand((1+alpha[1])^2)); $$\alpha_{1}^2+2\,\alpha_{1}+1$$ (D2) FALSE Here is another problem: Powers of subscripted functions don't TeX with the exponent immediately following the function (C3) f[5](x)^9; (D3) f[5](x)^9 (C4) tex(%); $$f_{5}(x)^9$$ Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=834417&group_id=4933 
From: SourceForge.net <noreply@so...>  20031115 20:01:46

Bugs item #834417, was opened at 20031102 00:09 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=834417&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Luc Maisonobe (maisonobe) Assigned to: Nobody/Anonymous (nobody) Summary: wrong TeX output for indexed variables with exponent Initial Comment: This simple session shows the error: (C1) tex(expand((1+alpha[1])^2)); $$\alpha^2\left(1\right)+2\,\alpha_{1}+1$$ (D1) FALSE Instead of \alpha^2\left(1\right), I would have expected \alpha_{1}^2  >Comment By: Barton Willis (willisbl) Date: 20031115 14:01 Message: Logged In: YES user_id=895922 The file texmexpt.lisp has a possible fix for this bug. (C1) load("l:/texmexpt.lisp")$ (C2) tex(expand((1+alpha[1])^2)); $$\alpha_{1}^2+2\,\alpha_{1}+1$$ (D2) FALSE Here is another problem: Powers of subscripted functions don't TeX with the exponent immediately following the function (C3) f[5](x)^9; (D3) f[5](x)^9 (C4) tex(%); $$f_{5}(x)^9$$ Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=834417&group_id=4933 
From: SourceForge.net <noreply@so...>  20031115 19:47:44

Bugs item #836773, was opened at 20031105 14:43 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: ntrig is broken Initial Comment: (C4) load("ntrig.mac"); (D4) ?C\:\/maxima\/Maxima\/share\/maxima\/5\.9\.0 \/share\/trigonometry\/ntrig\.mac (C5) sin(6*%pi/5); (D5) (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) (C6) float(%); (D6) 0.58778525229247 (C7) sin(float(6*%pi/5)); (D7) 0.58778525229247 Barton  >Comment By: Barton Willis (willisbl) Date: 20031115 13:47 Message: Logged In: YES user_id=895922 I wrote testing code for ntrig; both proposed fixes pass all the tests. Barton  Comment By: Stavros Macrakis (macrakis) Date: 20031110 16:07 Message: Logged In: YES user_id=588346 Please test the following replacement for ntrig. /* Some of these simplifications give results with radicals in the denominator. These could be presimplified here, but ratsimp(...),algebraic:true will also take care of it. */ eval_when([translate,batch,demo,load,loadfile], matchdeclare(n,integerp))$ tellsimpafter(sin(n*%pi/10),usin(n))$ tellsimpafter(cos(n*%pi/10),ucos(n))$ tellsimpafter(tan(n*%pi/10),usin(n)/ucos(n))$ tellsimpafter(cot(n*%pi/10),ucos(n)/usin(n))$ tellsimpafter(sec(n*%pi/10),1/ucos(n))$ tellsimpafter(csc(n*%pi/10),1/usin(n))$ usin_list: makelist( [ 0, (sqrt(5)1)/4, sqrt(2)*(sqrt(5)1)*sqrt(sqrt(5)+5)/8, (sqrt(5)+1)/4, sqrt(sqrt(5)+5)/(2*sqrt(2)), 1] [i], i,[1,2,3,4,5,6,5,4,3,2,1,2,3,4,5,6,5,4,3,2] ) /* In the pre11/2003 version, usin_list[3] was (sqrt(5)1)*sqrt(sqrt(5)+5)/(4*sqrt(2)). But the new version is better simplified, and simplifies just as well. */ $ usin(n):= (declare(n,integer), n : remainder(n,20), if n<0 then n:n+20, /* Workaround for bug in remainder */ (if n <= 10 then usin_list[n+1] /* 1origin indexing */ else usin_list[n+1]) )$ ucos(n):=usin(n+5)$  Comment By: Wolfgang Jenkner (wjenkner) Date: 20031107 11:51 Message: Logged In: YES user_id=581700 Perhaps simply rewrite USIN(N):= BLOCK([YUK:mod(N,20)], signum(YUK)*(YUK:abs(mod(YUK,10)), IF YUK=1 THEN (SQRT(5)1)/4 ELSE IF YUK=2 THEN (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) ELSE IF YUK=3 THEN (SQRT(5)+1)/4 ELSE IF YUK=4 THEN SQRT(SQRT(5)+5)/(2*SQRT(2))))$ UCOS(N):= USIN(5N)$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 
From: SourceForge.net <noreply@so...>  20031112 17:36:30

Bugs item #840848, was opened at 20031112 12:36 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=840848&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: trigexpand doesn't enter unknown functions Initial Comment: trigexpand(f(sin(x)^2)) returns f(sin(x)^2). That is, it doesn't recurse into the arguments of unknown functions. Compare with trigreduce(f(sin(2*x))). Why the discrepancy? If there is a good reason (which I doubt), this behavior should at least be documented, recommending the use of scanmap. I discovered this when trying to simplify an expression involving atan2(...sin(x)^2+cos(x)^2...). Not even atan2 is considered a 'known' function here....  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=840848&group_id=4933 
From: SourceForge.net <noreply@so...>  20031112 14:49:54

Bugs item #836780, was opened at 20031105 15:55 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(acosh(x),x,1,1) Initial Comment: The following is wrong: (C1) taylor(acosh(x),x,1,1); (D1) (%PI2*LOG(SQRT(2)+1))/2SQRT(2)*(x1)/2 Converting to log form, we do better  it's correct, I believe. But the result isn't simplified. (C2) taylor(logarc(acosh(x)),x,1,1); (D2) +SQRT(2)*SQRT(x1) (C3) ?print(%); ((MRAT SIMP (((MEXPT RATSIMP) 2 ((RAT) 1 2)) ((%LOG SIMP) ((MPLUS SIMP) $x ((MEXPT SIMP) ((MPLUS SIMP) 1 ((MEXPT SIMP) $x 2)) ((RAT) 1 2)))) ((MPLUS SIMP) 1 $x)) (#:2^(1/2)21846 #:LOG(SQRT(2)+1)21846 #:2^(1/2)21846 #:%PI21846 #:ACOSH(x)21846 #:x21846) (($x ((1 . 1)) 1 NIL #:2^(1/2)21846 . 2)) TRUNC) PS (#:2^(1/2)21846 . 2) ((1 . 1)) ((1 . 2) (#:2^(1/2)21846 1 1) . 1)) (D3) +SQRT(2)*SQRT(x1) (C4) Barton  >Comment By: Stavros Macrakis (macrakis) Date: 20031112 09:49 Message: Logged In: YES user_id=588346 The leading '+' in Taylor series is a property (bug?) of 1d display. Compare string(taylor(1,x,0,0)) => +1  Comment By: Barton Willis (willisbl) Date: 20031112 09:41 Message: Logged In: YES user_id=895922 No, I wasn't bothered by the unnecessary ratvars; rather, I was puzzled by the leading '+' in the result: +SQRT(2)*SQRT(x1) < why +sqrt(2) ...? Maybe the leading '+' is harmless; however, it's unusual. Barton  Comment By: Stavros Macrakis (macrakis) Date: 20031110 15:19 Message: Logged In: YES user_id=588346 This was already reported in bug report # 623165. I am not sure what you mean by "the result isn't simplified". Are you thinking that the unnecessary ratvars are a problem? That is the usual way that CREs work: in a fresh Maxima, try (?print(rat(x)), ?print(rat(y)), ?print(rat(x)) )$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 
From: SourceForge.net <noreply@so...>  20031112 14:41:54

Bugs item #836780, was opened at 20031105 14:55 Message generated for change (Comment added) made by willisbl You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(acosh(x),x,1,1) Initial Comment: The following is wrong: (C1) taylor(acosh(x),x,1,1); (D1) (%PI2*LOG(SQRT(2)+1))/2SQRT(2)*(x1)/2 Converting to log form, we do better  it's correct, I believe. But the result isn't simplified. (C2) taylor(logarc(acosh(x)),x,1,1); (D2) +SQRT(2)*SQRT(x1) (C3) ?print(%); ((MRAT SIMP (((MEXPT RATSIMP) 2 ((RAT) 1 2)) ((%LOG SIMP) ((MPLUS SIMP) $x ((MEXPT SIMP) ((MPLUS SIMP) 1 ((MEXPT SIMP) $x 2)) ((RAT) 1 2)))) ((MPLUS SIMP) 1 $x)) (#:2^(1/2)21846 #:LOG(SQRT(2)+1)21846 #:2^(1/2)21846 #:%PI21846 #:ACOSH(x)21846 #:x21846) (($x ((1 . 1)) 1 NIL #:2^(1/2)21846 . 2)) TRUNC) PS (#:2^(1/2)21846 . 2) ((1 . 1)) ((1 . 2) (#:2^(1/2)21846 1 1) . 1)) (D3) +SQRT(2)*SQRT(x1) (C4) Barton  >Comment By: Barton Willis (willisbl) Date: 20031112 08:41 Message: Logged In: YES user_id=895922 No, I wasn't bothered by the unnecessary ratvars; rather, I was puzzled by the leading '+' in the result: +SQRT(2)*SQRT(x1) < why +sqrt(2) ...? Maybe the leading '+' is harmless; however, it's unusual. Barton  Comment By: Stavros Macrakis (macrakis) Date: 20031110 14:19 Message: Logged In: YES user_id=588346 This was already reported in bug report # 623165. I am not sure what you mean by "the result isn't simplified". Are you thinking that the unnecessary ratvars are a problem? That is the usual way that CREs work: in a fresh Maxima, try (?print(rat(x)), ?print(rat(y)), ?print(rat(x)) )$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 
From: SourceForge.net <noreply@so...>  20031111 23:28:11

Bugs item #840360, was opened at 20031111 18:28 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=840360&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: qunit(4) internal errors Initial Comment: qunit of a square, e.g. qunit(4), gives the internal error "Floatingpoint exception." It should either give a friendly Maxima error, or possibly return 1  isn't the "quadratic" number field of sqrt(n^2) simply the rationals? qunit(2) goes into an infinite loop. It should give an error, since we're dealing with *real* quadratic number fields, but wouldn't it be nice if...  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=840360&group_id=4933 
From: SourceForge.net <noreply@so...>  20031111 15:26:05

Bugs item #814957, was opened at 20030930 08:17 Message generated for change (Comment added) made by pukruppa You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=814957&group_id=4933 Category: Xmaxima Group: Fix for 5.9.0 Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: No internet => No xMaxima Initial Comment: The windows compiled version of xMaxima refuses to work if a connection to internet is not available and working. A message saying "Error starting Maxima: Could not open a socket" is issued. After that the window xMaxima keeps mute, no prompt, no echo, nothing... If we have no Internet connection how can we keep working with Maxima through the xMaxima interface?  Comment By: Peter Ulrich Kruppa (pukruppa) Date: 20031111 15:26 Message: Logged In: YES user_id=778327 Some of my students have reported this problem to me, too, although everything on my laptop (WinXP Home Edition) works fine  even without internet connection. The only significant difference to their machines, I can think of, is that I once activated XP's compatibility mode for older programs (sorry, I don't know the exact words since I am running XP in german  do search for "compatibility" in the help menu and reinstall maxima with it). Just one possibility. Regards, Uli.  Comment By: Stavros Macrakis (macrakis) Date: 20030930 18:25 Message: Logged In: YES user_id=588346 An active Internet connection is not required by xMaxima. However, you must have sockets installed, working, and enabled, because Maxima uses them for interprocess communication. Firewall software may be closing off sockets, for example.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=814957&group_id=4933 
From: SourceForge.net <noreply@so...>  20031111 15:07:58

Bugs item #839926, was opened at 20031111 13:09 Message generated for change (Comment added) made by pukruppa You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=839926&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: save(...,all) can't be reloaded Initial Comment: (discovered by Peter Ulrich Kruppa) save("file",all)$ loadfile("file")$ => Load failed The problem is that the symbol DVA (used in NUSUM) has a property which is a compiled closure (a translated mmacro), which is not printable. The workaround is to kill DVA before saving. The problem is some interaction between Translate and Save which I don't have the time to track down.  Comment By: Peter Ulrich Kruppa (pukruppa) Date: 20031111 15:07 Message: Logged In: YES user_id=778327 Since I didn't have the time to find out what dva is and what it does, here another simplicistic workaround: open the session file with the editor of your choice and delete all the stuff that comes _after_ your last input, probably something like: (DSKSETQ $C6 '(($SAVE) &session02 $ALL)) (ADDLABEL '$C6). loadfile(...) will work correctly then. Regards, Peter Ulrich Kruppa  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=839926&group_id=4933 
From: SourceForge.net <noreply@so...>  20031111 13:09:43

Bugs item #839926, was opened at 20031111 08:09 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=839926&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: save(...,all) can't be reloaded Initial Comment: (discovered by Peter Ulrich Kruppa) save("file",all)$ loadfile("file")$ => Load failed The problem is that the symbol DVA (used in NUSUM) has a property which is a compiled closure (a translated mmacro), which is not printable. The workaround is to kill DVA before saving. The problem is some interaction between Translate and Save which I don't have the time to track down.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=839926&group_id=4933 
From: SourceForge.net <noreply@so...>  20031111 04:07:22

Bugs item #839732, was opened at 20031110 21: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=839732&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Robert Dodier (robert_dodier) Assigned to: Nobody/Anonymous (nobody) Summary: plot2d fails if function returns a rational number Initial Comment: I am running Maxima 5.9.0 (clisp 2.29) on Redhat linux 7.1. plot2d('(1/2),[x,0,1]); fails with: ***  argument to COMMONLISP:FLOAT should be a real number: ((RAT SIMP) 1 2) However, plot2d('(float(1/2)),[x,0,1]); and plot2d('(0.5),[x,0,1]); both succeed and plot a horizontal line segment, as expected. Other examples such as 1/4, 2/3, 3/5, all fail. However, 1, 3, 5 (whole numbers) succeed. The closest alreadyreported bug seems to be 572668, but the fix for that bug made it into version 5.9.0 (as version 1.17 of plot.lisp) so that didn't fix this new bug. I have also tried the same thing with GCL (sorry, forgot the version number) and there is again an error although the message is different.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=839732&group_id=4933 
From: SourceForge.net <noreply@so...>  20031110 22:07:03

Bugs item #836773, was opened at 20031105 15:43 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: ntrig is broken Initial Comment: (C4) load("ntrig.mac"); (D4) ?C\:\/maxima\/Maxima\/share\/maxima\/5\.9\.0 \/share\/trigonometry\/ntrig\.mac (C5) sin(6*%pi/5); (D5) (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) (C6) float(%); (D6) 0.58778525229247 (C7) sin(float(6*%pi/5)); (D7) 0.58778525229247 Barton  >Comment By: Stavros Macrakis (macrakis) Date: 20031110 17:07 Message: Logged In: YES user_id=588346 Please test the following replacement for ntrig. /* Some of these simplifications give results with radicals in the denominator. These could be presimplified here, but ratsimp(...),algebraic:true will also take care of it. */ eval_when([translate,batch,demo,load,loadfile], matchdeclare(n,integerp))$ tellsimpafter(sin(n*%pi/10),usin(n))$ tellsimpafter(cos(n*%pi/10),ucos(n))$ tellsimpafter(tan(n*%pi/10),usin(n)/ucos(n))$ tellsimpafter(cot(n*%pi/10),ucos(n)/usin(n))$ tellsimpafter(sec(n*%pi/10),1/ucos(n))$ tellsimpafter(csc(n*%pi/10),1/usin(n))$ usin_list: makelist( [ 0, (sqrt(5)1)/4, sqrt(2)*(sqrt(5)1)*sqrt(sqrt(5)+5)/8, (sqrt(5)+1)/4, sqrt(sqrt(5)+5)/(2*sqrt(2)), 1] [i], i,[1,2,3,4,5,6,5,4,3,2,1,2,3,4,5,6,5,4,3,2] ) /* In the pre11/2003 version, usin_list[3] was (sqrt(5)1)*sqrt(sqrt(5)+5)/(4*sqrt(2)). But the new version is better simplified, and simplifies just as well. */ $ usin(n):= (declare(n,integer), n : remainder(n,20), if n<0 then n:n+20, /* Workaround for bug in remainder */ (if n <= 10 then usin_list[n+1] /* 1origin indexing */ else usin_list[n+1]) )$ ucos(n):=usin(n+5)$  Comment By: Wolfgang Jenkner (wjenkner) Date: 20031107 12:51 Message: Logged In: YES user_id=581700 Perhaps simply rewrite USIN(N):= BLOCK([YUK:mod(N,20)], signum(YUK)*(YUK:abs(mod(YUK,10)), IF YUK=1 THEN (SQRT(5)1)/4 ELSE IF YUK=2 THEN (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) ELSE IF YUK=3 THEN (SQRT(5)+1)/4 ELSE IF YUK=4 THEN SQRT(SQRT(5)+5)/(2*SQRT(2))))$ UCOS(N):= USIN(5N)$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 
From: SourceForge.net <noreply@so...>  20031110 20:19:39

Bugs item #836780, was opened at 20031105 15:55 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(acosh(x),x,1,1) Initial Comment: The following is wrong: (C1) taylor(acosh(x),x,1,1); (D1) (%PI2*LOG(SQRT(2)+1))/2SQRT(2)*(x1)/2 Converting to log form, we do better  it's correct, I believe. But the result isn't simplified. (C2) taylor(logarc(acosh(x)),x,1,1); (D2) +SQRT(2)*SQRT(x1) (C3) ?print(%); ((MRAT SIMP (((MEXPT RATSIMP) 2 ((RAT) 1 2)) ((%LOG SIMP) ((MPLUS SIMP) $x ((MEXPT SIMP) ((MPLUS SIMP) 1 ((MEXPT SIMP) $x 2)) ((RAT) 1 2)))) ((MPLUS SIMP) 1 $x)) (#:2^(1/2)21846 #:LOG(SQRT(2)+1)21846 #:2^(1/2)21846 #:%PI21846 #:ACOSH(x)21846 #:x21846) (($x ((1 . 1)) 1 NIL #:2^(1/2)21846 . 2)) TRUNC) PS (#:2^(1/2)21846 . 2) ((1 . 1)) ((1 . 2) (#:2^(1/2)21846 1 1) . 1)) (D3) +SQRT(2)*SQRT(x1) (C4) Barton  >Comment By: Stavros Macrakis (macrakis) Date: 20031110 15:19 Message: Logged In: YES user_id=588346 This was already reported in bug report # 623165. I am not sure what you mean by "the result isn't simplified". Are you thinking that the unnecessary ratvars are a problem? That is the usual way that CREs work: in a fresh Maxima, try (?print(rat(x)), ?print(rat(y)), ?print(rat(x)) )$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 
From: SourceForge.net <noreply@so...>  20031110 02:13:17

Bugs item #801244, was opened at 20030905 12:13 Message generated for change (Comment added) made by mcelrath You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=801244&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: MAX(a+b,c) is NOT equal to MAX(c,a+b) Initial Comment:  Maxima version: 5.9.0 Maxima build date: 19:1 8/5/2003 host type: i586pclinuxgnu lispimplementationtype: CLISP lispimplementationversion: 2.27 (released 20010717) (built 3223390905) (memory 3269088151)  (C1) assume(a+b>c)$ (C2) MAX(a+b,c); MAX(c,a+b); (D2) b + a (C3) (D3) MAX(c, b + a) #why it is not simplified ? (C4) is(MAX(a+b,c)=MAX(c,a+b)); (D4) FALSE #why FALSE if must be TRUE ? P.S. Appreciate if you exclude javascripts from the site. Alexander VIDYBIDA vidybida@...  Comment By: Bob McElrath (mcelrath) Date: 20031109 20:13 Message: Logged In: YES user_id=10643 It seems this is due to: (C1) assume(a+b>c)$ (C2) is(a+b>c); (D2) TRUE (C3) is(a+b<=c); MACSYMA was unable to evaluate the predicate: b + a <= C  an error. Quitting. To debug this try DEBUGMODE(TRUE);) Maxima should be able to figure out that a>b implies !(b<=a), but it doesn't... looks like the relevant file is compar.lisp.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=801244&group_id=4933 
From: SourceForge.net <noreply@so...>  20031108 02:17:43

Bugs item #838301, was opened at 20031107 21:17 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=838301&group_id=4933 Category: Share Libraries Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stephen Leake (stephen_leake) Assigned to: Nobody/Anonymous (nobody) Summary: vect negate cross product simplification Initial Comment: load (vect); declare (n, nonscalar)$ declare (r, nonscalar)$ declare (a, scalar)$ a*r~n; The current result is  a ( n) ~ r In vect.mac, there is this declaration: TELLSIMP(ETRUE~LESSP, LESSP~ETRUE) $ If I change this to: TELLSIMP(ETRUE~LESSP, (LESSP~ETRUE)) $ the above expression reduces to: a n ~ r which is better.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=838301&group_id=4933 
From: SourceForge.net <noreply@so...>  20031107 17:51:03

Bugs item #836773, was opened at 20031105 21:43 Message generated for change (Comment added) made by wjenkner You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: ntrig is broken Initial Comment: (C4) load("ntrig.mac"); (D4) ?C\:\/maxima\/Maxima\/share\/maxima\/5\.9\.0 \/share\/trigonometry\/ntrig\.mac (C5) sin(6*%pi/5); (D5) (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) (C6) float(%); (D6) 0.58778525229247 (C7) sin(float(6*%pi/5)); (D7) 0.58778525229247 Barton  >Comment By: Wolfgang Jenkner (wjenkner) Date: 20031107 18:51 Message: Logged In: YES user_id=581700 Perhaps simply rewrite USIN(N):= BLOCK([YUK:mod(N,20)], signum(YUK)*(YUK:abs(mod(YUK,10)), IF YUK=1 THEN (SQRT(5)1)/4 ELSE IF YUK=2 THEN (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) ELSE IF YUK=3 THEN (SQRT(5)+1)/4 ELSE IF YUK=4 THEN SQRT(SQRT(5)+5)/(2*SQRT(2))))$ UCOS(N):= USIN(5N)$  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 
From: SourceForge.net <noreply@so...>  20031107 17:44:32

Bugs item #706562, was opened at 20030319 23:07 Message generated for change (Comment added) made by wjenkner You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=706562&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: mod(2,4) => 2 not 2 Initial Comment: mod(2,4) => 2 It should be 2 to be consistent with mod(2,4), mod(6,4), etc. Otherwise there are FIVE distinct values for mod (x,4) instead of four. Maxima 5.9.0 GCL 2.5.0 mingw Windows 2000  >Comment By: Wolfgang Jenkner (wjenkner) Date: 20031107 18:44 Message: Logged In: YES user_id=581700 GCL imports some symbols for modular arithmetic from the SYSTEM package (see maximapackage.lisp). In particular, >(lispimplementationversion) "GCL25.2000000000000002" >(let ((si:modulus 4)) (si:cmod 2)) 2 >(let ((si:modulus 4)) (si:cmod 2)) 2 The other lisps define this in rat3a.lisp * (let ((modulus 4)) (cmod 2)) 2 * (let ((modulus 4)) (cmod 2)) 2  Comment By: Wolfgang Jenkner (wjenkner) Date: 20030322 16:42 Message: Logged In: YES user_id=581700 For both Clisp and SBCL I get (C1) makelist(mod(i,4),i,4,1); (D1) [0, 1, 2,  1] Wolfgang  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=706562&group_id=4933 
From: SourceForge.net <noreply@so...>  20031107 17:38:51

Bugs item #836708, was opened at 20031105 20:00 Message generated for change (Comment added) made by wjenkner You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836708&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: rat/tellrat/modulus:2 can return 1 Initial Comment: algebraic:true$ modulus:2$ tellrat(a^2+a)$ rat(a^2) =>  a With modulus=2, 1 is supposed to simplify to 1. This can be fixed (bizarrely) by ev: ev(%) => a  >Comment By: Wolfgang Jenkner (wjenkner) Date: 20031107 18:38 Message: Logged In: YES user_id=581700 On both SBCL and CLISP I get the expected result (C4) rat(a^2); (D4)/R/ a I'd guess this is the same GCL specific misfeature as #706562 (mod(2,4) => 2 not 2).  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836708&group_id=4933 
From: SourceForge.net <noreply@so...>  20031106 21:49:22

Bugs item #836704, was opened at 20031105 19:56 Message generated for change (Comment added) made by wjenkner You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836704&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: gendiff is all bugs: should be deprecated Initial Comment: The gendiff package only gets the most trivial cases right, and should probably be deprecated until it is fixed. load("gendif") gendiff(x,x,n) => 'diff(x,x,n) which is correct, except that gendiff is supposed to give an explicit form *without* 'diff's. After all, it does it for x^2: genfact(2,n,1)*x^(2n) so why doesn't it give genfact(1,n,1)*x^(1n) in this case? Sure, that's clumsy, but it's explicit.  gendiff((x1)^2,x,n) => 0 (!!) The ratexpand in gendiffpow is presumably supposed to be ratexpand(x1^x2); but though that fix works in this particular case, it doesn't work for sin(x)^2.  q: gendiff(x*(x+1),x,n) => 'DIFF(x,x,n)*'SUM(BINOMIAL(n,I)*'DIFF(x,x,I),I,0,n) which looks very impressive, but is wrong. In particular, subst(0,n,q) => x^2 (!!) and subst(1,n,q) => (x+1) (!!). Correct answers are of course x*(x+1) and 2*x+1.  gendiff(x^x,x,n) => x^(xn)*GENFACT(x,n,1) (!!) but diff(x^x,x,1) => x^x*(log(x)+1) (OK) gendiff is assuming that all powers are free of the differentiation variable. To fix, add a freeof clause in gendiffpow: ...if x1=x <<<and freeof(x,x2)>>> then ...  gendiff is also pretty limited, not even handling sin(x) ===> sin(x+n*%pi/2) sinh(x) ===> (%e^x+(1)^n*%e^x)/2 (also expressible as pure hyperbolics)  >Comment By: Wolfgang Jenkner (wjenkner) Date: 20031106 22:49 Message: Logged In: YES user_id=581700 > q: gendiff(x*(x+1),x,n) => > 'DIFF(x,x,n)*'SUM(BINOMIAL(n,I)*'DIFF(x,x,I),I,0,n) > > which looks very impressive, but is wrong. Just a case sensitivity bug. Here's a simpler example: (C2) gendiff(x+1,x,k); n d x (D2)  n dx $N is bound to $k but $n, which is passed to the recursive calls to $GENDIFF, is not. I propose to just upcase the whole file. Downcasing would also be possible, but in this case we would probably want to keep the name GENDIFF itself uppercase, such that gendiff(...) and GENDIFF(...) both work. I prefer upcasing also because this way people may spot more easily what changed (most of the file content being uppercase).  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836704&group_id=4933 
From: SourceForge.net <noreply@so...>  20031105 20:55:56

Bugs item #836780, was opened at 20031105 14:55 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=836780&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(acosh(x),x,1,1) Initial Comment: The following is wrong: (C1) taylor(acosh(x),x,1,1); (D1) (%PI2*LOG(SQRT(2)+1))/2SQRT(2)*(x1)/2 Converting to log form, we do better  it's correct, I believe. But the result isn't simplified. (C2) taylor(logarc(acosh(x)),x,1,1); (D2) +SQRT(2)*SQRT(x1) (C3) ?print(%); ((MRAT SIMP (((MEXPT RATSIMP) 2 ((RAT) 1 2)) ((%LOG SIMP) ((MPLUS SIMP) $x ((MEXPT SIMP) ((MPLUS SIMP) 1 ((MEXPT SIMP) $x 2)) ((RAT) 1 2)))) ((MPLUS SIMP) 1 $x)) (#:2^(1/2)21846 #:LOG(SQRT(2)+1)21846 #:2^(1/2)21846 #:%PI21846 #:ACOSH(x)21846 #:x21846) (($x ((1 . 1)) 1 NIL #:2^(1/2)21846 . 2)) TRUNC) PS (#:2^(1/2)21846 . 2) ((1 . 1)) ((1 . 2) (#:2^(1/2)21846 1 1) . 1)) (D3) +SQRT(2)*SQRT(x1) (C4) Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&group_id=4933 
From: SourceForge.net <noreply@so...>  20031105 20:43:24

Bugs item #836773, was opened at 20031105 14:43 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=836773&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: ntrig is broken Initial Comment: (C4) load("ntrig.mac"); (D4) ?C\:\/maxima\/Maxima\/share\/maxima\/5\.9\.0 \/share\/trigonometry\/ntrig\.mac (C5) sin(6*%pi/5); (D5) (SQRT(5)1)*SQRT(SQRT(5)+5)/(4*SQRT(2)) (C6) float(%); (D6) 0.58778525229247 (C7) sin(float(6*%pi/5)); (D7) 0.58778525229247 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836773&group_id=4933 
From: SourceForge.net <noreply@so...>  20031105 19:01:15

Bugs item #836708, was opened at 20031105 14:00 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=836708&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: rat/tellrat/modulus:2 can return 1 Initial Comment: algebraic:true$ modulus:2$ tellrat(a^2+a)$ rat(a^2) =>  a With modulus=2, 1 is supposed to simplify to 1. This can be fixed (bizarrely) by ev: ev(%) => a  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836708&group_id=4933 
From: SourceForge.net <noreply@so...>  20031105 18:57:21

Bugs item #836704, was opened at 20031105 13:56 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=836704&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: gendiff is all bugs: should be deprecated Initial Comment: The gendiff package only gets the most trivial cases right, and should probably be deprecated until it is fixed. load("gendif") gendiff(x,x,n) => 'diff(x,x,n) which is correct, except that gendiff is supposed to give an explicit form *without* 'diff's. After all, it does it for x^2: genfact(2,n,1)*x^(2n) so why doesn't it give genfact(1,n,1)*x^(1n) in this case? Sure, that's clumsy, but it's explicit.  gendiff((x1)^2,x,n) => 0 (!!) The ratexpand in gendiffpow is presumably supposed to be ratexpand(x1^x2); but though that fix works in this particular case, it doesn't work for sin(x)^2.  q: gendiff(x*(x+1),x,n) => 'DIFF(x,x,n)*'SUM(BINOMIAL(n,I)*'DIFF(x,x,I),I,0,n) which looks very impressive, but is wrong. In particular, subst(0,n,q) => x^2 (!!) and subst(1,n,q) => (x+1) (!!). Correct answers are of course x*(x+1) and 2*x+1.  gendiff(x^x,x,n) => x^(xn)*GENFACT(x,n,1) (!!) but diff(x^x,x,1) => x^x*(log(x)+1) (OK) gendiff is assuming that all powers are free of the differentiation variable. To fix, add a freeof clause in gendiffpow: ...if x1=x <<<and freeof(x,x2)>>> then ...  gendiff is also pretty limited, not even handling sin(x) ===> sin(x+n*%pi/2) sinh(x) ===> (%e^x+(1)^n*%e^x)/2 (also expressible as pure hyperbolics)  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836704&group_id=4933 