Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

How to get the solution for the  given differential equation,in Variational Iteration Method.I tried but i got error.Is this maple code is correct?

Ics:  u[0](x)=sin (x) and take alpha=1, x belongs to R.

Here is my maple code

When adding singsol=all to this ode in the first example, dsolve returns nothing. But it should have returned the general solution if it can not find singular solution. This is what it does always as can be seen from the second example given below.

Why did dsolve return nothing in the first example? Is this a bug?



`Standard Worksheet Interface, Maple 2023.0, Windows 10, March 6 2023 Build ID 1689885`



ode:=y(x)=x*diff(y(x),x)+ a*diff(y(x),x)/(sqrt(1+diff(y(x),x)^2));

y(x) = x*(diff(y(x), x))+a*(diff(y(x), x))/(1+(diff(y(x), x))^2)^(1/2)

y(x) = x*c__1+a*c__1/(c__1^2+1)^(1/2)

"sol_2 := "

ode:=diff(y(x),x)^2-(1+2*x*y(x))*diff(y(x),x)+2*x*y(x) = 0;
sol_2:=dsolve(ode,y(x), singsol=all)

(diff(y(x), x))^2-(1+2*y(x)*x)*(diff(y(x), x))+2*y(x)*x = 0

y(x) = c__1*exp(x^2), y(x) = x+c__1

y(x) = c__1*exp(x^2), y(x) = x+c__1



Hello everyone, I have observed the trend between two line plot is decreasing means the peak of one value is large then peak of other value or we can say that the values is decreasing. I want to observe decreaing trend in percentage means how can we calculate that how much percent one value is decreases from other vlaue.

For example like this way... Figure demonstrates 17% or 30% or 40 % or 50 % decrease in the value of p1 as compared to the value of p0 by change the value of M.

I have a question regarding the combstruct package, to be more precise, the specification of a combinatorial species as it is described here:

Is there any available constructor to generate the functorial composition of two species? For example, a simple graph can be seen as a subset of the set of all 2-element subsets of a given node set. This is nothing else than the functorial composition of E*E with E_2*E, where E denotes the species of sets and E_2 the species of 2-element sets. Thus, the species of subsets can be expressed as E*E and the one of 2-element subsets as E_2*E. Written in the Magma syntax, it should be the functorial composition of Prod(Set(Z),Set(Z)) with Prod(Set(Z,card=2),Set(Z)). Applying the count-function, it should produce the series

How can I achieve that?

If you've seen Paulina's announcement then you know that we are once again holding a virtual Maple Conference this year.  As well, we are once again going to have a virtual gallery featuring artwork and creative projects submitted by the Maple community!

Last year we had a number of great submissions to our Maple Art Gallery and our Maple Learn Creative Showcase.  These were our excellent prize winners.

From left to right we have A visualization of all the primitive roots of 10037 created by Simon Plouffe, winner of the Judge’s Choice, Mother’s Day Rose created with Maple plots by Greg Wheaton, winner of the People’s Choice, and Mona Lisa in Maple Learn created by Paul DeMarco (with help from Leonardo DaVinci), the winner of the People’s Choice for the Maple Learn Showcase.

This year we are expanding the Gallery into two collections to encourage more people to submit.  They are

  • The Art Gallery - A small gallery to highlight high effort, mathematically interesting works (with stricter criteria)

  • The Creative Works Showcase - A larger showcase for nearly any interesting visual works created with Maplesoft products like Maple Learn and Maple

Feel free to submit nearly anything cool for the Creative Works Showcase, if we find it particularly impressive we might even ask you to let us consider it for the gallery.  Also, do not be intimidated by the title "Art Gallery" we're looking for anything that has taken some artistic effort and tells a mathematical story.

For more information on critera and how to submit, please visit our Call for Creative Works.  The important deadline to know is the September 14th deadline for submission of works with virtual gallery reception and awards ceremony durring the conference October 26-27.

I look forward to seeing all the submissions for the Maple community again this year!

How to find series values.I got this error.Please Help.

Maple code for the problem is

Hello everyone,

Has someone encountered this error while tempting

to install the Physics Package in Maple2023.

It is NOT Edgardo's fault but comes from Maple itself!

The same problem occured in Maple 2021 and it took...9 months

to be settled.

I already posted here around one month ago.

I also sent a mail to @Support twice.

Here is the attached file (in French but nobody will miss the point)

Kind regards to all,



Edit :Here's a Copy/Paste of the error (since I do not see my screenshot)

Récupération du package de l’infonuage...
ID: Physics Updates
Version: 5137472255164416
URL: 1460

La taille du fichier est 10326016 octets(9 MB).

Installation du package...
PackageTools:-Install("C:/Users/jm/AppData/Local/Temp/cloudDownload16751660190612963792/Physics Updates_1686091821883.maple",overwrite=true,pkgname="Physics Updates");

2, "failed to resolve filename C:\Maths\Maple2023\toolbox\2023"
ERREUR: Le package n’a pas pu être installé.
error PackageTools:-Install, "permission denied"

Thanks for your time.


i have a question. can maple recognize similar expressions in an equation? for example can maple recognize these two silimiar expressions in the picture below and give it ? I want maple itself recognize similar expressions, not i give an input to find similiar results. thnx in advance

