Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

When I type "with(combinat)" in Maple Linux, I got the following erros message:

"Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received combinat"

So how to make it work? 

 

Many thanks.

code_(1).mw
 

 

 

Hello! hope you all be fine. i m running a code for `Anew Technique of intial Boundary Value Problem Using ADM `(Elaf*Jaafar*Ali) but my code is for positive integer and i want to run it for fractional can anyone help me please code file is attached 

The pdf shows the derivation procedure of Image from 1.2.1 to 1.2.5. However, when we obtain the (7) in pdf, which is the 

 

expansion of 1.2.4, we hope to see the hamiltonian operator in the output.  We obtain the (10) in pdf (Maple). Actually, it is the expansion of 1.2.5. However, it is too cockamamie compared with that (1.2.5) in the textbook.

 

Hence, I hope to show the hamiltonian operator in the (10) in PDF. I hope to see the simplest type of equation.

How could I do? Appreciate the kind help from anyone!

The Maple code is also uploaded.
 

 

 

 

restart*with(Physics[Vectors])

Step(mathematicalnotation = true)*SetCoordinates('cartesian'[i, j, k])

Step(mathematicalnotation = true)*SetCoordinates(cartesian[i, j, k])

(1)

``

NULL``

"  UU(x,y,z):=ux(x,y,z)*(_i)+uy(x,y,z)*_j+uz(x,y,z)*_k; "

proc (x, y, z) options operator, arrow, function_assign; Physics:-Vectors:-`+`(Physics:-Vectors:-`+`(ux(x, y, z)*_i, uy(x, y, z)*_j), uz(x, y, z)*_k) end proc

(2)

"∇(pressure(x,y,z));  Laplacian(UU(x,y,z));  BB(x,y,z):=0*_i+0*_j+10*_k"

proc (x, y, z) options operator, arrow, function_assign; Physics:-Vectors:-`+`(Physics:-Vectors:-`+`(0*_i, 0*_j), 10*_k) end proc

(3)

-VectorCalculus[Nabla](pressure(x, y, z))+Laplacian(mu*UU(x, y, z))+BB(x, y, z) = 0*_i+0*_j+0*_k

_i*(-(diff(pressure(x, y, z), x))+mu*(diff(diff(ux(x, y, z), x), x))+mu*(diff(diff(ux(x, y, z), y), y))+mu*(diff(diff(ux(x, y, z), z), z)))+_j*(-(diff(pressure(x, y, z), y))+mu*(diff(diff(uy(x, y, z), x), x))+mu*(diff(diff(uy(x, y, z), y), y))+mu*(diff(diff(uy(x, y, z), z), z)))+_k*(-(diff(pressure(x, y, z), z))+mu*(diff(diff(uz(x, y, z), x), x))+mu*(diff(diff(uz(x, y, z), y), y))+mu*(diff(diff(uz(x, y, z), z), z))+10) = 0

(4)

``

Divergence(-VectorCalculus[Nabla](pressure(x, y, z))+Laplacian(mu*UU(x, y, z))+BB(x, y, z) = 0*_i+0*_j+0*_k)

-(diff(diff(pressure(x, y, z), x), x))+mu*(diff(diff(diff(ux(x, y, z), x), x), x))+mu*(diff(diff(diff(ux(x, y, z), x), y), y))+mu*(diff(diff(diff(ux(x, y, z), x), z), z))-(diff(diff(pressure(x, y, z), y), y))+mu*(diff(diff(diff(uy(x, y, z), x), x), y))+mu*(diff(diff(diff(uy(x, y, z), y), y), y))+mu*(diff(diff(diff(uy(x, y, z), y), z), z))-(diff(diff(pressure(x, y, z), z), z))+mu*(diff(diff(diff(uz(x, y, z), x), x), z))+mu*(diff(diff(diff(uz(x, y, z), y), y), z))+mu*(diff(diff(diff(uz(x, y, z), z), z), z)) = 0

(5)

 

algsubs(Divergence(UU(x, y, z)) = 0, %)

