MaplePrimes Questions

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. 

OEIS A034828 and OEIS A000292 (which give the Wiener index for the cycle graph and the path graph respectively) mention that 

the Wiener index of the cycle of length 19 is 855 and 
the Wiener index of the path with 19 edges is 1330

However, 

GraphTheory:-WienerIndex(GraphTheory:-CycleGraph(19));
 = 
                               38

GraphTheory:-WienerIndex(GraphTheory:-PathGraph(20));
 = 
                               38

So what happened here? 

I fail to understand from the manual how to do total differentiation.

e.g let's say I have a trivial transform defined as

u=t(x,y)-s(x,y);

then I should get

du=diff(t,x)dx+diff(t,y)dy ...etc

Apparently the operator D_Dx should be able to perform the total derivative, but somehow I cant see this is possible.

What do I miss ?

I was attempting to remove part of string using StringTools:-Remove()

But it causes server.exe crash each time.

Any idea why this happens?

Worksheet below

interface(version);

`Standard Worksheet Interface, Maple 2024.0, Windows 10, March 01 2024 Build ID 1794891`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1735. The version installed in this computer is 1732 created 2024, April 17, 19:33 hours Pacific Time, found in the directory C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib\`

restart;

3048

ode:=x*y(x)*diff(y(x), x) = (x + 1)*(y(x) + 1):
ode:=diff(y(x), x$2) = (y(x) + 1):
the_output:=Student:-ODEs:-ODESteps(ode,y(x)):
 

s:=latex(the_output,'output'=string):

StringTools:-Remove(s,"\\begin{array}{ccc}");


Download crash_during_string_tools_remove_maple_2024.mw

ps. reported to Maplesoft also.

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