A chordal graph is one in which all cycles of four or more vertices have a chord, which is an edge that is not part of the cycle but connects two vertices of the cycle. A perfect elimination ordering in a graph is an ordering of the vertices of the graph such that, for each vertex v, v and the neighbors of v that occur after v in the order form a clique.  A graph is chordal if and only if it has a perfect elimination ordering.

I  use IsChordal  to test whether the lexicographic product of  two graphs g1,g2 is a chordal graph. It returned true and provided a perfect elimination sequence 1, 2, ..., 30. However, vertices of s  are "1:1", "1:2", ..., "10:3", rather than using Arabic numerals. Therefore, it is difficult for me to extract useful information from the perfect elimination sequence.




true, [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]





["1:1", "1:2", "1:3", "2:1", "2:2", "2:3", "3:1", "3:2", "3:3", "4:1", "4:2", "4:3", "5:1", "5:2", "5:3", "6:1", "6:2", "6:3", "7:1", "7:2", "7:3", "8:1", "8:2", "8:3", "9:1", "9:2", "9:3", "10:1", "10:2", "10:3"]



Hi senior, i am trying to customize the range for horizontal axis means in my case y axis but got error, also trying to give different color for outcomes. Please help in this regards.

Is it possible to use a vector in the piecewise fuction? I have tried the following and other stuff but the the error below

timestep := <seq(1 .. 100)>;
piecewise(`<`~(30, timestep), 1);
Error, (in PiecewiseTools:-Is) wrong kind of parameters in piecewise

I have a time vector and would like to create a u(t), unit step fuction to "turn on" a function at a certain time. e.g.


Where u(t)=1 for t>turn_on_time and 0 otherwise.

TIA for help

What  the  error "could not store"means and how can fix it? a peace of file attached.

  Could you please give maple code for Shooting Method/Kellor box Method or any numerical method except RK Method to this system of Differential Equations.

In a recent answer I posted, I had a relation of the form

I*Int(f(x), x) = something - 2*I*Int(f(x), x)

and I wanted to isolate the term Int(f(x), x).
The function isolate failed to do it and I was forced to use some workaround to do the "isolation".

Trying to understand what happened here, it seems that isolate fails when the term to isolate is multiplied by the imaginary unit
Here are a few examples

expr := I*(Int(x^2*ln(-x+sqrt(x^2-1)), x)) = g(x) -(2*I)*(Int(x^2*ln(-x+sqrt(x^2-1)), x))

I*(Int(x^2*ln(-x+(x^2-1)^(1/2)), x)) = g(x)-(2*I)*(Int(x^2*ln(-x+(x^2-1)^(1/2)), x))


# no isolation

isolate(expr, lhs(expr))

I*(Int(x^2*ln(-x+(x^2-1)^(1/2)), x)) = g(x)-(2*I)*(Int(x^2*ln(-x+(x^2-1)^(1/2)), x))


# isolation

expr_1 := expand(expr / I)
isolate(expr_1, lhs(expr_1))

Int(x^2*ln(-x+(x^2-1)^(1/2)), x) = -((1/3)*I)*g(x)


# no isolation neither, so the problem is not related to "Int"

expr := I*diff(h(x), x) = g(x) -2*I*diff(h(x), x):
isolate(expr, lhs(expr))

I*(diff(h(x), x)) = g(x)-(2*I)*(diff(h(x), x))


# no isolation neither, so the problem comes from "I"

expr := I*A = g(x) -2*I*A:
isolate(expr, lhs(expr))

I*A = g(x)-(2*I)*A


# isolation (of course)

expr := c*A = g(x) -2*c*A:
isolate(expr, lhs(expr))

c*A = (1/3)*g(x)




I guess this is a known behavior, but why it is so?
Is there a way to force the "isolation" without using a trick like in result (3)

Thanks in advance

Hello community, I am new in this forum and sorry if the following seems rudimentary:

I am replicating a distribution function using Statistics[Distribution] as defined here:

Vose Software's Beta Subjective

I have tested numerically and the function works for this initial values: Min := 3; Mlikely := 8; Avg := 9; Max := 18;

But failed to complete with these values (which I tested works in another software) Min := 1000; Mlikely := 1400; Avg := 1500; Max := 2100;

What puzzles me is the inconsistency of handling the functions (be aware that the server kernel could be slow at times). I hope I don't have to give up and continue my custom distribution project in another software.

Here is the document: 


The Beta Subjective Distribution







Min := 3; Mlikely := 8; Avg := 9; Max := 18




Mid := (Min+Max)*(1/2)



alpha := 2*(Avg-Min)*(Mid-Mlikely)/((Avg-Mlikely)*(Max-Min))



beta := alpha*(Max-Avg)/(Avg-Min)





f := simplify(piecewise(Min <= x and x <= Max, (x-Min)^(alpha-1)*(Max-x)^(beta-1)/(Beta(alpha, beta)*(Max-Min)^(alpha+beta-1)), 0))

piecewise(x < 3, 0, x <= 18, (4/16875)*(x-3)*(-18+x)^2, 18 < x, 0)




MD := Distribution(PDF = unapply(f, x), Conditions = [`and`(Min < Mlikely and Mlikely < Max and Min < Avg and Avg < Max and Mlikely <> Avg, piecewise(Avg < Mlikely, Mlikely > Mid, Mid > Mlikely))])








X := RandomVariable(MD)












evalf(Quantile(X, .25))




A := Sample(X, 10^5)
















Quantile(A, .25)





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