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

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 







1
(1) 
2

3
(4) 
4
(1) 
5
(5) 
6
(4) 
7
(2) 
8

9
(100) 
10
(57) 
11
(1) 
12
(11) 
13
(2) 
14
(2) 
15
(2) 
16

17
(5) 
18
(1) 
19

20

21

22

23

24
(1) 
25

26
(2) 
27
(4) 
28
(1) 
29

30
(5) 






From: SourceForge.net <noreply@so...>  20060409 21:45:57

Bugs item #1467368, was opened at 20060409 17: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=1467368&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: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: logcontract returns unsimplified expr Initial Comment: expr: log(%e*k)log(%e^1*k)$ logcontract(expr) => log(%e^2) which is unsimplified: expand(logcontract(expr)) => 2 Maxima version: 5.9.3 Maxima build date: 0:52 3/20/2006 host type: i686pcmingw32 lispimplementationtype: GNU Common Lisp (GCL) lispimplementationversion: GCL 2.6.7  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1467368&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 21:06:52

Bugs item #708947, was opened at 20030324 13:46 Message generated for change (Comment added) made by macrakis You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=708947&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  Taylor Group: None >Status: Closed >Resolution: Fixed Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor of taylor => fatal error Initial Comment: taylor(taylor(a*b,a,0,1),b,1/a,1) => fatal error Maxima 5.9.0 GCL 2.5.0 Windows 2000 mingw  >Comment By: Stavros Macrakis (macrakis) Date: 20060409 17:06 Message: Logged In: YES user_id=588346 OK on Maxima version: 5.9.3 Maxima build date: 0:52 3/20/2006 host type: i686pcmingw32 lispimplementationtype: GNU Common Lisp (GCL) lispimplementationversion: GCL 2.6.7 Closing bug  Comment By: Robert Dodier (robert_dodier) Date: 20050410 23:13 Message: Logged In: YES user_id=501686 I've tried this with Maxima cvs rev 2005/04/09 on gcl 2.6.6, fedora linux. I get the "attempt to expand" for all three examples given here  taylor(taylor(a*b,a,0,1),b,1/a,1) => "attempt to expand ..." taylor(taylor(a*b,a,0,1),b,a,1) => "attempt to expand ..." taylor(a*b,a,0,1,b,a,1) => "attempt to expand ..." Looking at the log for src/hayat.lisp, there have been some code changes over the years (sloop replacement and cleanups contributed by Andreas Eder). Maybe this bug can be retested on Windows.  Comment By: Stavros Macrakis (macrakis) Date: 20030324 13:57 Message: Logged In: YES user_id=588346 Simpler example: taylor(taylor(a*b,a,0,1),b,a,1) Note that taylor(a*b,a,0,1,b,a,1) sensibly gives a user oriented error: TAYLOR: attempt to expand a b at a point depending on b  an error.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=708947&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:46:29

