MaplePrimes Questions

Recently I have wanted to add some questions/answers to "Favourites" . I added some posts a few years ago. Now I can't figure out how do do it. A good example is this evenings question on setting up a library. How to install a library on Maple? - MaplePrimes

The answer has useful links. And this is a particular topic I tend to have difficulty with.

On a Windows machine, I want to install a library in Maple? 

Could you please help me?

Thank you in advance.

Sincerely yours

Dear Collaboratory, as a former mathcad user I am trying to transfer slowly towards Maple. In the attached file I have a problem in solving the ODE numerical. I would appreciate it if someone would be willing to help me out. I tried to upload my maple file but I am not sure it is attached.


I am trying to integrate a function from -infinity to infinity.  When I enter the integrate command into Maple, I just get the command back.  Does anyone know how to properly code the command for an indefinite integral?

I have attached the relevant

When converting a list which includes lists to a set, duplicate entries are dropped.

However converting a list with vectors to a set, no duplicate entries are dropped.

Any explanation to that behaviour?

a := {Vector(2, [1, 0]), Vector(2, [2, 0])}

{Vector[column](%id = 36893490347685388884), Vector[column](%id = 36893490347685389004)}


b := [Vector(2, [1, 0]), Vector(2, [2, 0]), Vector(2, [1, 0])]

[Vector[column](%id = 36893490347747359308), Vector[column](%id = 36893490347747359428), Vector[column](%id = 36893490347747359548)]


c := convert(b, set)

{Vector[column](%id = 36893490347747359308), Vector[column](%id = 36893490347747359428), Vector[column](%id = 36893490347747359548)}


d := [[1, 0], [2, 0], [1, 0]]``

[[1, 0], [2, 0], [1, 0]]


e := convert(d, set)

{[1, 0], [2, 0]}





Please can you give me a hand with numerical solving and visualising the integral equation with stochastic term B(t), using quadrature method ?

I had a look at the "stochastic" package by Sasha Cyganowski

Look forward to your help.



I need help on maple code for solving both linear and non linear boudary condition for fractional order partial differential equation 


I don't understand the number of solutions of my trigonometric equation if i work with solve command ( see maple file)

Ideas ?

I cannot solve the equation. The problem is:

we have the integrals (1) and variable (2)
(1) F := (a1, b1) -> int(1/sqrt(1 + (-1)*b1^2*sin(teta)*sin(teta)), teta = 0 .. a1)

(2) Phi(Pi/2, teta) = 2*(F((Pi/2 - teta_l)/2, sqrt(2/(1 - sin(teta_l)))) - F1(Pi/4, sqrt(2/(1 - sin(teta_L)))))/sqrt(sin(teta_l) - 1)

After the solution for variable (2) I received: 

(-InverseJacobiAM(-Pi/4 + teta_l/2, 1/abs(cos(Pi/4 + teta_l/2))) - InverseJacobiAM(Pi/4, 1/abs(cos(Pi/4 + teta_l/2))))^2/(2*(sin(teta_l) - 1))

Next, I have the equation and (by the followed paper) should solve the equation respect to teta_l

M*B*(L/C)^2/G = 1/8*Phi(Pi/2, teta_l)^2

But I received just empty solution. The reference paper 10.1016/j.jmps.2020.104045

Can somebody recommend the way to solve the problem? 


Student here, and I'm currently using a few year old laptop (amd a4-9125  2-core, 8gb ram) running latest stable ubuntu, and maple keeps crashing on me seemingly out of nowhere, even on very simple tasks, like copy-pasting equations and calculating integrals. This doesn't happend all the time, but when it does, it just shuts the program. I checked the logs and found Maple got this error message:
: free(): invalid next size (normal)
: Aborted (core dumped)

I checked a bit online, and it might be a memory issue.
Is my computer just too old for this program, or is it something fixable?

What's the simplest way to partially substitute for expressions that contain operators from the Physics package? 
(hopefully my worksheet below shows what I mean)




`The "Physics Updates" version in the MapleCloud is 1092 and is the same as the version installed in this computer, created 2021, October 18, 14:11 hours Pacific Time.`


Setup(quantumoperators = {sigma__x, sigma__y, sigma__z, rho, H, comm}, mathematicalnotation = true);

[mathematicalnotation = true, quantumoperators = {H, comm, rho, sigma__x, sigma__y, sigma__z}]


