Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I would like to be able to use python's matplotlib from Maple. I am doing the following currently

with(Python)
Start('home' = "/opt/anaconda3/lib/python3.9", 'path' = "/opt/anaconda3/lib/python3.9/site-packages")
ImportModule("matplotlib as mpl")

and the error I get back is

Error, (in Python:-ImportModule) ModuleNotFoundError: No module named 'matplotlib'
['Traceback (most recent call last):\n', '  File "<string>", line 1, in <module>\n', "ModuleNotFoundError: No module named 'matplotlib'\n"]

The string I passed into path is the directory where matplotlib seems to be located as far as I can tell. I looked this up in my IDE Spyder, but using the command print(matplotlib__file__).

UPDATE: but actually if I open a new worksheet and execute 

with(Python)
Start('home' = "/opt/anaconda3/lib/python3.9", 'path' = "/opt/anaconda3/lib/python3.9/site-packages")

The error is that 

Error, (in Python:-Start) unable to start Python: could not create process, mpython

So Maple isn't finding Python apparently.

I tried for 'home' also "/opt/anaconda3/bin", did not work.
 

Hello guys.
Is there any way to frame the styles: Input 2D, Output 2D and Text?

Oliveira.

A robust largescale index-1 DAE (and stiff ODE) solver has been developed in Maple.
Details about our approach are provided in a paper just submitted. 
See arxiv at
https://arxiv.org/abs/2212.02630

We encourage everyone to test these codes and report bugs. All the examples can be run with Intel's Pardiso (provided users have access to libraries) by calling DAESolverP.txt and by calling "IMPDAEP" instead of "IMPDAE". This is useful for large-scale problems. The symbolic capability and ListTools search capability of Maple are very good and can be used for developing optimization solvers as well.

I would like to know if CPU time/memory usage can be reduced significantly. In particular,  for examples 5 and 6. Some ways to contribute include

(1) Running the code in evalhf or compiled form. This may be hard.
(2) Providing options to run other parallel open-source linear solvers (eg., MUMPS).

(3) Other examples that show the use of the developed solver. We are able to solve > 100,000 DAEs.

(4) Helping in converting the code to Maple 14 or earlier (by doing sparse LU Decomposition. Just using LinearSolve will slow down the code).

 

Please avoid ~,*, etc (shortcuts) unless it improves the speed of calculation.
Thanks 
Dr. Venkat Subramanian

Download examples_1_3.mws

Download example_4.mws

Download example_5.mws

Download example_6.mws
 

DAESolver.txt

DAESolverP.txt

Trying to use the following question dialog

with(Maplets[Elements]);
maplet := Maplet(QuestionDialog("Is x > 0?"));
result := parse(Maplets[Display](maplet));

Displays question with Yes and No buttons but whichever button is clicked the result is always false.  This was based closely on an example from the documentation.

So, how can I tell which button was clicked.  Also, the close button could have been clicked - how should that be dealt with?

Hello guys.

Is there any way to convert a worksheet to a document and vice versa?

Regards,

Oliveira.

xwhat is the fenchel conjugate transform of f(x)=x*x2*x3 over xi in R?  the general question is what is the fenchel transform of a polynomial that is not a quadratic?

Fyi, this looks like a bug.

But I do not know if it is new or not. I am using Maple 2022.2 on windows 10

Why it fails first time?

(ps. unable to paste content of worksheet here. looks like forum is not working, will paste an image)

Here is the code also in plain text

restart;
delMdelY:=-2*y^2/((x^2 - y^2)^2*(x - k*sqrt(x^2 + y^2))) - 1/((x^2 - y^2)*(x - k*sqrt(x^2 + y^2))) - y^2*k/((x^2 - y^2)*(x - k*sqrt(x^2 + y^2))^2*sqrt(x^2 + y^2));
delNdelX:=-2*x/(x^2 - y^2)^2;
if simplify(delMdelY - delNdelX)=0 then
    print("ok");
else
    print("not ok");
fi;

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/simplify_fails_in_proc.mw .

Download simplify_fails_in_proc.mw

please, help me..this is my model

I have a fairly simple worksheet with a 56x2 matrix M1 which I am trying to export to Excel with the following 

ExcelTools:-Export( M1, "C:/Temp/M1.xlsx", 1, "B2")

This causes the following message to appear

After that the program is frozen and totally unresponsive.  The only way to proceed is end the task in Task Manager..

I am trying to upload the worksheet but am not sure if I have succeeded

Download Calculation_of_pH_from_known_values_of_Vb.mwCalculation_of_pH_from_known_values_of_Vb.mw

I created a Github repository with files to reproduce the example that my question is based on.

I am using Grid:-Run to compute a 10000x3 matrix in each of 100 parallel runs. I then average the 100 matrices to obtain a final 10000x3 matrix, which I saved to an .m file.

It contains a variable called with the matrix.