-(diff(diff(pressure(x, y, z), x), x))+mu*(-(diff(diff(diff(uy(x, y, z), x), x), y))-(diff(diff(diff(uz(x, y, z), x), x), z)))+mu*(-(diff(diff(diff(uy(x, y, z), y), y), y))-(diff(diff(diff(uz(x, y, z), y), y), z)))+mu*(-(diff(diff(diff(uy(x, y, z), y), z), z))-(diff(diff(diff(uz(x, y, z), z), z), z)))-(diff(diff(pressure(x, y, z), y), y))+mu*(diff(diff(diff(uy(x, y, z), x), x), y))+mu*(diff(diff(diff(uy(x, y, z), y), y), y))+mu*(diff(diff(diff(uy(x, y, z), y), z), z))-(diff(diff(pressure(x, y, z), z), z))+mu*(diff(diff(diff(uz(x, y, z), x), x), z))+mu*(diff(diff(diff(uz(x, y, z), y), y), z))+mu*(diff(diff(diff(uz(x, y, z), z), z), z)) = 0

(6)

simplify(%)

-(diff(diff(pressure(x, y, z), x), x))-(diff(diff(pressure(x, y, z), y), y))-(diff(diff(pressure(x, y, z), z), z)) = 0

(7)

Laplacian(-VectorCalculus[Nabla](pressure(x, y, z))+Laplacian(mu*UU(x, y, z))+BB(x, y, z) = 0*_i+0*_j+0*_k)

(mu*(diff(diff(diff(diff(ux(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), z), z), z), z))-(diff(diff(diff(pressure(x, y, z), x), x), x))-(diff(diff(diff(pressure(x, y, z), x), y), y))-(diff(diff(diff(pressure(x, y, z), x), z), z)))*_i+(mu*(diff(diff(diff(diff(uy(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), z), z), z), z))-(diff(diff(diff(pressure(x, y, z), x), x), y))-(diff(diff(diff(pressure(x, y, z), y), y), y))-(diff(diff(diff(pressure(x, y, z), y), z), z)))*_j+(mu*(diff(diff(diff(diff(uz(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), z), z), z), z))-(diff(diff(diff(pressure(x, y, z), x), x), z))-(diff(diff(diff(pressure(x, y, z), y), y), z))-(diff(diff(diff(pressure(x, y, z), z), z), z)))*_k = 0

(8)

algsubs(-(diff(pressure(x, y, z), x, x))-(diff(pressure(x, y, z), y, y))-(diff(pressure(x, y, z), z, z)) = 0, %)

(mu*(diff(diff(diff(diff(ux(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(ux(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(ux(x, y, z), z), z), z), z)))*_i+(mu*(diff(diff(diff(diff(uy(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uy(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uy(x, y, z), z), z), z), z)))*_j+(mu*(diff(diff(diff(diff(uz(x, y, z), x), x), x), x))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), x), x), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), y), y), y), y))+2*mu*(diff(diff(diff(diff(uz(x, y, z), y), y), z), z))+mu*(diff(diff(diff(diff(uz(x, y, z), z), z), z), z)))*_k = 0

(9)

collect(%, mu)

((diff(diff(diff(diff(ux(x, y, z), x), x), x), x)+2*(diff(diff(diff(diff(ux(x, y, z), x), x), y), y))+2*(diff(diff(diff(diff(ux(x, y, z), x), x), z), z))+diff(diff(diff(diff(ux(x, y, z), y), y), y), y)+2*(diff(diff(diff(diff(ux(x, y, z), y), y), z), z))+diff(diff(diff(diff(ux(x, y, z), z), z), z), z))*_i+(diff(diff(diff(diff(uy(x, y, z), x), x), x), x)+2*(diff(diff(diff(diff(uy(x, y, z), x), x), y), y))+2*(diff(diff(diff(diff(uy(x, y, z), x), x), z), z))+diff(diff(diff(diff(uy(x, y, z), y), y), y), y)+2*(diff(diff(diff(diff(uy(x, y, z), y), y), z), z))+diff(diff(diff(diff(uy(x, y, z), z), z), z), z))*_j+(diff(diff(diff(diff(uz(x, y, z), x), x), x), x)+2*(diff(diff(diff(diff(uz(x, y, z), x), x), y), y))+2*(diff(diff(diff(diff(uz(x, y, z), x), x), z), z))+diff(diff(diff(diff(uz(x, y, z), y), y), y), y)+2*(diff(diff(diff(diff(uz(x, y, z), y), y), z), z))+diff(diff(diff(diff(uz(x, y, z), z), z), z), z))*_k)*mu = 0