Bugs item #649428, was opened at 20021206 02:42 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=649428&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Martin Rubey (kratt5) Assigned to: Nobody/Anonymous (nobody) Summary: new simpsum / taylor bug Initial Comment: This is a new simpsum / taylor bug. Note that I used a patched version of Maxima, as described in the simpsum bug thread: combin.lisp 933,935c927,929 < ;; nil ;; Kratt5 26.11.2002 < ; #cl < (let (*a *n (var *var*)) ; freevar expects "var", not "*var*"  > nil > #cl > (let (*a *n) (C1) build_info(); Maxima version: 5.9.0rc3 Maxima build date: 21:0 11/26/2002 host type: i686pclinuxgnu lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL25.0 (D1) (C2) factsum3(mt,ej):=sum((1)^(k+1)/(k*mt^k)*sum((1l)^kl^k,l,1,ej),k,1,INF)$ (C3) taylor(factsum3(mt,ej),[mt,0,3,ASYMP]); ej ej ej ==== ==== ==== \ \ \ 3 2 > ( 2 l + 1) > ( 2 l + 1) > ( 2 l + 3 l  3 l + 1) / / / ==== ==== ==== l = 1 l = 1 l = 1 (D3)/T/    +  mt 2 3 2 mt 3 mt + . . . (C4) taylor(factsum3(mt,ej),[mt,0,3,ASYMP]),simpsum; 2 2 4 2 l l l + l (D4)/T/   +    + . . . mt 2 3 2 mt 6 mt (C5) This is nearly correct, only l should be changed to ej... Martin  Comment By: Robert Dodier (robert_dodier) Date: 20051125 09:02 Message: Logged In: YES user_id=501686 Reopening this bug report; it was closed by mistake. The reported problem is still present in cvs Maxima. For the record here is a more readable presentation of the same example. (%i2) factsum3(mt,ej):=sum((1)^(k+1)/(k*mt^k)*sum((1l)^kl^k,l,1,ej),k,1,inf); (%o2) factsum3(mt,ej):=sum( (1)^(k+1)/(k*mt^k)*sum((1l)^kl^k,l,1,ej),k,1, inf) (%i3) taylor(factsum3(mt,ej),[mt,0,3,asymp]); (%o3) ('sum(2*l+1,l,1,ej))/mt('sum(2*l+1,l,1,ej))/(2*mt^2) +('sum(2*l^3+3*l^23*l+1,l,1,ej)) /(3*mt^3) (%i4) ''%, simpsum; (%o4) (6*ej^2*mt^23*ej^2*mt+ej^4+ej^2)/(6*mt^3) < OK (%i5) taylor(factsum3(mt,ej),[mt,0,3,asymp]), simpsum; (%o5) l^2/mt+l^2/(2*mt^2)(l^4+l^2)/(6*mt^3) < OOPS  Comment By: Robert Dodier (robert_dodier) Date: 20051124 01:43 Message: Logged In: YES user_id=501686 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=649428&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:49

Bugs item #1042702, was opened at 20041007 19:29 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042702&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor TAN(a*ATAN(q/a)) stack overflow Initial Comment: taylor( TAN(a*ATAN(q/a)), a, 1, 1) > stack overflow No stack overflow with q=1  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1042702&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:49

Bugs item #1203443, was opened at 20050517 04:15 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1203443&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: taylor about minf Initial Comment: Taylor expansions about minf (minus infinity) are broken; for example (%i1) taylor(1/x^2,x,minf,2); (%o1) +((1/x)^2) Three problems: %o1 is wrong, %o1 isn't simplified, and with display2d : false, Taylor expansions do not display as /T/ <stuff> + ... Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1203443&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:49

Bugs item #1220979, was opened at 20050615 01:18 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1220979&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: Taylor Initial Comment: M : MATRIX([a,b],[b,c]); L : eigenvalues(M) Tc : taylor(L[1][1],c,0,2); Ta : taylor(L[1][1],a,0,2); But "a" and "c" is simmetric. it work for Tc, but don't work for Ta.  Comment By: Barton Willis (willisbl) Date: 20050703 04:28 Message: Logged In: YES user_id=895922 This is a bug connected with the default algorithm for gcd (subres). To work around the bug, set 'gcd' to 'spmod' (%i1) M : MATRIX([a,b],[b,c])$ (%i2) L : eigenvalues(M)$ (%i3) Ta : taylor(L[1][1],a,0,2); Quotient by a polynomial of higher degree  an error. Quitting. To debug this try DEBUGMODE(TRUE); (%i4) gcd : spmod; (%o4) SPMOD (%i5) Ta : taylor(L[1][1],a,0,2); (%o5) (csqrt(c^2+4*b^2))/2+((c^2+sqrt(c^2+4*b^2)*c+4*b^2)*a)/(2*c^2+8*b^2)(sqrt(c^2+4*b^2)*b^2*a^2)/(c^4+8*b^2*c^2+16*b^4)+... Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1220979&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:48