Each row of the matrix is a 3d point, and the entire matrix represents a trajectory on a sphere. I'd like to display an animation of this trajectory. Now, my question is not about how to do this, but rather about the memory requirements necessary to do so.

I can plot trajectories with, say, 1000 points, but apparently 10000 requires too much memory.

The Github repository I linked to contains an example of all of this.

The worksheet matrixTest.mw reads the matrix from the .m file and reads some procedures defined in the .mpl file that are used for creating the animation data.

Then, a call to the procedure animateMatrices basically calls plots:-display to display the matrix as an animation.

I've been able to display the first 5000 rows of the matrix.

What type of memory is running out when I try to display the animation of the full 10000 rows?

Here is a screenshot of my Activity Monitor right before I get the

Error, out of memory error while processing result

Hi! I need to simplify a polynomial over an arbitrary field F on a radical variable ν, such that for some power j, ν is in F.

the polynomial takes the form:

k0+a1ν1/j+...+anνn/j

but, since νj is in F, this divides every element νi/j into a cotient group of order j, so this can be rewritten as a radical extension of F, in the form: 

k1+(a1+...)ν1/j+...+(ak+...)ν(j-1)/j

where ki, ai are in F. I feel like this is a very straightforward technique for handwritten algebra, but i can't see a command for this on Maple 18. Nor the simplify(..., radicals) Or the combine(...,radicals) seems to help here. Maybe there's an special command for this? Must i do it myself? Please, any help is aprecciated.

This question is based on the following worksheet I created specifically to show my question: Animation_Test.mw

Unfortunately, MaplePrimes gives me an error when I try to insert the contents here directly.

Basically my question is about the syntax of plots:-display.

I have a sphere

sphereDisplay := plottools:-sphere([0, 0, 0], 1, transparency = 0.9);

and something which is the result of calling plots:display on a sequence of two other other plots:-display.

arrowAndTrailDisplay := plots:-display(seq(plots:-display(computeArrowPlotFrame(m, t, red), computeTrailPlotFrame(m, t, windowSize, red)), t = times), insequence = true)

You can see in the attached worksheet what computeArrowPlotFrame and computeTrailPlotFrame are, but they are basically lists of plot structures. arrowAndTrailDisplay is thus an animation containing two different sequences of frames.

I wish to have Maple plot the sphere and the animation together. This works as expected (a static sphere and the animation of an arrow and a trail of the path of the tip of the arrow)

plots:-display([sphereDisplay, arrowAndTrailDisplay]):

but the following does not work as expected (instead it plots each frame separately)

plots:-display([arrowAndTrailDisplay])

adding insequence=true to the latter makes it work. But my question is why do I have to add that? I can't understand, from syntax standpoint why the former plots arrowAndTrailDisplay as an animation but the latter does not.

Hello everyone,

I am trying to solve a system of six equations through a matrix and a vector. The matrix is 6x6, so the function "LinearSolve" should find a solution 

The matrix is

Matrix(6, 6, [[1., -1., 1., -1., 1., -1.], [1., 1., 1., 1., 1., 1.], [-2., 1.618033989, 30.94427190, -153.8246851, 371.1559479, -572.9674774], [-2., 0.6180339876, 22.94427191, -41.15905356, -37.04759741, 149.3606798], [-2., -0.6180339876, 13.05572810, 22.82468509, -20.15594802, -81.03252254], [-2., -1.618033989, 5.055728096, 28.15905368, 68.04759752, 104.6393203]])

and the vector associated with the system is

Vector[column](6, [-1/3, -1, 0, 0, 0, 0])

I am trying to solve this system with

a := LinearSolve(M3, v);

but it stays evaluating.

I have solved a similar system (5x5) with this function (LinearSolve) and it took less than a second, so I dont understand why it takes so long in this case.

Thank you in advance.

Consider the following code in an .mpl file

$define MACRO1 "1"

$define MACRO2 "2"

 

$ifdef __FILE__

$define MACRO3 "3"

$endif

 

$ifdef MACRO1

$define MACRO4 "4"

$endif

 

f := proc()

print(MACRO1, MACRO2, MACRO3, MACRO4, __FILE__):

end:

If we run this from a worksheet and call the procedure f, the printed output is

"1", "2", MACRO3, "4", "/Users/zenterix/..."

Why does $ifdef __FILE__ evaluate to false? Ie, why doesn't MACRO3 get defined? I thought __FILE__ was a predefined macro.

I formatted all the input math expressions with a different style and then exported them to create a style set called Mystyle1. However, when I load Mystyle1 and apply it to a document, the functions that were typed partially obey the style, but the expressions that were inserted through the palette do not; font and size are different and bold is not applied.
Does anyone know why this happens?

If anyone wants to test it, Mystyle1 and a simple document are attached.

Regards,

Oliveira

Example.mw

Mystyle1.mw

1 2 3 4 5 6 7 Last Page 1 of 1967