You can subscribe to this list here.
2002 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}
(67) 
_{Jul}
(61) 
_{Aug}
(49) 
_{Sep}
(43) 
_{Oct}
(59) 
_{Nov}
(24) 
_{Dec}
(18) 

2003 
_{Jan}
(34) 
_{Feb}
(35) 
_{Mar}
(72) 
_{Apr}
(42) 
_{May}
(46) 
_{Jun}
(15) 
_{Jul}
(64) 
_{Aug}
(62) 
_{Sep}
(22) 
_{Oct}
(41) 
_{Nov}
(57) 
_{Dec}
(56) 
2004 
_{Jan}
(48) 
_{Feb}
(47) 
_{Mar}
(33) 
_{Apr}
(39) 
_{May}
(6) 
_{Jun}
(17) 
_{Jul}
(19) 
_{Aug}
(10) 
_{Sep}
(14) 
_{Oct}
(74) 
_{Nov}
(80) 
_{Dec}
(22) 
2005 
_{Jan}
(43) 
_{Feb}
(33) 
_{Mar}
(52) 
_{Apr}
(74) 
_{May}
(32) 
_{Jun}
(58) 
_{Jul}
(18) 
_{Aug}
(41) 
_{Sep}
(71) 
_{Oct}
(28) 
_{Nov}
(65) 
_{Dec}
(68) 
2006 
_{Jan}
(54) 
_{Feb}
(37) 
_{Mar}
(82) 
_{Apr}
(211) 
_{May}
(69) 
_{Jun}
(75) 
_{Jul}
(279) 
_{Aug}
(139) 
_{Sep}
(135) 
_{Oct}
(58) 
_{Nov}
(81) 
_{Dec}
(78) 
2007 
_{Jan}
(141) 
_{Feb}
(134) 
_{Mar}
(65) 
_{Apr}
(49) 
_{May}
(61) 
_{Jun}
(90) 
_{Jul}
(72) 
_{Aug}
(53) 
_{Sep}
(86) 
_{Oct}
(61) 
_{Nov}
(62) 
_{Dec}
(101) 
2008 
_{Jan}
(100) 
_{Feb}
(66) 
_{Mar}
(76) 
_{Apr}
(95) 
_{May}
(77) 
_{Jun}
(93) 
_{Jul}
(103) 
_{Aug}
(76) 
_{Sep}
(42) 
_{Oct}
(55) 
_{Nov}
(44) 
_{Dec}
(75) 
2009 
_{Jan}
(103) 
_{Feb}
(105) 
_{Mar}
(121) 
_{Apr}
(59) 
_{May}
(103) 
_{Jun}
(82) 
_{Jul}
(67) 
_{Aug}
(76) 
_{Sep}
(85) 
_{Oct}
(75) 
_{Nov}
(181) 
_{Dec}
(133) 
2010 
_{Jan}
(107) 
_{Feb}
(116) 
_{Mar}
(145) 
_{Apr}
(89) 
_{May}
(138) 
_{Jun}
(85) 
_{Jul}
(82) 
_{Aug}
(111) 
_{Sep}
(70) 
_{Oct}
(83) 
_{Nov}
(60) 
_{Dec}
(16) 
2011 
_{Jan}
(61) 
_{Feb}
(16) 
_{Mar}
(52) 
_{Apr}
(41) 
_{May}
(34) 
_{Jun}
(41) 
_{Jul}
(57) 
_{Aug}
(73) 
_{Sep}
(21) 
_{Oct}
(45) 
_{Nov}
(50) 
_{Dec}
(28) 
2012 
_{Jan}
(70) 
_{Feb}
(36) 
_{Mar}
(71) 
_{Apr}
(29) 
_{May}
(48) 
_{Jun}
(61) 
_{Jul}
(44) 
_{Aug}
(54) 
_{Sep}
(20) 
_{Oct}
(28) 
_{Nov}
(41) 
_{Dec}
(137) 
2013 
_{Jan}
(62) 
_{Feb}
(55) 
_{Mar}
(31) 
_{Apr}
(23) 
_{May}
(54) 
_{Jun}
(54) 
_{Jul}
(90) 
_{Aug}
(46) 
_{Sep}
(38) 
_{Oct}
(60) 
_{Nov}
(92) 
_{Dec}
(17) 
2014 
_{Jan}
(62) 
_{Feb}
(35) 
_{Mar}
(72) 
_{Apr}
(30) 
_{May}
(97) 
_{Jun}
(81) 
_{Jul}
(63) 
_{Aug}
(64) 
_{Sep}
(28) 
_{Oct}
(45) 
_{Nov}
(48) 
_{Dec}
(109) 
2015 
_{Jan}
(106) 
_{Feb}
(36) 
_{Mar}
(65) 
_{Apr}
(63) 
_{May}
(95) 
_{Jun}
(56) 
_{Jul}
(48) 
_{Aug}
(55) 
_{Sep}
(100) 
_{Oct}
(57) 
_{Nov}
(33) 
_{Dec}
(46) 
2016 
_{Jan}
(76) 
_{Feb}
(53) 
_{Mar}
(88) 
_{Apr}
(79) 
_{May}
(62) 
_{Jun}
(65) 
_{Jul}
(37) 
_{Aug}
(23) 
_{Sep}
(108) 
_{Oct}
(68) 
_{Nov}
(66) 
_{Dec}
(47) 
2017 
_{Jan}
(55) 
_{Feb}
(8) 
_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

From: Robert Dodier <robert_dodier@us...>  20170225 06:00:45