Bugs item #658850, was opened at 20021226 17:47 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=658850&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  Taylor Group: None Status: Open Resolution: Accepted Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: Taylor "invalid call to varexpand" Initial Comment: Maxima 5.5 on Windows GCL C15) taylor(exp(x)/x,x,0,5); Invalid call to varexpand  an error. This is an internal error. Though this expression does not have a Taylor/Laurent expansion, it should still give a useroriented error, not an internal error.  Comment By: Stavros Macrakis (macrakis) Date: 20030602 14:00 Message: Logged In: YES user_id=588346 expr: (nn^(nn+1)(nn2)^nn*nn+2*(nn2)^nn)*p/(2*nn^nn*ppy) taylor(expr,nn,inf,0) gives Invalid call to varexpand but taylor(expand(expr),nn,inf,0) returns a result which looks correct.  Comment By: Martin Rubey (kratt5) Date: 20030123 03:39 Message: Logged In: YES user_id=651552 It should read taylor(exp(1/x)/x,x,0,5); Invalid call to varexpand  an error. Quitting. To debug this try DEBUGMODE(TRUE);  Comment By: Stavros Macrakis (macrakis) Date: 20030122 06:48 Message: Logged In: YES user_id=588346 Sorry, I must have miscopied the error case (or perhaps some flags were set?), and I can't reconstruct it. The example I gave does work in 5.5, and does have a Laurent expansion. I am deleting this bug report and marking it as "invalid".  Comment By: Martin Rubey (kratt5) Date: 20030122 06:39 Message: Logged In: YES user_id=651552 I do not quite understand. On my (slightly patched) maxima I get:  Maxima version: 5.9.0rc3 Maxima build date: 21:0 11/26/2002 host type: i686pclinuxgnu lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL25.0  The above information is also available from the Maxima function build_info(). (D1) (C2) taylor(exp(x)/x,x,0,5); 2 3 4 5 1 x x x x x (D2)/T/   1 +    +    +  + . . . x 2 6 24 120 720 (C3) which seems to be correct... Martin  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=658850&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:48

Bugs item #803247, was opened at 20030909 09:51 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=803247&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Camm Maguire (yycamm) Assigned to: Nobody/Anonymous (nobody) Summary: Infinite loop in Taylor(asin(sqrt(1x^2)/(1+eps))...)/workar Initial Comment: forwarded 198466 maxima@... tags 198466 +upstream tags 198466 +confirmed thanks A Debian user submitted this report, which I'm also registering at the website. Any ideas on a fix? Take care, ============================================================================= When trying to do some complicated Taylor series expansions, maxima hangs; I believe it is in an infinite loop, as xload shows the cpu saturated. An example is: f(x):= asin(sqrt(1x^2)/(1+eps)); taylor(f(x),x,0,2); which produces this condition. While trying to isolate the bug, I tried: g(x):=asin(sqrt(1x^2)); taylor(g(x),x,0,2); which produces the error message: "Taylor encountered an unfamiliar singularity in ABS(X)" I am unable to find any documentation on this "unfamiliar singularity" error message. However, from looking at the second derivative of G(x), I think I can see what's causing the problem: the sqrt() in the above expression is symmetrical about x=0, but has a cusp there. However, this ought to cause taylor() to produce an error message in the first example, instead of looping; this is not very userfriendly. What's the difference between a familiar singularity and an unfamiliar one, anyway? =============================================================================  Camm Maguire camm@... ========================================================================== "The earth is but one country, and mankind its citizens."  Baha'u'llah  Comment By: Stavros Macrakis (macrakis) Date: 20030918 10:33 Message: Logged In: YES user_id=588346 In the first problem, the bug appears to be in the default gcd; you can get around this by changing gcd routine: gcd:'spmod$ The second case is a bug in Taylor. To get around it, first expand the inner expression, then plug it into asin: inner: taylor(sqrt(1x^2)+err*x^6,x,0,6); taylor(asin(inner),x,0,5); In general, you have to be careful about the orders of expansion  you may need more terms in the inner expression than the whole expression. That is the sort of thing Taylor is supposed to do for you, but obviously there is a problem in this case. That is why I have added an explicit error term  to be able to trace it through. The second derivative business is, I believe, a red herring  the absolute value has to do with returning the principal value of sqrt, but the Taylor expansion doesn't care about principal values, it cares about analytic continuation.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=803247&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:48

