I have a code in Maple that solves a system of polynomial equations (18 equations, 18 variables) using the command 'RootFinding[Isolate]' (it takes to Maple about 40 minutes to find all the solutions to the equations). Is it possible to see the implementation of the Isolate command? if not, is there any other open source code that implements the equivalent to Isolate command and also able to solve such equations? (I'm talking about the  numeric solutions. not symbolic).

I already tried to use the command 'solve_poly_system' in python (sympy package) which is open source. but the function didn't return a solution after a very long time (and I used a system of 9 equations and 9 variables which were much simpler than the equations Maple was able to solve using Isolate command).

I'm attaching the equations Maple solved using Isolate command. and also the equivalent in python with the simpler equations (that it wasn't able to solve).

solving_18_polynomial_equations.mw

trying_to_solve_polynomial_equations_python.txt

I'd appreciate any idea of how to use open source code that produces the same results as the code in Maple (if it's possible).

Thanks

Hey MaplePrimes

I have a problem with opening some of my documents in maple. When i try to open them, a box pops up with the text "How do you want to open this file?" with the options "Maple Text, Plain Text, Maple Inputs". None of these options work (document is blank) and "Plain Text" and "Maple Inputs" causes maple to crash. Is it possible to recover the documents?

ps. when i try to upload one the documents this happens:

 
Maple Worksheet - Error
Failed to load the worksheet /maplenet/convert/Maple_docs_til_mat.mw .

Download Maple_docs_til_mat.mw

Maple_docs_til_mat.mw
 

What is the meaning of the output of the command Weierstrassform?

E.g. what does mean x0^3 - 7*x0 + 88 +y0^2 ?

Maybe y0^2 = x0^3 + 7*x0 -88?

Hey everyone
 

I am trying to create a procedure, which includes a part, where a function is plotted. I want to be able to declare the range of the plot within the procedure parameters as x=range - example when executing: procedureName(function,x=range). I have isolated a small part of the procedure and the different approaches I have tried:

For simplicity's sake let's call the procedure test.

test:=proc(function,{x::range:=0..1})

#I would then like to be able to plot the following:

plot(function,'x'=x);

end proc:
 

However, the single quotes do not prevent x from being evaluated within the procedure, and there interprets the second argument of 'plot' as 0..1=0..1. Another approach I tried was the following:

test:=proc(function,{x::range:=0..1})

local xRange:=x;

unassign('x');

plot(function,x=xRange);

end proc:

However, I run into the same problem as above - the elements within the single quotes are still evaluted within the procedure. A workaround would be to simply rename the range parameter, for example:

test:=proc(function,{xInr::range:=0..1})

plot(function,x=xInr);

end proc:

I would really like to be able use x as the parameter, though.

Thanks in advance!

Dear Maple users,

I have a big exporession in the form f(s)/g(s) where f and g are polynomial involving irrational coefficients. Both f and g consist of hundreds of terms. One way to get the inverse laplace is to factor the expression and then take the invlaplace.

The code is as follows:

with(inttrans):

p:=factor(f(s)/g(s));

sol:=invlaplace(p,s,t);

I get sol as a function of time, which gives satisfactory result for t=0 but  grows bigger and bigger with increasing value of time. From the limit theorum I see that  Limit (s*f(s)/g(s)) as s--> 0 is finite, which means that the invlaplace is not correct at very large value of time.

What is the error in the process?

Thanks,

Hi all,

I attached a program here and the desire is the calculation of Nu.As you could see through the attached file, the F2(r) function contains an 'integral' which makes it difficult if I want to calculate F2(r=1) and it goes the same for D(F2)(r=\phi). So I firstly corrected the F2(r) function as FF2(r) in which the parameter R is calculated in the top of the file. Afte I replace the values of \phi and 1 in the FF2(r) not in F2(r) just because of the integral. The Nu is calculated finally, but it differs from the initial guess. How could I make a loop in order to correct the initial guess by replacing the first Nu calculated in the end of the program??

Ther is a problem and that is the "k" which must be entered in this part in each loop:

Digits:= 10:

K:=Nu->( k )

I do not know whether or not it is possible to change it so that the program identify the 'k' and replace it in the "K:=Nu->( k )" for k.

SolveCode.mw
 

Hi,

I am collecting first order terms in psi and phi in the following expression.

Basically I want to keep coefficients for which (power of phi + power of psi) < 2 and ignore the rest.

Here, H, a, psi, phi are all functions, not variabels.

I tried following answer but it doesn't work due to derivatives.

https://www.mapleprimes.com/questions/37228-Ignore-Terms-Above-Particular-Orderpower

(Error, selecting function must return true or false)

Thanks a lot,

Rahul

Dear all

How can I finish the attached code and I plot   the set of points (x[i],y[j], U[i,j]) 

U[i,j] represent the value of U at the position  (x[i],y[j])

Many thanks for your help

Plotsetofpoint.mw

Before Maple 2015, units were enlosed in [[ double brackets ]] when displaying 2D input and output. In later versions these brackets are displayed only when editing input, and never in output, which in my opinion detracts from readability. Does anybody know if it is possible to revert to the old behaviour?