we could just use Physics:-Psigma, etc... but here let's define these manually:

%Commutator(sigma__x, sigma__y) = 2*I*sigma__z,
%Commutator(sigma__y, sigma__z) = 2*I*sigma__x,
%Commutator(sigma__z, sigma__x) = 2*I*sigma__y

[algebrarules = {%Commutator(sigma__x, sigma__y) = (2*I)*sigma__z, %Commutator(sigma__y, sigma__z) = (2*I)*sigma__x, %Commutator(sigma__z, sigma__x) = (2*I)*sigma__y}]


some expression (normally long and complicated that contains many operators):

expr:=(sigma__y*sigma__x*rho - sigma__y*rho*sigma__x);

Physics:-`*`(sigma__y, sigma__x, rho)-Physics:-`*`(sigma__y, rho, sigma__x)


say want to make this type of substitution:


Physics:-`*`(sigma__y, sigma__x) = H


this doesn't work

algsubs(subs_vars, expr);

Physics:-`*`(sigma__y, sigma__x, rho)-Physics:-`*`(sigma__y, rho, sigma__x)


this also doesn't work:

subs(subs_vars, expr);

Physics:-`*`(sigma__y, sigma__x, rho)-Physics:-`*`(sigma__y, rho, sigma__x)


this works (here we explicitly specify the whole thing, that includes rho)

subs(sigma__y*sigma__x*rho=H*rho, expr);

Physics:-`*`(H, rho)-Physics:-`*`(sigma__y, rho, sigma__x)


but I would like for the substitution to also work for *parts* of the expression, as in the normal case, when not using operators.





I need to declare a whole set of variables as local. The variable names are generates algorithmically using assign. Like so:


Stand-alone, this works and creates all these Vectors for later use. But this:

local seq(seq(assign(cat(S,i,j)=Vector(datatype=float)),i=1..9),j=1..9);

does not work; I get an "error; '(' unexpected".

I really do not want to type all these by hand... on the other hand, if I do not declare these as local I get 99 warnings about implicit local declaration; not nice.

Is there a way to do this?



PS: I do not upload as the one line really is all that is needed. At the lowest level one does not get the implicit-declaration warning, but with "local" it still fails.

My forced spring mass system is 4x"+4x'+3x=sin(wt). I calculated my w=w* value that maximizes the amplitude (0.5) and my initial conditions are x(0)=x'(0)=0. I need to graph x(t) when w= w* and when w=w*/2. How am I supposed to input this information into maple to create a graph? 

I drew a polygon using some point coordinates combined with the pointplot.

c1 := <0, -2>:
c2 := <1, 2>:
c3 := <2, 2>:
c4 := <0.5, 6>:
c5:=<1, 4>:
p1:=pointplot([c1, c2, c3, c4,c5, c1], color = red, connect = true);

I ‘d like to add  label  of endpoints to polygons, although I did this via textplot, but it didn't feel very neat. I hope when the polygon is drawn, labels of  every endpoints appear . These labels are the names of the assignment variable for the coordinates of the endpoints.

c1 := <0, -2>:
c2 := <1, 2>:
c3 := <2, 2>:
c4 := <0.5, 6>:
c5:=<1, 4>:

t1 :=textplot([0, -2, 'typeset'("c1")], 'align' = 'above'):
t2 :=textplot([1, 2, 'typeset'("c2")], 'align' = 'above'):
t3 :=textplot([2, 2, 'typeset'("c3")], 'align' = 'above'):
t4 :=textplot([0.5, 6, 'typeset'("c4")], 'align' = 'above'):
t5 :=textplot([1, 4, 'typeset'("c5")], 'align' = 'above'):
p1:=pointplot([c1, c2, c3, c4,c5, c1], color = red, connect = true):
display({p1, t1,t2,t3,t4,t5});


I think  DrawGraph command  in graph theory package  is very good to achieve this in a sense. 

Hello everyone.

Please help me solve this equation for u using Maple.

Div(u) = d(ru)/(dr)*1/r = 2*a = const; u = u(r);


Divergence(u) = d(ru)/(r*dr) and d(ru)/(r*dr) = 2*a and 2*a = const

u = u(r)

Thanks in advance for your help.

First 222 223 224 225 226 227 228 Last Page 224 of 2281