Bugs item #788933, was opened at 20030814 14:35 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=788933&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor unhandled multivar datum comparison Initial Comment: taylor(x^y,x,0,1,y,0,1) and taylor(taylor(x^y,y,0,1),x,0,1) => Break: Unhandled multivar datum comparison though taylor(x^y,y,0,1,x,0,1); => 1+(LOG(x)+...)*y+... and TAYLOR(TAYLOR(x^y,x,0,1),y,0,1) => 1+(LOG(x)+...)*y+... and taylor(x^y,[x,y],[0,0],1); = taylor(x^y,[y,x],[0,0],1); => 1+LOG(x)*y+...  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=788933&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #932302, was opened at 20040409 07:45 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=932302&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: partfrac/ratnumer of taylor bad: taychk2rat/FIX Initial Comment: expr: 1/(x^21)$ texpr: taylor(expr,x,1,1)$ ptexpr: partfrac(texpr,x) => (x+4/(x1)3)/8 NO! This is algebraically correct, but not in partfrac form. The correct answer is given by: partfrac(ratdisrep(texpr),x) == partfrac(ptexpr,x) => 1/(2*(x1))+(x3)/8 The immediate fix is to replace (DESETQ (RATFORM . EXP) (TAYCHK2RAT EXP)) with (DESETQ (RATFORM . EXP) (RATF (TAYCHK2RAT EXP))) however, I wonder if TAYCHK2RAT shouldn't be doing this. Compare: ratnumer(taylor(x+1/x,x,0,1)) => x+1/x ???  Comment By: Stavros Macrakis (macrakis) Date: 20040409 08:05 Message: Logged In: YES user_id=588346 taychk2rat is used in four places in Maxima: ratnumer ratdenom partfrac horner In all of them, it would be better if taychk2rat performed the ratf. So the fix should be in taychk2rat, not in the callers. Or maybe in $taytorat. I don't understand what srrat does  is that the problem? Maybe look at this later.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=932302&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #817516, was opened at 20031003 19:27 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=817516&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(xxx)xxx incorrect Initial Comment: Define q: q:(LOG(n)n)/(LOG(n)1)$ Get a series representation around infinity to 1 term: qt: taylor(q,n,inf,1); ((1/LOG(n))+...)*n+(1+1/LOG(n)+...)+... OK so far. Now take the difference of the original q and the series: qdiff: qqt; (11/LOG(n)1/LOG(n)^21/LOG(n)^3 1/LOG(n)^4+...) * n + (1+1/LOG(n)+...)+... taylorinfo(qt) => [[1/LOG(n),ZEROA,1],[n,INF,1]] taylorinfo(qdiff) => [[1/LOG(n),ZEROA,4],[n,INF,4]] The result should have been the same as: qdiffr: taylor(qratsimp(qt),n,inf,1) which gives (almost correctly, cf bug #774065) ZEROA*N + (ZEROA * LOG(N)  ZEROA + ...) +... i.e. 0+... Huh? First of all, the answer is not correct. The initial terms should have cancelled. Secondly, why did it calculate three more terms?  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=817516&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #974734, was opened at 20040617 07:30 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=974734&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisbl) Assigned to: Nobody/Anonymous (nobody) Summary: derivatives of Taylor polynonials Initial Comment: Derivatives of Taylor polynomials sometimes are not simplified; an example (%i1) taylor(x/(1+x),x,inf,4); (%o1) 11/x+1/x^21/x^3+1/x^4 (%i2) diff(%,x); (%o2) 1/x^22/(x^2*x)+3/(x^2*x^2)4/(x^2*x^3) Notice the x^2 * x, x^2 * x^2, ... terms. Sometimes derivatives of Taylor polynomials work fine (%i3) taylor(x/(1+x),x,1,4); (%o3) 1/2+(x1)/4(x1)^2/8+(x1)^3/16(x1)^4/32 (%i4) diff(%,x); (%o4) 1/4(x1)/4+3*(x1)^2/16(x1)^3/8 %i5) build_info(); Maxima version: 5.9.0.1cvs Maxima build date: 12:16 5/28/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.7.0 Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=974734&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #939022, was opened at 20040420 20:37 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=939022&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: matrix(...taylor...)^^1 wrong Initial Comment: Consider m: matrix([taylor(1+a*x,x,0,1),0], [0,taylor(1+d*x,x,0,1)]); Now, m^^1 => matrix([1/(d*x+1),0],[0,1/(d*x+1)]) There are two problems with this. First, the answer is incorrect. Compare: matrixmap(ratdisrep,m)^^1 => matrix([1/(a*x+1),0],[0,1/(d*x+1)]) Second, the answer is not in terms of taylor series. A silent ratdisrep was done in the middle, losing truncation information.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=939022&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #774065, was opened at 20030718 23:36 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=774065&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor returns zeroa Initial Comment: taylor(x/log(x),x,inf,1) returns (zeroa+1/log(x)+...)*x+... I believe that 'zeroa' is intended to be strictly internal. Substituting 0 for zeroa gives the correct result.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=774065&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #734851, was opened at 20030508 14:50 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=734851&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: pade interfered with by taylor Initial Comment: (C1) taylor(sin(x),x,0,3); (D1) xx^3/6 (C2) pade(d1,2,2); (D2) [6*x/(x^2+6)] Fine so far. (C3) taylor(x,x,0,3); (D3) +x (C4) pade(d1,2,2); (D4) [6/7] What's this?! Pade is being affected by the Taylor calculation in C3?! Presumably there is some global flag being set but not reset.... (C5) taylor(exp(x),x,0,3); (D5) 1+x+x^2/2+x^3/6 (C6) pade(d1,2,2); (D6) [6*x/(x^2+6)] Somehow that restored the thing that was causing the problem.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=734851&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #807275, was opened at 20030916 11:19 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=807275&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: Taylor Illegal log kernel: log(cos(th)) @ %pi/2 Initial Comment: taylor(log(cos(th)),th,%pi/2,2) gives the internal error Illegal log kernel but taylor(log(sin(th)),th,0,2), which is equivalent, gives a perfectly reasonable result.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=807275&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #752332, was opened at 20030610 20:47 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=752332&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor: "sign called on imaginary arg Initial Comment: expr: atan(1/sqrt(1b))/sqrt(1b); taylor(expr,b,1,2) => SIGN called on an imaginary argument: %I (bizarrely, it says this twice....) I suspected that this was because of %i introduced by using logarc on expr internally, BUT: taylor(logarc(expr),b,1,2) => works just fine Maxima 5.9.0 gcl 2.5 Windows 2000  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=752332&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #1037916, was opened at 20040930 11:44 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1037916&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor 2^n/(7^n+1) n>inf internal error Initial Comment: taylor(2^n/(7^n+1),n,inf,3) => Maxima encountered a Lisp error: two equal vars generated 5.9.0.9beta2 Maxima build date: 10:50 7/27/2004 host type: i686pcmingw32 lispimplementationtype: Kyoto Common Lisp lispimplementationversion: GCL 2.6.3  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1037916&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #836780, was opened at 20031105 13:55 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=836780&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  Taylor 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 07: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 07: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 13: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...>  20060409 20:43:47

