MaplePrimes Questions

I want to solve or try to solve this equation 

PDE := diff(G(a, H, phi, PI), a)(aH) + diff(G(a, H, phi, PI), H)(k/a^2 - kappa^2/2*PI^2/a^6) + diff(G(a, H, phi, PI), phi)(PI/a^3) = diff(G(a, H, phi, PI), PI)(a^3*diff(V(phi), phi))

with pdsolve(PDE, G)

and maple answer me the next

Error, (in pdsolve/info) first argument does not have a differentiated function with name G

I nw in maple, maybe I´m make a mistake, but I can't find what

How can I  get the result for the integration when m =n or m is not equal to n and How can I add assumption that m, or n can be even or odd?

Thanks in advance for your help.

 

restart;
phi := proc(k,x,L)
  if (type(k,even)) then sqrt(2)*sin(Pi*k*x/L)/sqrt(L)
  else sqrt(2)*cos(Pi*k*x/L)/sqrt(L)
  end if;
end proc:

Int(phi(m,x,L)*_h^2/m2*diff(phi(n,x,L),x,x),x=-L/2..L/2);

Int(-2*cos(Pi*m*x/L)*_h^2*Pi^2*n^2*cos(Pi*n*x/L)/(L^3*m2), x = -(1/2)*L .. (1/2)*L)

(1)
 

 

Download test1.mw

 

I have  4 worksheets with derived equations. So I export the equations and  possibly some procedures (but they can be handled seperately if needed)  from each worksheet as a .mpl file. 

I want to combine the .mpl files  together without using copy/paste. Then I can open that single file in the VS code editor.
There may be other ways to achieve this so I am open to suggestions.

Good evening!

I am Athanasios Paraskevopoulos, a graduate student specializing in applied mathematics. Recently, I've started exploring Maple through its trial version and I'm considering making a purchase. My question for you all is: Am I restricted to buying only the graduate student package, or am I free to choose from any of the Maple packages available? Any guidance or personal experiences with different packages would be greatly appreciated!

Thank you in advance for your help!

I chased down a problem to factoring a square that has sqrt in the coefficients. All numbers are real,
The code is inside a procedure in a package. Iso I could do with something robust.

expand((sqrt(A+B)*x+sqrt(7-K)*y)^2)
     2      2            (1/2)          (1/2)        2      2
  A x  + B x  + 2 (A + B)      x (7 - K)      y - K y  + 7 y 

factor(%) 

 

Using a command in succession can coax maple into different outputs. 

eq:=(a*x+b)/(c*x+d)=1
                                
Getting rid of the denominator is probably the most deciding factor in how maple displays an output.

isolate(eq,b)
                      

isolate(%,x,1)
                      

Using the isolate command, you couldn't arrive at that output without using it twice. 

I was trying to solve a system of polynomial equations, which contains three equations and six variables $a_0,a_1,a_2,b_0,b_1,b_2$. However, as I swap the variable name, Maple solve function gives me a totally different solutions. Only the solutions before swapping the variables are useful for the problem I study. I have already attached the file. Could anyone tell me if the choice of variable name really matters, or if i just misuse this function?

Choice_of_name_infolevel.mw

Hello, everyone,

I am new to Maple and I am trying to get use of it.

I tried to plot the following linear systems in different ways. I realized that the Student Linear Algebra is not as flexible as Linear Algebra. My question is the following. Is there any other way to create a plot without defining the implicit plots?

with(Student[LinearAlgebra])

A := Matrix([[1, 1], [12, 16]]); b := Vector([10, 136]); sol := LinearSolve(A, b)

Matrix(2, 2, {(1, 1) = 1, (1, 2) = 1, (2, 1) = 12, (2, 2) = 16})

 

Vector(2, {(1) = 10, (2) = 136})

 

Vector[column](%id = 36893488153382714652)

(1)

LinearSystemPlot({x+y = 10, 12*x+16*y = 136}, axes = normal)

 

restart

 

 

with(Student[LinearAlgebra])

A := Matrix([[2, -1, 1], [0, 1, 3], [0, 0, 1]]); b := Vector([-5, 7, 2]); sol := LinearSolve(A, b)

Matrix(3, 3, {(1, 1) = 2, (1, 2) = -1, (1, 3) = 1, (2, 1) = 0, (2, 2) = 1, (2, 3) = 3, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1})

 

Vector(3, {(1) = -5, (2) = 7, (3) = 2})

 

LinearSolve(Matrix(%id = 36893488151878044716), Vector[column](%id = 36893488151878030628))

(2)

LinearSystemPlot({w = 2, y+3*w = 7, 2*x-y+w = -5}, axes = normal)

 

restart

