Hello, 

I've created a model of a rose flower that I'd like to export as a high quality picture without the white space.  Is there a way to do this?  

deq := FunctionAdvisor(DE, hypergeom([a, b], [c], z))[2, 1];
          diff(f(z), z, z) = ((-a-b-1)*z+c)*(diff(f(z), z))/(z*(z-1))-f(z)*a*b/(z*(z-1))

dsolve(deq, f(z));
Error, (in dsolve) Not an ODE w.r.t. f(z)

Apparently FunctionAdvisor is doing something special for the variable names in order to avoid collisions with user variables, but if I want to use the output programmatically, without copying and pasting it, is there a better way than something like this:

dsolve(subs(op(0, indets(deq)[-1]) = f, deq), f(z));
          f(z) = _C1*hypergeom([a, b], [c], z)+_C2*z^(1-c)*hypergeom([b-c+1, a-c+1], [2-c], z)

Hi there,

I was using maple 2016 and would like to plot the following:
Question.mw

However,when I run the code,only the axis are shown.How should I fix it?

Thanks!

For a numerical integration in Maple, one only adds evalf(...) to the integration commend. Now my question is that is it also the case for the multi-integration? Also consider that the multi-integration of the Maple doen't give any answer for some examples. So I think putting evalf around them, doesn't help. How can someone compute a numerical approximation of a multi-integration by Maple? Is there any prepared commend or package for it?

i have a PDE of heat equation, i can solve this using two boundary and initial condition, but i can not handle it using second boundary condition and maple also can not ! hoc can i solve it ?
 

Download pde.mw

Dear 

Hope everyone is good. I am face to attaine the converges solution of the attached problem. Please have a look and fix my problem. I am waiting your response

diverges.mw

I have 6 equations, eq1,...,eq6, with 6 unknowns a,b,c,d,f,g. I can find all the solutions with "solve".
If "solve" has found only one solution, then the ith solution number is given by sol[i]. For example,
c:=sol[3]. This comes from

sol:=solve({eq1,eq2,eq3,eq4,eq5,eq6},{a,b,c,d,f,g});

However, if solve has found more than one solution,
then the jth element of the ith solution is sol[i][j]. It is thus imperative to know how many solutions
the command "solve" has found. Also, how do we restrict "solve" to look only for real soultions?
Thank you very much!

mapleatha

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