Soegtrop, Michael | 11 Feb 13:39 2016
Picon

wxMaxima file format: Why is it not plain XML?

Dear wxMaxima Developers,

 

first I am not sure if this is the right forum. is there a separate one for wxMaxima?

 

My question/comment is about the XML format of wxMaxima. I very much appreciate that now the output is embedded into the file, so that I can look at the results without reprocessing it.

 

But I think the format chosen has some disadvantages. It is a XMLish format, but not plain XML, because it has a binary header and binary footer. This has 2 major drawbacks:

 

-          It cannot be processed with standard XML processing tools. E.g. a XSLT processor is convenient to create an overview file for a library or automated documentation. But as is, I first have to strip away the header and footer.

-          I cannot edit the file in a normal text editor. When I do, I get errors on loading the file. Editing the file in another editor makes sense, e.g. for spell checking. And it also would be convenient to use standard XML tools to take a wxMaxima file apart into two parts, one I want to spell check, one I won’t, then correct the part I want to correct and then to put them together again.

-          I cannot automatically create wxMaxima files. I can imagine applications, where the document format and plotting capabilities of wxMaxima  would be convenient e.g. to automatically create and process a measurement report.

 

In short, as the format is, I cannot use any standard XML tools to handle the files, but to a certain extent it is the point of XML to be able to do this.

 

How about ASCII encoding the PNG images? This would result in standard conformant XML files and blow up the file size only marginally (20%).

 

Best regards,

 

Michael

Intel Deutschland GmbH
Registered Address: Am Campeon 10-12, 85579 Neubiberg, Germany
Tel: +49 89 99 8853-0, www.intel.de
Managing Directors: Christin Eisenschmid, Christian Lamprechter
Chairperson of the Supervisory Board: Nicole Lau
Registered Office: Munich
Commercial Register: Amtsgericht Muenchen HRB 186928

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Soegtrop, Michael | 11 Feb 10:21 2016
Picon

Why is pdf_binomial called a Probability Density Function?

Dear Maxima Users,

 

I wonder why pdf_binomial from the distrib package is called a Probability DENSITY Function.

 

I would think a probability density function is something I need to integrate over a variable in order to get a probability. Like when I integrate a normal distribution from x1 to x2 I get the probability that x is between x1 and x2.

 

But the binomial distribution pdf_binomial(k,n,p) returns the probability that in a Bernoulli sequence of n trials and a pass probability of p, I see k passes. I don’t integrate it to get a probability, it already is a probability. So why is it called a PDF then?

 

Best regards

 

Michael

Intel Deutschland GmbH
Registered Address: Am Campeon 10-12, 85579 Neubiberg, Germany
Tel: +49 89 99 8853-0, www.intel.de
Managing Directors: Christin Eisenschmid, Christian Lamprechter
Chairperson of the Supervisory Board: Nicole Lau
Registered Office: Munich
Commercial Register: Amtsgericht Muenchen HRB 186928

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
Jinsong Zhao | 11 Feb 05:16 2016
Picon

No link for ChangeLog in Maxima web page

Hi maintainers,

There are no links for the ChangeLog of recent releases, i.e., from 
5.37.1, in the right pane of Maxima web page. Please add for them.

Best,
Jinsong

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
Gunter Königsmann | 10 Feb 19:54 2016
Picon

Implicit plot resolution

Dear all,

Is there a way to make the circles produced by the following command
have a few more edges?

load("draw");draw2d(
    makelist
    (
        implicit(x^2+y^2=i/10,x,-1,1,y,-1,1),
        i,1,10
    ),
    ip_grid=[100,100],
    ip_grid_in=[100,100]
);

The command seems to produce polygons that are quite good approximations
of perfect circles. But if you zoom in a little (on my computer it is
the + key) it looks like the points the polygon connects look like their
coordinates were somewhat rounded.

Thanks a lot,

and kind regards,

 Gunter.

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
Marco Verpelli | 10 Feb 17:56 2016
Picon

Bernoulli numbers

I was looking at the PSLQ code by Thomas Baruchel, and in his example there is a list of the first 10 nonzero Bernoulli numbers.

 1 => -5.0b-1
 2 =>  1.666666666666666666666666666666666666666666666666666666666666667b-1
 3 => -3.333333333333333333333333333333333333333333333333333333333333333b-2
 4 =>  7.575757575757575757575757575757575757575757575757575757575757576b-2
 5 => -2.531135531135531135531135531135531135531135531135531135531135531b-1
 6 =>  1.166666666666666666666666666666666666666666666666666666666666667b0
 7 => -7.092156862745098039215686274509803921568627450980392156862745098b0
 8 =>  5.497117794486215538847117794486215538847117794486215538847117795b1
 9 => -5.291242424242424242424242424242424242424242424242424242424242424b2