Bugs item #904522, was opened at 20040225 13:30 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=904522&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylorinfo fatal err in multivar case Initial Comment: taylorinfo(taylor(x,[x,y],0,1)) fatal error Maxima 5.9.0 gcl 2.5.0 mingw32 W2k Athlon  Comment By: Barton Willis (willisbl) Date: 20040226 08:12 Message: Logged In: YES user_id=895922 If you fix this bug, also fix (the easy to fix) bug 867310. Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=904522&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #1041148, was opened at 20041005 20:46 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1041148&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  Taylor Group: None Status: Open Resolution: None Priority: 3 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(x^a,[x],0,1) unsimplified Initial Comment: taylor(x^a,[x],0,1) => 1^2*x^a+... Note unsimplified 1^2.  Comment By: Barton Willis (willisbl) Date: 20041006 10:59 Message: Logged In: YES user_id=895922 Maybe there is a connection between this bug and bug 974734. Barton  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1041148&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #991622, was opened at 20040715 08:34 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=991622&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(taylor(exp(x^2)...)) errors Initial Comment: taylor(taylor(exp(x^2),x,inf,2),x,inf,2); Error: NIL is not of type CONS. Error signalled by PUSHPW. Also, bizarrely: taylor(taylor(exp(x^2)/sqrt(%pi),x,inf,2),x,inf,2) => 1/sqrt(%pi)+... which is incorrect, but not a fatal error.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=991622&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #1358239, was opened at 20051116 09:52 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1358239&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(exp(sqrt(log(x)*log(log(x))))@inf > stack overflow Initial Comment: Error in PROGN [or a callee]: Bind stack overflow. This is a different stack overflow than some of the other cases I've reported.  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1358239&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:47

