Maple 2020 Questions and Posts

These are Posts and Questions associated with the product, Maple 2020

Can I specify a unit for a symbol using the Units package, so that when I substitute a value for that symbol into an expression, the units for the value are already known?


Problem description

1) Define a function in terms of variables and parameters,

2) FInd the optimal solution

3) Plot the optimal solution in terms of a single parameter



What I have done is very basic since I haven't use Maple for many years. Here is an example

1) f(x;a,b,c):= ax^2+bx+c; (I don't knpw how to tell Maple that {a,b,c} is a set of parameters)

2) solve(f(x)=0,x); (I could use something like solve(f(x)=0,x,'parametric','real') but I am not interested in so detailed solution)

3) Here is my main problem. I want to save the first solution of x in optx(a;b,c) and the second in opttx(a;b,c) but I don't know how to do it (again, here a is my variable and b,c my parameters)

4) I also don't know hot to plot optx(a) as a increases, whereas the values of {b,c}={e.g., 2,3}


I would appreciate any resources/guidance 


I purchased the license of Maple 2019 and I am using the same version up to now. I am wondering if I download the latest version, can I use the former license? When I click on "check for updates" there is nothing for updating but when I want to run a Maple file written in the latest version, I face an error. 

Hello colleagues

I use Maple 2020 and I try to create Histogram of Poisson with a mean 200 [ Histogram(Sample(Poisson(200),10000))]. but I get an incorrect result in the form of a histogram as uniform distribution. But In my Maple 2015 I get a correct result.what is the problem?


I use Syrup to solve for the transfer function of a circuit, and in particular, the voltage at a given node.   If I use the "trans" option, the solution generates the expected result, in the expected format.  Syrup solves the differential equations and the resulting time domain equation is as expected.

If I then use the 'AC'  option, I get the result in the s-domain.  I then solve for the step response, and take the inverse Laplace transform.  The solution is given in terms of sinh and cosh.   

If I plot the voltage in either solution as a function of time, the plots are identical.  So, the sinh/cosh solution must be somehow simplifiable so that it "looks" like the time domain solution I get from the transient analysis.





interface('displayprecision' = 8)



assum := R1::real, C1::real, R2::real, C2::real, R1 > 0, R2 > 0, C1 > 0, C2 > 0, t::realvin::real, vin > 0

vin::real, 0 < vin


ckt := [vin(5), Rm(1.0*10^7), Cm(1.8*10^(-11), ic = 0), `&//`(R1(9.*10^6), C1(1.5*10^(-11), ic = 0)), `&//`(R2(1.*10^6), C2(1.5*10^(-11), ic = 0))]

[vin(5), Rm(10000000.0), Cm(0.1800000000e-10, ic = 0), `&//`(R1(9000000.), C1(0.1500000000e-10, ic = 0)), `&//`(R2(1000000.), C2(0.1500000000e-10, ic = 0))]




ckt1 := [vin(0), Rm(0.10e8), Cm(0.18e-10, ic = 2.5), `&//`(R1(0.9e7), C1(0.15e-10, ic = 2.25)), `&//`(R2(0.1e7), C2(0.15e-10, ic = .25))]

[vin(0), Rm(0.10e8), Cm(0.18e-10, ic = 2.5), `&//`(R1(0.9e7), C1(0.15e-10, ic = 2.25)), `&//`(R2(0.1e7), C2(0.15e-10, ic = .25))]



Solution doing Transient Analysis

deqs, rest := Solve(ckt, 'tran', 'returnall')

sol1 := dsolve(deqs):

Reassign rest to include the solution.

rest := eval(rest,sol1):

Compute the voltage across the capacitor Cm.

tmp1:=simplify(eval(v[Cm](t), sol1));







p10 := plot(tmp1, t = 0 .. 0.00050000, color = 'blue', linestyle = 'dash')

p10 := plot(tmp1, t = 0 .. 0.5e-3, color = 'blue', linestyle = 'dash')


Solution from AC analysis:

sol, rest := Solve(ckt, 'returnall')

vnode2 := eval(v[2], sol)







Impulse Response

temp1 := invlaplace(vnode2, s, t)






Step Response

temp2 := invlaplace(vnode2/s, s, t)



p1 := plot({temp2, (1-exp(-1))*2.5}, t = 0 .. 0.5e-3)

temp3 := solve(temp2 = (1-exp(-1))*2.5, t)



p2 := plot([temp3, temp3], [0, 2.5])

display({p1, p10, p2})



What is going on?  Where are the sinh and cosh terms in the solution coming from and why wont it simplify so that it looks like the other solution if they are identical?  Is there something special I need to do when using the invlaplace function?