10 =>  6.192123188405797101449275362318840579710144927536231884057971015b3

When I try to reproduce the results I get a different list ?!?!

Maxima version: "5.37.3"
Maxima build date: "2016-01-11 20:53:44"
Host type: "i686-pc-mingw32"
Lisp implementation type: "GNU Common Lisp (GCL)"
Lisp implementation version: "GCL 2.6.11"

 1 => -0.5
 2 =>  0.1666666666666667
 3 => -0.03.333333333333333
 4 =>  0.02.380952380952381
 5 => -0.03.333333333333333
 6 =>  0.07575757575757576
 7 => -0.2531135531135531
 8 =>  1.166666666666667
 9 => -7.092156862745098
10 => 54.97117794486216

The Maxima 5.37.3 results agree with Wolfram Alpha.

Best regards

Marco Verpelli
------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
David Scherfgen | 10 Feb 13:02 2016
Gravatar

Antiderivative of sin(x^3) is wrong by a factor of -1/2

Hello,

(I will open a bug report, but I'm also sending this to the mailing list)

Maxima gives an antiderivative of sin(x^3) that's wrong by a factor of -1/2:

(%i1) integrate(sin(x^3),x);
(%o1) -(gamma_incomplete(1/3,%i*x^3)+gamma_incomplete(1/3,-%i*x^3))/12

(%i2) diff(%,x);
(%o2) -(3*%i*%e^-(%i*x^3)-3*%i*%e^(%i*x^3))/12

(%i3) demoivre(%);
(%o3) -(3*%i*(cos(x^3)-%i*sin(x^3))-3*%i*(%i*sin(x^3)+cos(x^3)))/12

(%i4) trigsimp(%);
(%o4) -sin(x^3)/2

Any ideas how this can happen?

I think this is quite serious. Even much older Maxima versions give
this wrong answer.

The correct answer would be:
(gamma_incomplete(1/3,%i*x^3)+gamma_incomplete(1/3,-%i*x^3))/6

Best regards,
David Scherfgen

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
Roland Salz | 9 Feb 23:25 2016
Picon

Testing the index notation

Hi Robert,

 

annexed you find my rudimentary test script frame for the index notation on the console.

 

Having integrated verify_display_2d and verify_display_1d as you suggested, I got a failed test with the first one where it should not fail. From your examples in rtest9a.mac I noticed that you regularly inserted one blank after each fraction bar. Then I found out that I had to insert two blanks after the ‘foo’ in order to make the test pass.

 

What is the underlying principle in having to append such blanks at the end of certain lines? That is: how do I know where and how many of them to append?

 

Best regards,

Roland

 

Attachment (rtest_index_notation_console.mac): application/octet-stream, 1723 bytes
------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
有坂富雄 | 9 Feb 21:28 2016

Draw + VTK

Hi Mario,

I found some bugs in 'vtk.lisp'.
Function vtk2d-explicit and vtk2d-parametric do not return values.
So, my patch is as follows.

Thanks,

Tomio Arisaka


$ diff -u /sw/share/maxima/5.37post/share/draw/vtk.lisp ./vtk.lisp
--- /sw/share/maxima/5.37post/share/draw/vtk.lisp    2016-02-08 18:42:07.000000000 +0900
+++ ./vtk.lisp    2016-02-10 03:49:41.000000000 +0900
<at> <at> -1749,8 +1749,8 <at> <at>
     (when (not (subsetp (append (rest ($listofvars xfun)) (rest ($listofvars yfun))) (list par)))
        (merror "draw2d (parametric): non defined variable"))
     (setq *plot-realpart* (get-option '$draw_realpart))
-    (setq f1 (coerce-float-fun xfun `((mlist), par)))
-    (setq f2 (coerce-float-fun yfun `((mlist), par)))
+    (setq f1 (coerce-float-fun xfun `((mlist) ,par)))
+    (setq f2 (coerce-float-fun yfun `((mlist) ,par)))
     (setf result
        (loop
           do (setf xx ($float (funcall f1 tt)))
<at> <at> -1789,7 +1789,8 <at> <at>
       (case linetype ; translate some gnuplot codes into vtk codes
         (0 3)
         (6 4)
-        (otherwise linetype) ) ) ))
+        (otherwise linetype) ) )
+    str ))
 
 
 