I'm not sure this is a bug. 0 has the appropriate property that `x + 0` is `x` even when `x` is a list. If you change your interpretation of 0 to mean "the additive identity", does anything work incorrectly? Are there algebraic identities which do not hold?  ** [bugs:#3287] The cross product is wrong when I use the same vector** **Status:** open **Group:** None **Created:** Mon Feb 20, 2017 05:11 PM UTC by Daniel Volinski **Last Updated:** Mon Feb 20, 2017 05:11 PM UTC **Owner:** nobody The cross product is wrong when I use the same vector: load(vect)$ express([0,1,1]~[0,1,2]);=> [1,0,0] /* Correct! */ express([0,1,1]~[0,1,1]);=> 0 /* Incorrect, it should be a vector! */ Daniel  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170223 20:33:56

It appears from the log file you posted that wxMaxima has crashed (with a segfault). Does command line Maxima work OK?  ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** open **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Tue Feb 21, 2017 06:45 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Van <lvwarren@us...>  20170221 06:46:00

After each Maxima crash, a common lisp process 'sbcl' is left running on the machine and this process has very high CPU utilization. If not detected several of these processes taken together can overheat some machines, Stuxnet style. I have to keep on an icepack to keep the CPU temperature below 210 degrees. Bringing up the activity manager or doing a ps au shows the zombie common lisps which can be killed one at a time.  ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** open **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Tue Feb 21, 2017 06:32 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Van <lvwarren@us...>  20170221 06:32:32

 ** [bugs:#3288] Maxima is dead in the water on the MacOS after Apple security update.** **Status:** open **Group:** None **Created:** Tue Feb 21, 2017 06:32 AM UTC by Van **Last Updated:** Tue Feb 21, 2017 06:32 AM UTC **Owner:** nobody **Attachments:**  [MacOSMaximaCrash.rtf](https://sourceforge.net/p/maxima/bugs/3288/attachment/MacOSMaximaCrash.rtf) (63.8 kB; text/rtf)  [stringLoad.wxm](https://sourceforge.net/p/maxima/bugs/3288/attachment/stringLoad.wxm) (348 Bytes; application/octetstream) I use Maxima almost every day. As of a week ago Maxima is dead in the water on the MacOS Platform running El Capitan. I have redownloaded most recent version, scrubbed and cleaned and I get a reproducible Maxima crash with the attached error log. and .wxm file that contains only two lines kill(all)$ and 2+2; All other codes on this machine function properly. It is also stable otherwise. I downloaded Apple security updates and the problem started around the time of that update. Recent versions of MacOS have protected certain directories that were formerly available for install. My instinct is that some access to these formerly open directories is the cause of the crash.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Daniel Volinski <danielvolinski@us...>  20170220 17:11:09

 ** [bugs:#3287] The cross product is wrong when I use the same vector** **Status:** open **Group:** None **Created:** Mon Feb 20, 2017 05:11 PM UTC by Daniel Volinski **Last Updated:** Mon Feb 20, 2017 05:11 PM UTC **Owner:** nobody The cross product is wrong when I use the same vector: load(vect)$ express([0,1,1]~[0,1,2]);=> [1,0,0] /* Correct! */ express([0,1,1]~[0,1,1]);=> 0 /* Incorrect, it should be a vector! */ Daniel  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Gilles Schintgen <gschint@us...>  20170205 18:33:35

 ** [bugs:#3286] implicit_plot treats some discontinuities as sign changes** **Status:** open **Group:** None **Created:** Sun Feb 05, 2017 06:33 PM UTC by Gilles Schintgen **Last Updated:** Sun Feb 05, 2017 06:33 PM UTC **Owner:** nobody It seems that implicit_plot treats a change from positive to complex as a sign change. Here's a simple example: load(implicit_plot)$ implicit_plot(1/sqrt(x^2+y^21),[x,2,2],[y,2,2])$ Result: the unit circle is drawn even though the expression to be plotted is either positive or nonreal. It's not even approaching zero when (x,y) is close to the unit circle. A cursory look at the code seems to confirm that implicit_plot tries to detect sign changes but is actually looking for changes from larger than ip_epsilon to "everything else" (in function sampledata).  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: David Scherfgen <tomasriker@us...>  20170203 08:05:15

 Description has changed: Diff: ~~~~  old +++ new @@ 1,4 +1,4 @@ `radcan` fails to simplify an expression to zero that contains no variables, only logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: +`radcan` fails to simplify an expression to zero that contains no variables, only exponentials and logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: ~~~ (%i1) X : 2^(log(3)/(log(2)log(3)))3^(log(3)/(log(2)log(3))+1)$ (%i2) radcan(X); ~~~~  ** [bugs:#3285] radcan: failure to simplify to zero depending on slight change and/or ?expsumsplit** **Status:** open **Group:** None **Labels:** radcan **Created:** Fri Feb 03, 2017 08:04 AM UTC by David Scherfgen **Last Updated:** Fri Feb 03, 2017 08:04 AM UTC **Owner:** nobody `radcan` fails to simplify an expression to zero that contains no variables, only exponentials and logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: ~~~ (%i1) X : 2^(log(3)/(log(2)log(3)))3^(log(3)/(log(2)log(3))+1)$ (%i2) radcan(X); (%o2) (%e^(log(3)^2/(log(3)log(2)))*(%e^(log(3)^2/(log(3)log(2)))3*2^(log(3)/(log(3)log(2)))))/2^(log(3)/(log(3)log(2))) (%i3) radcan(logcontract(X)); (%o3) 0 (%i4) Y : 2^(log(3)/(log(2)log(3)))1/3^(log(2)/(log(3)log(2)))$ (%i5) radcan(Y); (%o5) 0 (%i6) radcan(X), ?expsumsplit:false; (%o6) 0 ~~~ Thanks to Richard Fateman for investigating.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: David Scherfgen <tomasriker@us...>  20170203 08:04:37

 ** [bugs:#3285] radcan: failure to simplify to zero depending on slight change and/or ?expsumsplit** **Status:** open **Group:** None **Labels:** radcan **Created:** Fri Feb 03, 2017 08:04 AM UTC by David Scherfgen **Last Updated:** Fri Feb 03, 2017 08:04 AM UTC **Owner:** nobody `radcan` fails to simplify an expression to zero that contains no variables, only logarithms of natural numbers. When we slightly rewrite the expression to an equivalent one, `radcan` succeeds. Applying `logcontract` also helps, as does setting the internal flag `?expsumsplit` to `false`: ~~~ (%i1) X : 2^(log(3)/(log(2)log(3)))3^(log(3)/(log(2)log(3))+1)$ (%i2) radcan(X); (%o2) (%e^(log(3)^2/(log(3)log(2)))*(%e^(log(3)^2/(log(3)log(2)))3*2^(log(3)/(log(3)log(2)))))/2^(log(3)/(log(3)log(2))) (%i3) radcan(logcontract(X)); (%o3) 0 (%i4) Y : 2^(log(3)/(log(2)log(3)))1/3^(log(2)/(log(3)log(2)))$ (%i5) radcan(Y); (%o5) 0 (%i6) radcan(X), ?expsumsplit:false; (%o6) 0 ~~~ Thanks to Richard Fateman for investigating.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Jinsong Zhao <jszhao@us...>  20170131 09:59:05

To the issue 4. I have wrote a Lisp script to replace the python script used during compiling the document. I attach the extract_categories1.awk and extract_categories1.sed here, FYI. The only issue with this approach is I don't know how to set correct Lisp program name in extract_categories.sh. ~~~  > tmpmakecategories.py + > tmpmakecategories.lisp ${PYTHONBIN:python} tmpmakecategories.py +#echo ${SBCL_NAME} +/d/Program/Steel\ Bank\ Common\ Lisp/1.3.10/sbcl.exe script tmpmakecategories.lisp ~~~ I meant here, now, it's hard coded in the script, it's not the right way. Best regards, Jinsong Attachments:  [extract_categories1.awk](https://sourceforge.net/p/maxima/bugs/_discuss/thread/9d22061f/886c/attachment/extract_categories1.awk) (1.3 kB; application/octetstream)  [extract_categories1.sed](https://sourceforge.net/p/maxima/bugs/_discuss/thread/9d22061f/886c/attachment/extract_categories1.sed) (617 Bytes; application/octetstream)  ** [bugs:#3284] maxima build system** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 03:04 PM UTC by Jinsong Zhao **Last Updated:** Sun Jan 29, 2017 03:04 PM UTC **Owner:** nobody There are several issues that related with the maxima building system on Windows platform. I don't know if there are bugs, however, it makes trouble or confusion when I try to make a maxima installer for Windows. I posted them in the maximadiscuss mailing list. 1) After './configure...' and 'make', 'make install' gives the following error message: /usr/bin/install: cannot create regular file 'C:/maxima5.39.0sbcl/bin/xmaxima': File exists 'xmaxima' in interface/xmaxima seems to be a script, and not need by XMaxima on Windows platform. On Windows, 'xmaxima.exe' is the correct executable file. This error prevents the running of make install. I have to changed the Makefile.in in interface/xmaxima from $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  exit $$?; \ to $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  echo "xxx" $$?; \ according to https://sourceforge.net/projects/maxima/files/MaximaWindows/5.XX.XWindows%20Build%20System/README.md Then make install could be run successfully. If possible, could the 'xmaxima' script be removed from makefile for Windows platform? 2) Could the 'xmaxima' script be removed from the 'maxima.iss' on Windows platform? 3) There appears a directory named ".cirtcl" outside the maxima toplevel directory during the making process. In the directory, there are two files, i.e., mgnuplot.bat and xmaixma.bat. I don't know if they are needed, and why there are be put outside the toplevel directory. It seems they are related with xmaxima interface. 4) Could we not compile the document again? In the doc/, there are correct html and info files, why we generate them again? I encountered the most difficult thing is the python script 'tmpmakecategories.py does not work in MSYS2 system. I don't know why. So, if possible, could configure provide a option that disable doc generation. Best regards, Jinsong  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170130 01:48:23

 **labels**: gamma, gamma_greek, gamma_incomp > gamma, gamma_greek, gamma_incomplete  ** [bugs:#3277] no numerical evaluation for gamma_greek** **Status:** open **Group:** None **Labels:** gamma gamma_greek gamma_incomplete **Created:** Thu Jan 19, 2017 12:50 AM UTC by kofoid **Last Updated:** Mon Jan 30, 2017 01:47 AM UTC **Owner:** nobody I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170130 01:47:21

 **labels**: > gamma, gamma_greek, gamma_incomp  **summary**: gammarelated functions under OS X > no numerical evaluation for gamma_greek  **Comment**: Adjusted title to focus on gamma_greek and omit OS.  ** [bugs:#3277] no numerical evaluation for gamma_greek** **Status:** open **Group:** None **Labels:** gamma gamma_greek gamma_incomp **Created:** Thu Jan 19, 2017 12:50 AM UTC by kofoid **Last Updated:** Sun Jan 22, 2017 10:23 PM UTC **Owner:** nobody I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170129 20:04:39

 **labels**: > lsquares, share  **status**: open > closed  **Comment**: Fixed by commit b58a547. Thanks for the pointing it out. Closing this report as fixed.  ** [bugs:#3282] lsquares and lists of list of data to be fitted on** **Status:** closed **Group:** None **Labels:** lsquares share **Created:** Tue Jan 24, 2017 06:58 PM UTC by Gunter Königsmann **Last Updated:** Thu Jan 26, 2017 11:46 PM UTC **Owner:** nobody I have a list of data sets that are huge => running lsquares_estimates() on them to fit them on an equation would need too big amounts of time. Normally lsquares_estimates_approximate() helps in this case as it skips the step that tries to find an exact solution. But as I have not a single data set, but a list of data sets that somehow backfires: ~~~ (%i34) load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; lsquares_mse(data[1],[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o30) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) [matrix( [0, 0], [1, 1], [2, 2] ),matrix( [2, 2], [1, 1], [2, 2] )] (eqtn) y=a*x^2+b*x+c (%o33) sum((m6[i,2]a*m6[i,1]^2b*m6[i,1]c)^2,i,1,1) Maxima encountered a Lisp error: The value ((MTIMES SIMP) 2.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2) ((MPLUS SIMP) 1.0 ((MTIMES SIMP) 1.0 ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1)) ((MTIMES SIMP) 1.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2)) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 2))) is not of type DOUBLEFLOAT. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. ~~~ Only running lsquares_estimates on a single entry of the list works: ~~~ (%i29) load("lsquares"); data:matrix([0,0],[1,1],[2,2]); eqtn:y=a*x^2+b*x+c; lsquares_mse(data,[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o25) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) matrix( [0, 0], [1, 1], [2, 2] ) (eqtn) y=a*x^2+b*x+c (%o28) sum((data[i,2]a*data[i,1]^2b*data[i,1]c)^2,i,1,3)/3 ************************************************* N= 3 NUMBER OF CORRECTIONS=25 INITIAL VALUES F= 1.000000000000000D+01 GNORM= 1.753726191728787D+01 ************************************************* I NFN FUNC GNORM STEPLENGTH 1 2 2.796937667198332D01 1.961076084276657D+00 5.702144409522791D02 2 3 1.464269617620298D01 4.983306664572200D01 1.000000000000000D+00 ~~~ snip ~~~ 13 14 3.397773998012277D06 4.742450181229817D03 1.000000000000000D+00 14 15 2.643557109421608D08 2.597101638571121D04 1.000000000000000D+00 THE MINIMIZATION TERMINATED WITHOUT DETECTING ERRORS. IFLAG = 0 (%o29) [[a=3.405898581949557*10^4,b=1.000676542954337,c=1.342533507359797*10^4]] ~~~ I have been able to ship around this problem by putting the data set I am working on into a temporary variable before running lsquares: ~~~ load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; makelist(block([tmp], tmp:data[i], lsquares_estimates_approximate(lsquares_mse(tmp,[x,y],eqtn),[a,b,c]) ), i,1,length(data) ); ...but I am convinced that the fact that I needed this is a bug. ~~~  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Stavros Macrakis <macrakis@us...>  20170129 19:48:03

More, simpler, examples: taylor(rat(1/z+1),z,0,1) => unfamiliar singularity error taylor( 1/z+1,z,0,1) => 1/z+1 taylor(rat(z+1),z,inf,0) => unfamiliar singularity error taylor( z+1 ,z,inf,0) => z+1 Maxima branch_5_37_base_528_g802ae7f using Lisp SBCL 1.2.11  ** [bugs:#3068] taylor of CRE fails** **Status:** open **Group:** None **Created:** Tue Jan 12, 2016 07:50 PM UTC by Stavros Macrakis **Last Updated:** Tue Jan 12, 2016 07:50 PM UTC **Owner:** nobody taylor(rat(x/(1+x),x,inf,2) => taylor: encountered an unfamiliar singularity though taylor((x/(1+x),x,inf,2) => 1  1/x + 1/x^2 + ...  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Jinsong Zhao <jszhao@us...>  20170129 15:04:46

 ** [bugs:#3284] maxima build system** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 03:04 PM UTC by Jinsong Zhao **Last Updated:** Sun Jan 29, 2017 03:04 PM UTC **Owner:** nobody There are several issues that related with the maxima building system on Windows platform. I don't know if there are bugs, however, it makes trouble or confusion when I try to make a maxima installer for Windows. I posted them in the maximadiscuss mailing list. 1) After './configure...' and 'make', 'make install' gives the following error message: /usr/bin/install: cannot create regular file 'C:/maxima5.39.0sbcl/bin/xmaxima': File exists 'xmaxima' in interface/xmaxima seems to be a script, and not need by XMaxima on Windows platform. On Windows, 'xmaxima.exe' is the correct executable file. This error prevents the running of make install. I have to changed the Makefile.in in interface/xmaxima from $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  exit $$?; \ to $(INSTALL_SCRIPT) $$files "$(DESTDIR)$(bindir)$$dir"  echo "xxx" $$?; \ according to https://sourceforge.net/projects/maxima/files/MaximaWindows/5.XX.XWindows%20Build%20System/README.md Then make install could be run successfully. If possible, could the 'xmaxima' script be removed from makefile for Windows platform? 2) Could the 'xmaxima' script be removed from the 'maxima.iss' on Windows platform? 3) There appears a directory named ".cirtcl" outside the maxima toplevel directory during the making process. In the directory, there are two files, i.e., mgnuplot.bat and xmaixma.bat. I don't know if they are needed, and why there are be put outside the toplevel directory. It seems they are related with xmaxima interface. 4) Could we not compile the document again? In the doc/, there are correct html and info files, why we generate them again? I encountered the most difficult thing is the python script 'tmpmakecategories.py does not work in MSYS2 system. I don't know why. So, if possible, could configure provide a option that disable doc generation. Best regards, Jinsong  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Dima Pasechnik <dimpase@us...>  20170129 11:49:27

Here is how to get this with Maxima only (5.39.0 with SBCL or ECL): ~~~ s1: 2*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)*(cos(z)*sin(y) + cos(y)*sin(z))*((cos(1/2*x)*cos(z)  cos(y)*cos(z) + sin(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)*cos(z)*sin(y) + cos(y)^2*sin(z)  (cos(z)*sin(y) + cos(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))/((cos(1/2*x)*cos(y)  1)*((cos(y)*cos(z)  sin(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)*cos(y)  1) + (cos(z)*sin(y) + cos(y)*sin(z))*((cos(1/2*x)*cos(z)  cos(y)*cos(z) + sin(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)*cos(z)*sin(y) + cos(y)^2*sin(z)  (cos(z)*sin(y) + cos(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))*sin(1/2*x)/((cos(1/2*x)*cos(y)  1)*((cos(z)*sin(y)  cos(1/2*x)*sin(z) + cos(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)^2*cos(z)  cos(y)*sin(y)*sin(z)  (cos(y)*cos(z)  sin(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))))*((cos(z)*sin(y)  cos(1/2*x)*sin(z) + cos(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)^2*cos(z)  cos(y)*sin(y)*sin(z)  (cos(y)*cos(z)  sin(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))^2*sin(y))  2*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/(((cos(z)*sin(y)  cos(1/2*x)*sin(z) + cos(y)*sin(z))*sin(1/2*x)/(cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1) + (cos(y)^2*cos(z)  cos(y)*sin(y)*sin(z)  (cos(y)*cos(z)  sin(y)*sin(z))*cos(1/2*x))*(cos(1/2*x)*cos(y) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(y)^2  2*cos(1/2*x)*cos(y) + 1)*sin(1/2*x)*sin(y)))*sin(1/2*x)*sin(y))$ trigsimp(s1); s2: 2*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)*(cos(b)*sin(a) + cos(a)*sin(b))*((cos(1/2*x)*cos(b)  cos(a)*cos(b) + sin(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)*cos(b)*sin(a) + cos(a)^2*sin(b)  (cos(b)*sin(a) + cos(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))/((cos(1/2*x)*cos(a)  1)*((cos(a)*cos(b)  sin(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)*cos(a)  1) + (cos(b)*sin(a) + cos(a)*sin(b))*((cos(1/2*x)*cos(b)  cos(a)*cos(b) + sin(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)*cos(b)*sin(a) + cos(a)^2*sin(b)  (cos(b)*sin(a) + cos(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))*sin(1/2*x)/((cos(1/2*x)*cos(a)  1)*((cos(b)*sin(a)  cos(1/2*x)*sin(b) + cos(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)^2*cos(b)  cos(a)*sin(a)*sin(b)  (cos(a)*cos(b)  sin(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))))*((cos(b)*sin(a)  cos(1/2*x)*sin(b) + cos(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)^2*cos(b)  cos(a)*sin(a)*sin(b)  (cos(a)*cos(b)  sin(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))^2*sin(a))  2*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/(((cos(b)*sin(a)  cos(1/2*x)*sin(b) + cos(a)*sin(b))*sin(1/2*x)/(cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1) + (cos(a)^2*cos(b)  cos(a)*sin(a)*sin(b)  (cos(a)*cos(b)  sin(a)*sin(b))*cos(1/2*x))*(cos(1/2*x)*cos(a) + sin(1/2*x)^2  1)/((cos(1/2*x)^2*cos(a)^2  2*cos(1/2*x)*cos(a) + 1)*sin(1/2*x)*sin(a)))*sin(1/2*x)*sin(a))$ trigsimp(s2); ~~~ gives ~~~ Maxima 5.39.0 http://maxima.sourceforge.net using Lisp SBCL 1.3.12 Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) (%i2) x (%o2) (2 cos(y)  2 cos()) sin(y) sin(z) 2 x 2 + (2 cos() cos(y)  2 cos (y)) cos(z) 2 (%i3) (%i4) Quotient by a polynomial of higher degree  an error. To debug this try: debugmode(true); ~~~ Expression s2 is obtained from s1 by renaming y>a, z>b.  ** [bugs:#3283] trigsimp fails for some order of the variables** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 10:26 AM UTC by Eric Gourgoulhon **Last Updated:** Sun Jan 29, 2017 10:26 AM UTC **Owner:** nobody Hi, trigsimp may fail depending on the (alphabetical?) order of the variables; please see this SageMath bug report: https://groups.google.com/forum/#!topic/sagedevel/bJBgJhvRX4  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Eric Gourgoulhon <egourgoulhon@us...>  20170129 10:26:26

 ** [bugs:#3283] trigsimp fails for some order of the variables** **Status:** open **Group:** None **Created:** Sun Jan 29, 2017 10:26 AM UTC by Eric Gourgoulhon **Last Updated:** Sun Jan 29, 2017 10:26 AM UTC **Owner:** nobody Hi, trigsimp may fail depending on the (alphabetical?) order of the variables; please see this SageMath bug report: https://groups.google.com/forum/#!topic/sagedevel/bJBgJhvRX4  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170126 23:46:18

Yes, this is a bug in `lsquares_mse`. The problem is that it introduces a temporary variable and tries to assign the data to the temporary variable, but it is incorrect. I think I see how to fix it.  ** [bugs:#3282] lsquares and lists of list of data to be fitted on** **Status:** open **Group:** None **Created:** Tue Jan 24, 2017 06:58 PM UTC by Gunter Königsmann **Last Updated:** Tue Jan 24, 2017 06:58 PM UTC **Owner:** nobody I have a list of data sets that are huge => running lsquares_estimates() on them to fit them on an equation would need too big amounts of time. Normally lsquares_estimates_approximate() helps in this case as it skips the step that tries to find an exact solution. But as I have not a single data set, but a list of data sets that somehow backfires: ~~~ (%i34) load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; lsquares_mse(data[1],[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o30) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) [matrix( [0, 0], [1, 1], [2, 2] ),matrix( [2, 2], [1, 1], [2, 2] )] (eqtn) y=a*x^2+b*x+c (%o33) sum((m6[i,2]a*m6[i,1]^2b*m6[i,1]c)^2,i,1,1) Maxima encountered a Lisp error: The value ((MTIMES SIMP) 2.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2) ((MPLUS SIMP) 1.0 ((MTIMES SIMP) 1.0 ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1)) ((MTIMES SIMP) 1.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2)) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 2))) is not of type DOUBLEFLOAT. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. ~~~ Only running lsquares_estimates on a single entry of the list works: ~~~ (%i29) load("lsquares"); data:matrix([0,0],[1,1],[2,2]); eqtn:y=a*x^2+b*x+c; lsquares_mse(data,[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o25) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) matrix( [0, 0], [1, 1], [2, 2] ) (eqtn) y=a*x^2+b*x+c (%o28) sum((data[i,2]a*data[i,1]^2b*data[i,1]c)^2,i,1,3)/3 ************************************************* N= 3 NUMBER OF CORRECTIONS=25 INITIAL VALUES F= 1.000000000000000D+01 GNORM= 1.753726191728787D+01 ************************************************* I NFN FUNC GNORM STEPLENGTH 1 2 2.796937667198332D01 1.961076084276657D+00 5.702144409522791D02 2 3 1.464269617620298D01 4.983306664572200D01 1.000000000000000D+00 ~~~ snip ~~~ 13 14 3.397773998012277D06 4.742450181229817D03 1.000000000000000D+00 14 15 2.643557109421608D08 2.597101638571121D04 1.000000000000000D+00 THE MINIMIZATION TERMINATED WITHOUT DETECTING ERRORS. IFLAG = 0 (%o29) [[a=3.405898581949557*10^4,b=1.000676542954337,c=1.342533507359797*10^4]] ~~~ I have been able to ship around this problem by putting the data set I am working on into a temporary variable before running lsquares: ~~~ load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; makelist(block([tmp], tmp:data[i], lsquares_estimates_approximate(lsquares_mse(tmp,[x,y],eqtn),[a,b,c]) ), i,1,length(data) ); ...but I am convinced that the fact that I needed this is a bug. ~~~  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Peter Gasparovic <buteo@us...>  20170126 10:03:36

I confirm, that your example is correctly solved in Maxima 5.21.1. Problem with solve inconsistency is present with newest Maxima 5.38.1. I have the same problem with: eq: y=(x^2/2+1)^k; sol_a1: solve([eq], [x]); (I.e. correctly solved in 5.21.1, problem in 5.38.1.) In new version you can get the correct result by repeating solve sommand. In my case the first solve command works when I assign as a first command: simp:false; However that doesn't work on your example. I described it on Stackexchange: [Why “simp:false” helps Maxima to solve equation?](http://math.stackexchange.com/questions/2106603/whysimpfalsehelpsmaximatosolveequation)  ** [bugs:#3070] solve() inconsistency** **Status:** open **Group:** None **Created:** Wed Jan 13, 2016 11:35 AM UTC by Jaime E. Villate **Last Updated:** Thu Jan 14, 2016 08:19 PM UTC **Owner:** nobody The result of the following two commands should be the same: ~~~~ (%i2) solve(sqrt((x^2+1)/2) = 1); (%o2) [x = 1,x = 1] (%i3) solve(sqrt(x^2/2+1/2) = 1); (%o3) [sqrt(x^2+1)/sqrt(2) = 1] ~~~~ %o3 used to be the same as %o2 in older versions (at least in 5.21).  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Gunter Königsmann <peterpall@us...>  20170124 18:58:05

 ** [bugs:#3282] lsquares and lists of list of data to be fitted on** **Status:** open **Group:** None **Created:** Tue Jan 24, 2017 06:58 PM UTC by Gunter Königsmann **Last Updated:** Tue Jan 24, 2017 06:58 PM UTC **Owner:** nobody I have a list of data sets that are huge => running lsquares_estimates() on them to fit them on an equation would need too big amounts of time. Normally lsquares_estimates_approximate() helps in this case as it skips the step that tries to find an exact solution. But as I have not a single data set, but a list of data sets that somehow backfires: ~~~ (%i34) load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; lsquares_mse(data[1],[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o30) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) [matrix( [0, 0], [1, 1], [2, 2] ),matrix( [2, 2], [1, 1], [2, 2] )] (eqtn) y=a*x^2+b*x+c (%o33) sum((m6[i,2]a*m6[i,1]^2b*m6[i,1]c)^2,i,1,1) Maxima encountered a Lisp error: The value ((MTIMES SIMP) 2.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2) ((MPLUS SIMP) 1.0 ((MTIMES SIMP) 1.0 ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1)) ((MTIMES SIMP) 1.0 ((MEXPT SIMP) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 1) 2)) ((MQAPPLY SIMP ARRAY) (($DATA SIMP ARRAY) 1) 1 2))) is not of type DOUBLEFLOAT. Automatically continuing. To enable the Lisp debugger set *debuggerhook* to nil. ~~~ Only running lsquares_estimates on a single entry of the list works: ~~~ (%i29) load("lsquares"); data:matrix([0,0],[1,1],[2,2]); eqtn:y=a*x^2+b*x+c; lsquares_mse(data,[x,y],eqtn); lsquares_estimates_approximate(%,[a,b,c]); (%o25) "/usr/local/share/maxima/branch_5_39_base_66_gbbb452f/share/lsquares/lsquares.mac" (data) matrix( [0, 0], [1, 1], [2, 2] ) (eqtn) y=a*x^2+b*x+c (%o28) sum((data[i,2]a*data[i,1]^2b*data[i,1]c)^2,i,1,3)/3 ************************************************* N= 3 NUMBER OF CORRECTIONS=25 INITIAL VALUES F= 1.000000000000000D+01 GNORM= 1.753726191728787D+01 ************************************************* I NFN FUNC GNORM STEPLENGTH 1 2 2.796937667198332D01 1.961076084276657D+00 5.702144409522791D02 2 3 1.464269617620298D01 4.983306664572200D01 1.000000000000000D+00 ~~~ snip ~~~ 13 14 3.397773998012277D06 4.742450181229817D03 1.000000000000000D+00 14 15 2.643557109421608D08 2.597101638571121D04 1.000000000000000D+00 THE MINIMIZATION TERMINATED WITHOUT DETECTING ERRORS. IFLAG = 0 (%o29) [[a=3.405898581949557*10^4,b=1.000676542954337,c=1.342533507359797*10^4]] ~~~ I have been able to ship around this problem by putting the data set I am working on into a temporary variable before running lsquares: ~~~ load("lsquares"); data:[ matrix([0,0],[1,1],[2,2]), matrix([2,2],[1,1],[2,2]) ]; eqtn:y=a*x^2+b*x+c; makelist(block([tmp], tmp:data[i], lsquares_estimates_approximate(lsquares_mse(tmp,[x,y],eqtn),[a,b,c]) ), i,1,length(data) ); ...but I am convinced that the fact that I needed this is a bug. ~~~  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: David Scherfgen <tomasriker@us...>  20170124 17:45:39

Another interesting observation: With `domain : complex`, the integral from `a` to `0` is `(asin(2/(sqrt(4/a^2)*a))*a)/2`, while the integral from `0` to `a` is `(%pi*a)/4`.  ** [bugs:#3281] Exchanging limits in integral of sqrt(1(x/a)^2) does not always result in opposite sign** **Status:** open **Group:** None **Labels:** integrate defint **Created:** Mon Jan 23, 2017 08:30 PM UTC by Torben **Last Updated:** Tue Jan 24, 2017 05:43 PM UTC **Owner:** nobody The resultes below should be of opposite sign. ~~~~ Maxima 5.38.1 http://maxima.sourceforge.net using Lisp CLISP 2.49 (20100707) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) integrate(sqrt(1(x/a)^2),x,0,a); Is a positive or negative? pos; %pi a (%o1)  4 (%i2) integrate(sqrt(1(x/a)^2),x,a,0); Is a positive or negative? pos; %pi a (%o2)  4 (%i3) bug_report(); Please report bugs to: http://sourceforge.net/p/maxima/bugs To report a bug, you must have a Sourceforge account. Please include the following information with your bug report:  Maxima version: "5.38.1" Maxima build date: "20161019 00:27:24" Host type: "x86_64suselinuxgnu" Lisp implementation type: "CLISP" Lisp implementation version: "2.49 (20100707) (built on cloud112 [127.0.0.1])"  The above information is also reported by the function 'build_info()'. (%o3) ~~~~ Regards Torben  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: David Scherfgen <tomasriker@us...>  20170124 17:43:36

Interestingly, if the answer is `neg` instead of `pos`, the results differ in sign as expected.  ** [bugs:#3281] Exchanging limits in integral of sqrt(1(x/a)^2) does not always result in opposite sign** **Status:** open **Group:** None **Labels:** integrate defint **Created:** Mon Jan 23, 2017 08:30 PM UTC by Torben **Last Updated:** Tue Jan 24, 2017 03:44 PM UTC **Owner:** nobody The resultes below should be of opposite sign. ~~~~ Maxima 5.38.1 http://maxima.sourceforge.net using Lisp CLISP 2.49 (20100707) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) integrate(sqrt(1(x/a)^2),x,0,a); Is a positive or negative? pos; %pi a (%o1)  4 (%i2) integrate(sqrt(1(x/a)^2),x,a,0); Is a positive or negative? pos; %pi a (%o2)  4 (%i3) bug_report(); Please report bugs to: http://sourceforge.net/p/maxima/bugs To report a bug, you must have a Sourceforge account. Please include the following information with your bug report:  Maxima version: "5.38.1" Maxima build date: "20161019 00:27:24" Host type: "x86_64suselinuxgnu" Lisp implementation type: "CLISP" Lisp implementation version: "2.49 (20100707) (built on cloud112 [127.0.0.1])"  The above information is also reported by the function 'build_info()'. (%o3) ~~~~ Regards Torben  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170124 15:44:58

 **labels**: > integrate, defint  Description has changed: Diff: ~~~~  old +++ new @@ 1,5 +1,5 @@ The resultes below should be of opposite sign.  +~~~~ Maxima 5.38.1 http://maxima.sourceforge.net using Lisp CLISP 2.49 (20100707) Distributed under the GNU Public License. See the file COPYING. @@ 35,7 +35,7 @@ The above information is also reported by the function 'build_info()'. (%o3)  +~~~~ Regards Torben ~~~~  ** [bugs:#3281] Exchanging limits in integral of sqrt(1(x/a)^2) does not always result in opposite sign** **Status:** open **Group:** None **Labels:** integrate defint **Created:** Mon Jan 23, 2017 08:30 PM UTC by Torben **Last Updated:** Mon Jan 23, 2017 08:30 PM UTC **Owner:** nobody The resultes below should be of opposite sign. ~~~~ Maxima 5.38.1 http://maxima.sourceforge.net using Lisp CLISP 2.49 (20100707) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) integrate(sqrt(1(x/a)^2),x,0,a); Is a positive or negative? pos; %pi a (%o1)  4 (%i2) integrate(sqrt(1(x/a)^2),x,a,0); Is a positive or negative? pos; %pi a (%o2)  4 (%i3) bug_report(); Please report bugs to: http://sourceforge.net/p/maxima/bugs To report a bug, you must have a Sourceforge account. Please include the following information with your bug report:  Maxima version: "5.38.1" Maxima build date: "20161019 00:27:24" Host type: "x86_64suselinuxgnu" Lisp implementation type: "CLISP" Lisp implementation version: "2.49 (20100707) (built on cloud112 [127.0.0.1])"  The above information is also reported by the function 'build_info()'. (%o3) ~~~~ Regards Torben  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Torben <torhans1971@us...>  20170123 20:30:54

 ** [bugs:#3281] Exchanging limits in integral of sqrt(1(x/a)^2) does not always result in opposite sign** **Status:** open **Group:** None **Created:** Mon Jan 23, 2017 08:30 PM UTC by Torben **Last Updated:** Mon Jan 23, 2017 08:30 PM UTC **Owner:** nobody The resultes below should be of opposite sign. Maxima 5.38.1 http://maxima.sourceforge.net using Lisp CLISP 2.49 (20100707) Distributed under the GNU Public License. See the file COPYING. Dedicated to the memory of William Schelter. The function bug_report() provides bug reporting information. (%i1) integrate(sqrt(1(x/a)^2),x,0,a); Is a positive or negative? pos; %pi a (%o1)  4 (%i2) integrate(sqrt(1(x/a)^2),x,a,0); Is a positive or negative? pos; %pi a (%o2)  4 (%i3) bug_report(); Please report bugs to: http://sourceforge.net/p/maxima/bugs To report a bug, you must have a Sourceforge account. Please include the following information with your bug report:  Maxima version: "5.38.1" Maxima build date: "20161019 00:27:24" Host type: "x86_64suselinuxgnu" Lisp implementation type: "CLISP" Lisp implementation version: "2.49 (20100707) (built on cloud112 [127.0.0.1])"  The above information is also reported by the function 'build_info()'. (%o3) Regards Torben  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: kofoid <eric_kofoid@us...>  20170123 19:47:50

Hi David, Thanks very much for your answer. It helps a lot. Cheers, Eric Eric Kofoid, Ph.D. Senior Project Scientist 318 Briggs Hall Dept. Microbiology University of California at Davis One Shields Avenue Davis, CA 95616 5307526788 eckofoid at ucdavis.edu CONFIDENTIALITY NOTICE: This communication with its contents may contain confidential or legally privileged information. It is solely for the use of the intended recipients. Unauthorized interception, review, use or disclosure is prohibited and may violate applicable laws including the Electronic Communications Privacy Act. If you are not the intended recipient, please contact the sender and destroy all copies of this communication. > On Jan 22, 2017, at 6:36 AM, David Scherfgen <tomasriker@...> wrote: > > This has nothing to do with OS X. Seems like Maxima cannot numerically evaluate gamma_greek. You can work around it by using gamma_incomplete, which Maxima does know how to evaluate numerically. If I'm not mistaken, gamma_greek(10,1) = gamma_incomplete(10,0)  gamma_incomplete(10,1). > > [bugs:#3277] <https://sourceforge.net/p/maxima/bugs/3277/>; gammarelated functions under OS X > > Status: open > Group: None > Created: Thu Jan 19, 2017 12:50 AM UTC by kofoid > Last Updated: Thu Jan 19, 2017 12:50 AM UTC > Owner: nobody > > I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output. > > Sent from sourceforge.net because you indicated interest in https://sourceforge.net/p/maxima/bugs/3277/ <https://sourceforge.net/p/maxima/bugs/3277/>; > To unsubscribe from further messages, please visit https://sourceforge.net/auth/subscriptions/ <https://sourceforge.net/auth/subscriptions/>;  ** [bugs:#3277] gammarelated functions under OS X** **Status:** open **Group:** None **Created:** Thu Jan 19, 2017 12:50 AM UTC by kofoid **Last Updated:** Sun Jan 22, 2017 10:23 PM UTC **Owner:** nobody I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 
From: Robert Dodier <robert_dodier@us...>  20170122 22:23:25

OP, there are at least 4 ways to compute numerical values for `gamma_greek`. ~~~~ (%i1) gamma_greek(10, 1), numer; (%o1) gamma_greek(10, 1) (%i2) gamma(10)  gamma_incomplete(10, 1), numer; (%o2) 0.04043407755671069 (%i3) gamma_incomplete(10, 0)  gamma_incomplete(10, 1), numer; (%o3) 0.04043407755671069 (%i4) gamma(10)*(1  gamma_incomplete_regularized(10, 1)), numer; (%o4) 0.04043407756135764 (%i5) gamma_incomplete_generalized(10, 0, 1), numer; (%o5) 0.04043407755671069 ~~~~ It is a bug that `gamma_greek` doesn't have its own numerical evaluation. It can easily be resolved by punting to one of the methods above, although it's not clear to me which is most accurate. Aside from `gamma_greek`, did you find other gamma functions which do not evaluate numerically? I can't find any examples other than `gamma_greek`.  ** [bugs:#3277] gammarelated functions under OS X** **Status:** open **Group:** None **Created:** Thu Jan 19, 2017 12:50 AM UTC by kofoid **Last Updated:** Sun Jan 22, 2017 02:36 PM UTC **Owner:** nobody I am running maxima 5.38.0 under OS X 10.10.5. Of the various gammaassociated functions, only "gamma" itself produces output. As an example, attempting to run "gamma_greek(10,1)" produces itself as output.  Sent from sourceforge.net because maximabugs@... is subscribed to https://sourceforge.net/p/maxima/bugs/ To unsubscribe from further messages, a project admin can change settings at https://sourceforge.net/p/maxima/admin/bugs/options. Or, if this is a mailing list, you can unsubscribe from the mailing list. 