This worksheet attempts to implement figures in Section 3 of the book "Indra's pearls" by Mumford, Series and Wright using the  command complexplot.

The worksheet seems to display Figure 3.1 on page 63 accurately but fails to correctly display Figure 3.4 on page 67.

Can complexplot be coded to accurately display this figure?

Initially I have a system of w(x,y,t) and f(x,y,t). I have a set of two partial differential equations with this system. Later on, the system is modified into w(x,y,t,z) andf(x,y,t,z) where x and z are related due to the presence of boundary layers. Because of this, the partial differential expressions with respect to x has to change.

For example, initially: diff(w(x,y,t),x)=diff(w(x,y,t),x)
Later on, when the system gets modified in the presence of boundary layers on application of chain rule: diff(w(x,y,t,z),x)=diff(w(x,y,t,z),x)+diff(w(x,y,t,z),z)*diff(z,x)

This change has to happen for all higher order derivatives of x. In my case, I have upto 4th order derivative in x. How do I achieve these changes in the expressions in Maple in an elegant manner?

I want to create a random polynomial using randpoly() with coefficients in a given list A := [1,1,2,5,6].

For that, I need to create a function that works as rand(1..9) but that when called returns a random element of A instead of a random number between 1 and 9.

What I have tried so far is

r := rand(1..numelems(A))

However, that returns a random element of A, not a function that returns a random element of A.

I've also tried


But that gives an error.

One way of solving this problem would be to find out how to create a `..` type with the desired elements, something that works as 1..9 but with the elements of A. But convert(A,`..`) doesn't work and the rest of the things I tried don't work either.


OK, so my real-life example is I have a 3m gap between the house and the garage which I want to hang a sunshade. 
The hanging point on the house is 2.6m high and on the garage 2.1m (above ground). 

I have a choice between 2 sun sails, and I want to know -approximately- what sag I can expect from each.

I realise this is not exactly the same as the hanging wire problem i posted before.. but it my attempt anyway. Can anyone improve?

there is this ref


Hi. following on from

Second case: unequal poles.

I tried to work with vv's solution, but I got a problem... the required formula is

y=10.85378553130*cosh(0.0921337534371039*x) - 10.85378553130


I am trying to follow this paper

Lets start with the easier problem, equal poles. Assume that the length of the cable is 120m and the two poles have equal height of 50m. Our goal is to determine the minimum distance between the two poles that will prevent the cable from touching the ground.

I was trying to get Maple to agree with their derived formula, namely

y(x)= 11*cosh(1/11*x) - 61,

but I think I have not set the IC's correctly. or provided for the length of the cable.


//can't solve the above directly, or maybe someone clever here can

//a is a constant
DIV := diff(y(x), x, x) = a*sqrt(1 + diff(y(x), x)^2);
RV := y(0) = 0, D(y)(0) = -50;
Opl := dsolve({DIV, RV}, y(x));



I'd like to know, if it is possible to define any sort of range for parameters in NonlinearFit. E. g. I know that one of parameters should be somewhere between 0.2 - 0.4. I know there is a possibility of initalvalues, but using it doesn't lead into this range.



I am so pleased to be here to look for help in using Maple.

I am trying to find solutions to a set of nonlinear systems of equations using fsolve.

This code I am tried to solve is shown as follow:

I can get one solution for three variables, however, it's clear there are some solutions missed.

My questions are listed as:

1) How do I get all solutions for a fixed variable of sigma?

2) How do I create a  loop to obtain the solutions for various values of sigma?

3) How to plot solution of each variables versus sigma? 

3) How do I determine the stability of each solution?

Many thanks in advance.




Hi Mapleprimes,

I am having difficulties trying to evaluate the following integral

int(-(C__A*K + 1)/(k*C__A), C__A = C__A0 .. C__A)

int(-(C__A*K + 1)/(k*C__A), C__A = C__A .. C__A0, numeric = false)

Where the limits are undefined(unknown) symbolic limits, C_A0 is intial concentration, CA is actual concentration. Maple returns the same equation when I try it. I haven't loaded any packages.

I could of course just evaluate the indefinite version and manually substitute but I think that defeats the purpose of using maple.

Any thoughts will be greatly appreciated.

Maple 2020.1

The worksheet below contains an example of the tiling of the hyperbolic plane.

I would like to produce this and other hyperbolic plane tilings from Maple worksheets, but I don't know the math technique for doing so.

Please direct me to a source of the requisite knowledge or an example worksheet which I can study and imitate.

First 8 9 10 11 12 13 14 Last Page 10 of 55