<at> <at> -2295,7 +2296,8 <at> <at>
       (case linetype ; translate some gnuplot codes into vtk codes
         (0 3)
         (6 4)
-        (otherwise linetype) ) ) ))
+        (otherwise linetype) ) )
+    str ))
 
 
 

On 2016-2-8 18:50, Mario Rodriguez <biomates <at> telefonica.net> wrote:



El 08/02/16 a las 01:41, Dimiter Prodanov escribió:


The following example fails in  

wxMaxima version: 15.8.2
Maxima version: 5.37.2
Maxima build date: 2015-09-18 21:03:17
Host type: i686-pc-mingw32
System type: Win32 6.1.7601 (Service Pack 1) X86
Lisp implementation type: SBCL
Lisp implementation version: 1.2.7

load("draw") $ draw_renderer : 'vtk $ draw3d( color = blue, triangle([0,0,0], [2,0,0], [3,1,2]), color = yellow, triangle([0,0,0], [2,2,0], [3,1,2]), color = gray, triangle([2,0,0], [2,2,0], [3,1,2]) ) $ I can only make it work in gnuplot.best regards,Dimiter


Hi Dimiter,

Version 5.37.2 does not contain the latest changes in draw. In order to test VTK you should:

a) Download the current draw version from the repository. In particular these two files:
    http://sourceforge.net/p/maxima/code/ci/master/tree/share/draw/vtk.lisp
    http://sourceforge.net/p/maxima/code/ci/master/tree/share/draw/grcommon.lisp
b) Place them wherever you want in your system and write:
    load(draw) $
    load(path_to_new/grcommon.lisp) $
    load(path_to_new/vtk.lisp) $
c) Make sure that your system has VTK6 installed with TCL support. (It is already bundled in the windows installer of Maxima 5.37.3.) This new version of draw is not compatible with the VTK5 series. Version VTK7 has been released short time ago, but I did not test it yet.

Of course, these steps will not be necessary once Maxima 5.38 has been released.

If you follow the steps described above and get into problems or find bugs, please let me know about them.

Thanks for your bug report.

--
Mario
------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
_______________________________________________
Maxima-discuss mailing list
Maxima-discuss <at> lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/maxima-discuss
sholman | 9 Feb 19:19 2016
Picon
Picon

is there a way to continue demo("") after a quick calculation in terminal window?

Hi, I love the demo("filename.mac")$ utility but was wondering if there 
was a way to continue doing the demo (hitting RETURN) after I do a quick 
calculation in the middle.

For instance, if I am hitting RETURN after demo("filename.mac"), and 
then I come across a line I want to check a calculation with, I'll type 
something at the _ prompt, get a result, and then want to continue - 
just like after a break(), but demo() won't continue.

Is there something simple I'm overlooking?

Thank you all!

Brett

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
David Scherfgen | 9 Feb 15:33 2016
Gravatar

Wrong implementation of "erfi"?

Hi,

according to my knowledge and the Maxima manual:

erfi(z) = -%i*erf(%i*z)

But plotting both expressions, I get a difference. The curves are
equal in the approximate interval [-0.476, 0.476], but they have
different sign outside that interval.

erfi suddenly switches signs at -0.476 and 0.476 ...

Here's the plot:
https://dl.dropboxusercontent.com/u/845597/erfi.png

This is current Git version of Maxima. The older one at
maxima-online.org doesn't have this bug.

Best regards,
David Scherfgen

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140
Andrei Zorine | 9 Feb 11:22 2016
Picon

how to plot a complex function?

Dear group!

Suppose two *complex* quantities, z1 and z2, are in relation F(z1, z2)=0. 
Further, suppose the equation can not be solved explicitly for either z1, or 
z2. Can Maxima plot z1 as a function of z2 (some complex branch possibly) in 
implicit way (e.g. a 3d plot with z-axis for abs(z1) and hue from carg(z1))?

Are there general algorithms for this task, not necesserely impliemented in 
Maxima? I haven't found any serching in google.

Thanks in advance,
Andrei Zorine

------------------------------------------------------------------------------
Site24x7 APM Insight: Get Deep Visibility into Application Performance
APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
Monitor end-to-end web transactions and take corrective actions now
Troubleshoot faster and improve end-user experience. Signup Now!
http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140

Gmane