(10)

``

``

 


``

 

Download 20200619.mw

 

 

Download 20200619-1.mw

20200619.pdf

 

Error, (in Compiler:-Setup) unable to write to initialization file, C:\Program Files\Maple 2020\bin.X86_64_WINDOWS/launch.ini, possibly due to insufficient permissions
NULL;

when choose bat file

but i see sufficient permission

any w should i s

choose a particular compiler or even visual studio i have in my system

 

where can i find if anything else

Hello,

I try import a big Excel file with the first row as column labels for performing statistical etc.

I use this command :

data := Import(FileTools:-JoinPath(["d:", "directory", "subset2.xlsx"], platform = "windows"))

i need a dataframe for performing statistical with column label as first row, but i have this error :

Error, (in Import) invalid input: MatrixToDataFrame uses a 3rd argument, columns, which is missing

I can not understand what is this error...so, what i can do ?

Thanks, regards


This is my notebook

 

restart

with(DifferentialGeometry)

DGsetup([t, rho, theta, phi, `#mover(mi("θ",fontstyle = "normal"),mo("&uminus0;"))`, `#mover(mi("φ",fontstyle = "normal"),mo("&uminus0;"))`, `#mover(mi("θ",fontstyle = "normal"),mo("˜"))`, `#mover(mi("φ",fontstyle = "normal"),mo("˜"))`], M)

M > 

``

M > 

evalDG(`&w`(dt, `dρ`))

_DG([["form", M, 1], [[[1], `dρ`]]])

(1)
M > 

evalDG(`&s`(dt, `dρ`))

_DG([["tensor", M, [["cov_bas"], []]], [[[1], `dρ`]]])

(2)
M > 

 

 

 

``

I don't understand why computing the wedge product of two differential forms gives me an output which prints the product of both. Why  don't it just return d rho ^ dt?

Dear Maple Users @acer @Carl Love @Kitonum @Preben Alsholm ,

Greetings.

How to plot a function "Am" for various values of "kt"(eg: kt=-1..1) at a point x=0.

Am vs kt.

The code has provided below.

Waiting for your response.

ktplot.mw

 

Hi experts

with(plots):
#epsilon:=5*Pi/2:
gammaa:=0.4:
gammao:=0.1:
gammab:=0.002:
omegab:=100:
X:=60:
Omega:=6*Pi:
g:=10:
sigmag:=5*gammaa:
a:=1/(sigmag*sqrt(2*Pi)):