Bugs item #769860, was opened at 20030711 14:49 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=769860&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor bind stack overflow: sin(2*atan(x))@q Initial Comment: taylor(sin(2*atan(x)),x,q,1) results in a bind stack overflow (internal error)  try expand appears to be in an infinite recursion. The expansion is wellbehaved, and can be derived using the TaylorMacLaurin formula: r1: taylor(f(x),x,q,1)$ r2: subst(lambda([x],sin(2*atan(x))),f,r1)$ r3: ev(r2,diff,at)$ Which gives: 2 COS(2 ATAN(q)) (x  q)  + SIN(2 ATAN(q)) 2 q + 1 This can be prettified: map(factor,trigexpand(r3)) ...=> 2 q 2 (q  1) (q + 1) (x  q)    2 2 2 q + 1 (q + 1)  Comment By: Stavros Macrakis (macrakis) Date: 20030711 18:46 Message: Logged In: YES user_id=588346 Untabified: 2 COS(2 ATAN(q)) (x  q)  + SIN(2 ATAN(q)) 2 q + 1 map(factor,trigexpand(r3)); 2 q 2 (q  1) (q + 1) (x  q)    2 2 2 q + 1 (q + 1)  Comment By: Stavros Macrakis (macrakis) Date: 20030711 18:44 Message: Logged In: YES user_id=588346 Did I forget to untabify? Or is sourceforge compressing spaces now? Testing: aaaTaaaTaaaTaaa aaa aaa aaa aaa aaa aaa aaa aaa aaaSSSSSaaaSSSSSaaaSSSSSaaa  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=769860&group_id=4933 
From: SourceForge.net <noreply@so...>  20060409 20:43:46

Bugs item #701628, was opened at 20030311 09:04 Message generated for change (Settings changed) made by robert_dodier You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=701628&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  Taylor Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: taylor(x^2,x,minf,2)  wrong sign Initial Comment: taylor(x^2,x,minf,2) => 1/(1/x)^2 + ... which is equal to x^2 instead of x^2.  Comment By: Wolfgang Jenkner (wjenkner) Date: 20030313 20:15 Message: Logged In: YES user_id=581700 Maybe I am logged in this time. I have concocted this patch. Wolfgang  Comment By: Nobody/Anonymous (nobody) Date: 20030313 19:53 Message: Logged In: NO What do you think about the following patch? Index: hayat.lisp =================================================================== RCS file: /cvsroot/maxima/maxima/src/hayat.lisp,v retrieving revision 1.2 diff C2 r1.2 hayat.lisp *** hayat.lisp 24 Feb 2001 02:21:59 0000 1.2  hayat.lisp 14 Mar 2003 01:48:45 0000 *************** *** 2325,2332 **** ((e> (setq exp (rcminus exp)) (currenttrunc temp)) (rczero)) (t (makeps (intvar temp) (ncons (if exactpoly (inf) (currenttrunc temp))) (ncons (term exp ! (if (eq (exppt temp) '$MINF) (rcmone) (rcone))))))))  2325,2343  ((e> (setq exp (rcminus exp)) (currenttrunc temp)) (rczero)) + ;; Say VAR is x and the original EXP is a=p/q with + ;; gcd(p,q)=1. Following the convention of taking + ;; real roots if available, we want to express x^a + ;; by (1/x)^a (for MINF, hence we can assume x<0) + ;; resp. (1/x)^a (for INF, hence we can assume + ;; x>0; INFINITY is treated like INF here). (t (makeps (intvar temp) (ncons (if exactpoly (inf) (currenttrunc temp))) (ncons (term exp ! (if (and (eq (exppt temp) '$MINF) ! ;; For odd q and x<0, ! ;; x^a=(1)^p/(1/x)^a ! (oddp (car exp)) ! (or (oddp (cdr exp)) ! (merror "Ambiguous root ~M in expansion at MINF." (m^t var (m//t 1 (cdr exp)))))) (rcmone) (rcone))))))))  You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=701628&group_id=4933 