with(plots); with(LinearAlgebra); A := Matrix([[1, 1], [12, 16]]); b := Vector([10, 136]); sol := LinearSolve(A, b); eq1 := x+y = 10; eq2 := 12*x+16*y = 136; plot1 := implicitplot(eq1, x = -5 .. 10, y = -50 .. 50, color = "red", thickness = 2, labels = ["x", "y"]); plot2 := implicitplot(eq2, x = -5 .. 10, y = -50 .. 50, color = "blue", thickness = 2); display(plot1, plot2, title = "Plot of Linear System", legend = [x+y = 10, 12*x+16*y = 136])

Matrix(2, 2, {(1, 1) = 1, (1, 2) = 1, (2, 1) = 12, (2, 2) = 16})

 

Vector(2, {(1) = 10, (2) = 136})

 

Vector[column](%id = 36893488153330030820)

 

x+y = 10

 

12*x+16*y = 136

 

 

A := Matrix([[2, -1, 1], [0, 1, 3], [0, 0, 1]]); b := Vector([-5, 7, 2]); solution := LinearSolve(A, b); eq1 := 2*x-y+z = -5; eq2 := y+3*z = 7; eq3 := z = 2; plot1 := implicitplot3d(eq1, x = -10 .. 10, y = -10 .. 10, z = -10 .. 10, color = "red", style = surface); plot2 := implicitplot3d(eq2, x = -10 .. 10, y = -10 .. 10, z = -10 .. 10, color = "blue", style = surface); plot3 := implicitplot3d(eq3, x = -10 .. 10, y = -10 .. 10, z = -10 .. 10, color = "green", style = surface); display(plot1, plot2, plot3, title = "3D Plot of Linear System", axes = boxed)

Matrix(3, 3, {(1, 1) = 2, (1, 2) = -1, (1, 3) = 1, (2, 1) = 0, (2, 2) = 1, (2, 3) = 3, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1})

 

Vector(3, {(1) = -5, (2) = 7, (3) = 2})

 

Vector[column](%id = 36893488152610632156)

 

 

NULL


 

Download linear_systems.mw

 

How can I get the roots of this equation: besselJ,  J0(xn)=0?

Multiplication is possible by a numerical constant

a < b;
-1*%;
                             a < b

                            -b < -a

However, using a name with assumptions

with(RealDomain):# just to make sure that a and b are real, probably not required
ineq := a/c < b/c;
assume(c < 0);
c*ineq;
expand(%);
eval(%);
                                 a   b
                         ineq := - < -
                                 c   c

                           /a   b   \
                          *|- < -, c|
                           \c   c   /

                           /a   b   \
                          *|- < -, c|
                           \c   c   /

                           /a   b   \
                          *|- < -, c|
                           \c   c   /

does not work. Have I missed an essential assumption/trick?

What else can I do?

Or: is there a mathematical reasons not to offer this possibilty.

Context: For manual case studies of parameters in inequalties I find it sometimes helpfull to remove rational expressions. Doing this manually with lhs/rhs and denom and flipping < or > is prone to errors

How can numbers be displayed inside the contour plot?

 restart;
with(plots);
contourplot(x*exp(-x^2 - y^2), x = -2 .. 2, y = -2 .. 2, axes = boxed);
like this

Dear All,

I solved the following partial differential equation numerically using pdsolve. As mentioned in the help pages, plotting the solution versus x at a given time is possible. 
Can anyone help me extract the plot of the solution versus time at a given space variable x? Also, how can I extract the data of the mentioned plot, for example, u_numeric (0.75, 0.5)?
pd_numeric:=(D[2,2])(u_numeric)(x,t)+(D[1,1,1,1])(u_numeric)(x,t)-0*h(x,t)=0;
bc_numeric[1]:=u_numeric(0,t)=0;
bc_numeric[2]:=(u_numeric)(1,t)=0; 
bc_numeric[3]:=D[1,1](u_numeric)(0,t)=0;
bc_numeric[4]:=D[1](u_numeric)(1,t)=0; 
ic_numeric[1]:=u_numeric(x,0)=0.1*x*(x-1)^2;
ic_numeric[2]:=D[2](u_numeric)(x,0)=0;

sol:=pdsolve(pd_numeric, {seq(bc_numeric[i], i=1..4), seq(ic_numeric[i],i=1..2)}, u_numeric(x,t), time=t, range=0..1, numeric, spacestep=1/2000, timestep=1/2000);
sol:-plot(t=1);
sol:-animate(t=1, frames=2000, title="time=%f");
Best wishes

I'm working on simulating a triple pendulum. I have the numeric solution of the ODEs and some nice plots that indicate I'm on the right track. I would to animate these solutions. 
It seems like there is a way to plot the positions over time, as well as the lines between the points. If anyone can put me on the right track towards this, I would really appreciate it. 

First 109 110 111 112 113 114 115 Last Page 111 of 2426