#deltao:=1:
D7:=evalf(a*int(g*Omega*(gammao/2-I*Delta)*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2),Delta=-infinity..infinity)):
D8:=evalf(0.5*a*int(g^2*gammao*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
D9:=evalf(a*int(g^2*Delta*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
b1 := (deltao-2*X^2*omegab*Y/(gammab^2+omegab^2)-D9)^2:

 f:=Y*((gammaa+D8)^2+(b1))=Im(D7)^2+(epsilon-Re(D7))^2:
#P1:=implicitplot(f,deltao=-70..50,Y=0..1,numpoints=100,axes=boxed,thickness=2,color=black,font=[1,1,20],tickmarks=[4, 3],linestyle=1);

implicitplot3d(f,Y = 0 .. 1,deltao=-40..40,epsilon=0..0.4,numpoints=100,labels=[Y,deltao,E],tickmarks=[3,3,3], style = surface, color = "Niagara Azure");

The above code to plot the implicitplot3d ,

first the figures are not smooth

second thing is how to make contourplot at epsilon= 5*Pi/2.

I appreciated any comments


 

Hello everyone,

I would have a question about how to colour a region of complicated inequalities.

I'm curious about the region where my 'lambdas' are all negatives regarding the parameters 'pm2' and 'am2', but I always find difficulties even if I try with implicitplot or inequal commands. Also, I would like to put into the same graph with different colours when they are separately smaller than 0.

I've already tried plenty of things, I hopefully attached some pics by snipping tool about them.  And here is an example of my problem:

p61w:= implicitplot([lambda61 <= 0], pm2 = 0 .. 1, am2 = 0 .. 1)

in61 := inequal(lambda61 <= 0, pm2 = 0 .. 1, am2 = 0 .. 1)

 

display(in61,p61w)

 

 

So, the main issue that it is not just coloured where my condition is true, but elsewhere too. In plot 'in61' as you can see, there is a region approx. pm2=0..0.6 where it shouldn't be coloured since my condition is not true there. In the graph; 'display(in61, p61w)',  it is visible that under the red line should be just coloured however my 'inequal' command coloured it outside too. 

Also, I tried with implicitplot, for example:

p6 := implicitplot([lambda61 <= 0], pm2 = 0 .. 1, am2 = 0 .. 1, filled = true, coloring = [blue, white], transparency = .5)

But the result is similarly wrong. And I'm not sure why it happens and how can I fix it.

(It's okay if for example for these particular lambdas (lambda61,lambda62) never negatives together, but I have more systems and I don't think it works properly)

I read similar questions too but I couldn't find an answer so, I would very appreciate any suggestions.

This is my first post in Maple, so please let me know if I have mistaken something.

Mapple doesn' work at document mode.How can i solve the problem.(My Os is Win 10 )

Hi!
I am trying to solve ODE with bvp[midrich], and I need to prove the validity of this method. 
Where can i read about realization of bvp[midrich]. What scheme does it use? Аre there any articles on this subject?
 

 

Hello everyone,

Please how can I caluclate an integrale based on data set (Q) usiing any integral method

thanks


 

restart:

NULL

NULL

NULL

Q := [21521.7, 11966.6, 7859.78, 7947.54, 9667.45, 10206.8, 10934.0, 10925.0, 11151.5, 10964.6, 11009.0, 11087.0, 11164.3, 11215.7, 11272.7, 11319.0, 11354.6, 11386.2, 11413.5, 11436.3, 11456.0, 11473.2, 11488.2, 11501.1, 11512.7, 11523.0, 11532.2, 11540.3, 11547.7, 11554.3, 11560.3, 11565.8, 11570.7, 11575.3, 11579.4, 11583.2, 11586.7, 11590.0, 11592.9, 11595.7, 11598.2, 11600.6]

[21521.7, 11966.6, 7859.78, 7947.54, 9667.45, 10206.8, 10934.0, 10925.0, 11151.5, 10964.6, 11009.0, 11087.0, 11164.3, 11215.7, 11272.7, 11319.0, 11354.6, 11386.2, 11413.5, 11436.3, 11456.0, 11473.2, 11488.2, 11501.1, 11512.7, 11523.0, 11532.2, 11540.3, 11547.7, 11554.3, 11560.3, 11565.8, 11570.7, 11575.3, 11579.4, 11583.2, 11586.7, 11590.0, 11592.9, 11595.7, 11598.2, 11600.6]

(1)

``

NULL


 

Download simpson.mw

Is it possible to import 3rd party python packages like pyGPGO (A  Bayesian Optimization package) available in github. Please reply. Thanks a lot 

The uploaded worksheet contains two examples of the use of VectorCalculus[VectorSpace]. The first example seems explicable but the second does not.

I have tried and failed to find a full, clear explanation of how a vector describing a simple vector, spacecurve, or surface in the default vector space is transformed to appear in a user defined vector space.

Can anyone direct me to such an explanation, so that I can understand Maple's processing within the uploaded worksheet and enable me to use this Maple feature to future advantage?

VectorSpaceTest.mw

How can I generate a code to plot the optimals; h, chi and psi?

2 3 4 5 6 7 8 Last Page 4 of 1708