Here is my code:

with(DEtools): with(plots): with(linearAlgebra):

DE1 := diff(y(x), x) = d^3*y(x)/dx^3+3*d^2*y(x)/dx^2+4*d*y(x)/dx+12*y(x);

DEplot(DE1, [y(x), x = 0 .. 5][[y(0) = -3, (D(y))(0) = 0, ((D^(2)(y))(0) = 0, ((D^(3)(y))(0)]]);

I get the Error message: Error, (in DEtools/DEplot) called with too few arguments

Any help would be greatly appreciated, thanks.

Many of you enjoyed our profile on one of our developers, Paulina Chin, so we’re happy to bring you another one!

Today, we’ll be talking with John May, Senior Developer of Maple. Let’s get started.

1. What do you do at Maplesoft?
   Until recently I was consulting on-site at the NASA Jet Propulsion Laboratory helping people there more effectively solve their engineering problems using Maplesoft products.  But my main job that I am back to full time now is the development and maintenance of various parts of the Maple library.
    
2. What did you study in school?
   I studied both Pure and Applied Mathematics at the University of Oregon,  focusing a lot on Abstract Algebra.  In graduate school, I specialized more in computation mathematics like computer algebra and numerical analysis.  My Ph.D. work focused on effective numerical algorithms for problems in polynomial algebra – with implementations in Maple!
    
3. What area(s) of Maple are you currently focusing on in your development?
   Right now I am focused on addressing complaints I’ve gotten from engineers about the usability of units with other parts of the math library.
    
4. What’s the coolest feature of Maple that you’ve had a hand in developing?
   A lot of the cool things I’ve built live pretty deep in the internals of Maple.  I’ve done a lot of meta-heuristic tuning to seamlessly integrate high-performance libraries into top-level Maple commands.
   
   I had a lot of fun developing a lot of the stuff for manipulation and visualization of colors in the ColorTools package.
    
5. What do you like most about working at Maplesoft? How long have you worked here?
   I started working at Maple in 2007, but I’ve been a Maple user since 1997.  I love being part of the magic that brings powerful algorithmic mathematics to everyone.  The R&D team is also full of eccentric nerds who are great fun to work with.
    
6. Favourite hobby?
   It varies by the season, but right now it is prime for mountain biking in southern California.  I ride my local trails a couple times a week, and when I get I chance, I love to get away on epic bikepacking adventures (like this one: https://www.bikemag.com/features/two-wheeled-escape-one-hour-from-l-a/  this is me: https://cdn.bikemag.com/uploads/2016/05/16File.jpg ).
    
7. What do you like on your pizza?
   Anything and everything. Something different every time. My all-time favorite pie my from grad school days is the “Rio Rancho” from the dearly departed That’s Amore Pizza (which was next to the comic book store and across the street from North Carolina State University).  It was an olive oil and mozzarella pizza with chopped bacon that was covered in sliced fresh roma tomatoes and drizzled with ranch dressing when it came out of the oven. 
    
8. What’s your favourite movie?
   It’s really hard to pick just one.  So, I’ll go with the safe answer and say the greatest movie of all time, and “Weird Al” Yankovic’s only foray into movies, UHF, is my favorite.
   http://www.imdb.com/title/tt0098546/
    
9. What skill would you love to learn? (That you haven’t already) Why?
   Another hard one.  I feel like I’ve dabbled in lots of things that I would like to get better at.  At the top of the list is probably unicycling.  I’d love to get good enough to play Unicyle Football or do Muni (mountain unicyling).
   https://en.wikipedia.org/wiki/Mountain_unicycling
   http://www.unicyclefootball.com/
    
10. Who’s your favourite mathematician?
    Batman. https://youtu.be/AcMEckOyoaM

Hi I have a function i d live to plot and integrate but maple tells me there is a probleme with the range when i want to plot and will not give me a numerical value of the integral :

What am i doing wrong ??

I am trying to repeat the 'first example' in Stephani & MacCullum, Differential equations see chapter 16 and (16.5).  The differential equation is,

pde := b1(x, y)*(diff(u(x, y), x))+b2(x, y)*(diff(u(x, y), y)) = 0

and I would like to compute DeterminingPDE.  The textbook answer is linear in b1 and b2 (as my 'by hand calculation' is) but Maple's,

DeterminingPDE(pde)

contains quadratic terms such as b1^2, b1*b2 etc.  I don't understand the appearance of the quadratic terms.  Is it possible for Maple to return an answer which is linear in b1 and b2?

In the link below to my worksheet I have attempted to construct the Fourier series for the function 1-x.  I should mention that the I often interchnge the variables x & t (probably a bad habit).  It appears that I am getting something close, but my term for n=0 appears to be undefined.  Given the plot I generated the term for n=0 should be approximately 1/4 to get the base of the triangle to intersect the t or x axis.

So I thought if I can get another set of eyes on my work maybe someone can point out what I have incorrect.

Appreciate any assistance

asymmetric_triangle_fourier_series.mw

I'm learning Maple with the provided user manual, and it clearly says to use ctrl+= to display inline results. It simply doesn't work on  my system. Running Maple 2017.3. Is this a bug?