Maple Questions and Posts

Does Maple provide bold notation for vectors...

Asked by: C_R 3387

May 01 2025

2 17

I am interested in having bold notation for vectors and matrices.

Any commands or packages that output expressions this way?

Verifying ODE Solution and Issues with Assumptions...

Asked by: salim-barzani 1555

April 30 2025

3 2

I’m trying to test a specific function as a solution to a nonlinear ODE in Maple. The equation is of the Riccati type, and my candidate solution involves parameters A, B, and C.

I've used assuming to specify the condition (4AC−B2)>0 and (4AC - B^2) <0, but when I use odetest to verify the solution, I still get a nonzero result. Additionally, when I apply the assumption, Maple sometimes introduces a negation sign in the output (e.g., changing sqrt(...) into -sqrt(...)), which wasn't part of the original solution.

>

(1)

>

(2)

>

G(xi) = (1/2)*((4*A*C-B^2)^(1/2)*tan((1/2)*(4*A*C-B^2)^(1/2)*(d[0]+xi))-B)/C

(3)

>

(4)

>

G(xi) = -(1/2)*((4*A*C-B^2)^(1/2)*cot((1/2)*(4*A*C-B^2)^(1/2)*(d[0]+xi))+B)/C

(5)

>

(6)

>

G(xi) = -(1/2)*((4*A*C-B^2)^(1/2)*tanh((1/2)*(4*A*C-B^2)^(1/2)*(d[0]+xi))+B)/C

(7)

>

(8)

Download A2.mw

Why the subs command not apply ?...

Asked by: salim-barzani 1555

April 29 2025

3 9

i did every thing coreectly but nothing happen not apply where is my mistake?

>

(1)

>

(2)

>

(3)

>

(4)

>

diff(diff(G(xi), xi), xi) = (1/2)*(2*r^2*G(xi)*(a+b*G(xi)+l*G(xi)^2)*(diff(G(xi), xi))+r^2*G(xi)^2*(b*(diff(G(xi), xi))+2*l*G(xi)*(diff(G(xi), xi))))/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(5)

>

diff(diff(G(xi), xi), xi) = (1/2)*r^2*G(xi)*(diff(G(xi), xi))*(4*l*G(xi)^2+3*b*G(xi)+2*a)/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(6)

>

diff((r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2), xi) = (1/2)*r^2*G(xi)*(4*l*G(xi)^2+3*b*G(xi)+2*a)

(7)

>

diff(diff(G(xi), xi), xi) = (1/2)*r^2*G(xi)*(diff(G(xi), xi))*(4*l*G(xi)^2+3*b*G(xi)+2*a)/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(8)

>

diff(diff(G(xi), xi), xi) = (1/2)*r^2*G(xi)*(diff(G(xi), xi))*(4*l*G(xi)^2+3*b*G(xi)+2*a)/(r^2*G(xi)^2*(a+b*G(xi)+l*G(xi)^2))^(1/2)

(9)

Download subs.mw

Computing the left and right eigenvectors...

Asked by: Elisha 50

April 29 2025

2 1

Someone please help me with the computation of the right and left eigenvectors. my system of equation is attached below

>

:

>

(2)

>

(3)

>

(4)

>

Limit(N(t) = N__0*exp(-mu*t)+exp(mu*t)*K/mu, t = infinity)

limit(N(t), t = infinity) = limit(N__0*exp(-mu*t)+exp(mu*t)*K/mu, t = infinity)

(5)

>

>	#to calculate the disease free equilibrium,

>

(6)

>

(7)

>

(8)

>

(9)

>

(10)

>

>	$solve({E1 = 0, E3 = 0}, {S, S__A})$

${S = Lambda__p/Âµ__C, S__A = Lambda__A/Âµ__A}$

(11)

>

>	#to calculate the Endemic Equilibrium state,

(12)

>

:

>

(14)

>

(15)

>

(16)

>

E1 := `Λ__p`-(`ϕ`*`θ__B`*I__A/N__p+µ__C)*S+`ω__B`*I__B

Lambda__p-(varphi*theta__B*I__A/N__p+Âµ__C)*S+omega__B*I__B

(17)

>

E2 := `ϕ`*`θ__B`*I__A*S/N__p-`ω__B`*I__B-(`σ__B`+µ__C)*I__B

varphi*theta__B*I__A*S/N__p-omega__B*I__B-(sigma__B+Âµ__C)*I__B

(18)

(19)

>

E3 := `Λ__A`-(µ__A+`ϕ`*`α__B`*I__B/N__p)*S__A+`δ__A`*I__A

Lambda__A-(Âµ__A+varphi*alpha__B*I__B/N__p)*S__A+delta__A*I__A

(20)

>

E4 := `ϕ`*`α__B`*I__B*S__A/N__p-(µ__A+`δ__A`)*I__A

varphi*alpha__B*I__B*S__A/N__p-(Âµ__A+delta__A)*I__A

(21)

>

(22)

>

(23)

>	$solve({E1 = 0, E2 = 0, E3 = 0, E4 = 0}, {I__A, I__B, S, S__A})$

${I__A = 0, I__B = 0, S = Lambda__p/Âµ__C, S__A = Lambda__A/Âµ__A}, {I__A = -(N__p^2*Âµ__A^2*Âµ__C^2+N__p^2*Âµ__A^2*Âµ__C*omega__B+N__p^2*Âµ__A^2*Âµ__C*sigma__B+N__p^2*Âµ__A*Âµ__C^2*delta__A+N__p^2*Âµ__A*Âµ__C*delta__A*omega__B+N__p^2*Âµ__A*Âµ__C*delta__A*sigma__B-varphi^2*Lambda__A*Lambda__p*alpha__B*theta__B)/(varphi*Âµ__A*theta__B*(N__p*Âµ__A*Âµ__C+N__p*Âµ__A*sigma__B+N__p*Âµ__C*delta__A+N__p*delta__A*sigma__B+varphi*Lambda__p*alpha__B)), I__B = -(N__p^2*Âµ__A^2*Âµ__C^2+N__p^2*Âµ__A^2*Âµ__C*omega__B+N__p^2*Âµ__A^2*Âµ__C*sigma__B+N__p^2*Âµ__A*Âµ__C^2*delta__A+N__p^2*Âµ__A*Âµ__C*delta__A*omega__B+N__p^2*Âµ__A*Âµ__C*delta__A*sigma__B-varphi^2*Lambda__A*Lambda__p*alpha__B*theta__B)/(alpha__B*(N__p*Âµ__A*Âµ__C^2+N__p*Âµ__A*Âµ__C*omega__B+N__p*Âµ__A*Âµ__C*sigma__B+varphi*Âµ__C*Lambda__A*theta__B+varphi*Lambda__A*sigma__B*theta__B)*varphi), S = (N__p*Âµ__A*Âµ__C+N__p*Âµ__A*sigma__B+N__p*Âµ__C*delta__A+N__p*delta__A*sigma__B+varphi*Lambda__p*alpha__B)*Âµ__A*N__p*(Âµ__C+omega__B+sigma__B)/(alpha__B*varphi*(N__p*Âµ__A*Âµ__C^2+N__p*Âµ__A*Âµ__C*omega__B+N__p*Âµ__A*Âµ__C*sigma__B+varphi*Âµ__C*Lambda__A*theta__B+varphi*Lambda__A*sigma__B*theta__B)), S__A = N__p*(N__p*Âµ__A^2*Âµ__C^2+N__p*Âµ__A^2*Âµ__C*omega__B+N__p*Âµ__A^2*Âµ__C*sigma__B+N__p*Âµ__A*Âµ__C^2*delta__A+N__p*Âµ__A*Âµ__C*delta__A*omega__B+N__p*Âµ__A*Âµ__C*delta__A*sigma__B+varphi*Âµ__A*Âµ__C*Lambda__A*theta__B+varphi*Âµ__A*Lambda__A*sigma__B*theta__B+varphi*Âµ__C*Lambda__A*delta__A*theta__B+varphi*Lambda__A*delta__A*sigma__B*theta__B)/(varphi*Âµ__A*theta__B*(N__p*Âµ__A*Âµ__C+N__p*Âµ__A*sigma__B+N__p*Âµ__C*delta__A+N__p*delta__A*sigma__B+varphi*Lambda__p*alpha__B))}$

(24)

>

Matrix(%id = 18446746854857131062)

(25)

>

$J := Matrix(4, 4, {(1, 1) = -`ϕ`*`θ__B`*I__A/N__p-µ__C, (1, 2) = `ω__B`, (1, 3) = 0, (1, 4) = -`ϕ`*`θ__B`*S/N__p, (2, 1) = `ϕ`*`θ__B`*I__A/N__p, (2, 2) = -`ω__B`-`σ__B`-µ__C, (2, 3) = 0, (2, 4) = `ϕ`*`θ__B`*S/N__p, (3, 1) = 0, (3, 2) = -`ϕ`*`α__B`*S__A/N__p, (3, 3) = -µ__A-`ϕ`*`α__B`*I__B/N__p, (3, 4) = `δ__A`, (4, 1) = 0, (4, 2) = `ϕ`*`α__B`*S__A/N__p, (4, 3) = `ϕ`*`α__B`*I__B/N__p, (4, 4) = -µ__A-`δ__A`})$

Matrix(%id = 18446746579340105118)

(26)

S := `Λ__p`/µ__C

Lambda__p/Âµ__C

(27)

S__A := `Λ__A`/µ__A

Lambda__A/Âµ__A

(28)

(29)

(30)

(31)

$J := Matrix(4, 4, {(1, 1) = -`ϕ`*`θ__B`*I__A/N__p-µ__C, (1, 2) = `ω__B`, (1, 3) = 0, (1, 4) = -`ϕ`*`θ__B`*S/N__p, (2, 1) = `ϕ`*`θ__B`*I__A/N__p, (2, 2) = -`ω__B`-`σ__B`-µ__C, (2, 3) = 0, (2, 4) = `ϕ`*`θ__B`*S/N__p, (3, 1) = 0, (3, 2) = -`ϕ`*`α__B`*S__A/N__p, (3, 3) = -µ__A-`ϕ`*`α__B`*I__B/N__p, (3, 4) = `δ__A`, (4, 1) = 0, (4, 2) = `ϕ`*`α__B`*S__A/N__p, (4, 3) = `ϕ`*`α__B`*I__B/N__p, (4, 4) = -µ__A-`δ__A`})$

Matrix(%id = 18446746579340107518)

(32)

>	$J := Matrix(4, 4, {(1, 1) = -µ__C, (1, 2) = `ω__B`, (1, 3) = 0, (1, 4) = -`β__1`, (2, 1) = 0, (2, 2) = -`ω__B`-`σ__B`-µ__C, (2, 3) = 0, (2, 4) = -`β__1`, (3, 1) = 0, (3, 2) = -`β__2`, (3, 3) = -µ__A, (3, 4) = `δ__A`, (4, 1) = 0, (4, 2) = `β__2`, (4, 3) = 0, (4, 4) = -µ__A-`δ__A`})$

Matrix(%id = 18446746579417403630)

(33)

>

Matrix(%id = 18446746579305905318)

(34)

>

Vector[column](%id = 18446746579445964182)

(35)

>

(36)

>

Vector[column](%id = 18446746579340117046)

(37)

# to find the trace we

$Matrix(7, 7, {(1, 1) = -beta*lambda-v__1-µ, (1, 2) = v__2, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (2, 1) = v__1, (2, 2) = beta*(w-1)*lambda-µ-v__2-alpha, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (3, 1) = 0, (3, 2) = alpha, (3, 3) = -µ, (3, 4) = 0, (3, 5) = `ρ__A`, (3, 6) = `ρ__F`, (3, 7) = -(-1+k)*`ρ__Q`, (4, 1) = beta*lambda, (4, 2) = -beta*(w-1)*lambda, (4, 3) = 0, (4, 4) = -q__E-delta-µ, (4, 5) = 0, (4, 6) = 0, (4, 7) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = a*delta, (5, 5) = -`ρ__A`-q__A-µ, (5, 6) = 0, (5, 7) = k*`ρ__Q`, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = -delta*(-1+a), (6, 5) = 0, (6, 6) = -`ρ__F`-q__F-`δ__F`-µ, (6, 7) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = q__E, (7, 5) = q__A, (7, 6) = q__F, (7, 7) = -`ρ__Q`-`δ__Q`-µ})$

Matrix(%id = 36893490965935089652)

(38)

(39)

#this shows that trace is negative

#to Achieve stability, the value below must be less than zero

(-q__E-delta-Âµ)*(-rho__F-q__F-delta__F-Âµ)*(-k*q__A*rho__Q+q__A*rho__Q+q__A*Âµ+q__A*delta__Q+rho__A*rho__Q+rho__A*Âµ+rho__A*delta__Q+rho__Q*Âµ+Âµ^2+Âµ*delta__Q)*Âµ < 0

(40)

>

(41)

>

${N(t) = Pi/Âµ+exp(-Âµ*t)*_C1}$

(42)

>	$eval({N(t) = Pi/Âµ+exp(-Âµt)_C1}, [t = infinity])$

${N(infinity) = Pi/Âµ+exp(-Âµ*infinity)*_C1}$

(43)

>

${N(infinity) = Pi/Âµ+exp(-Âµ*infinity)*_C1}$

(44)

>

(45)

>

(46)

>

${S(t) = _C1*exp(-(beta*lambda+v__1+Âµ)*t)}$

(47)

Download Cotton_DFE_and_Jacobian.mw

Help with Solving ODE and Representing ± Sign...

Asked by: salim-barzani 1555

April 29 2025

3 8

I tried solving this ODE, but my result is very different from the expected one. How can I correctly obtain the solution? Also, is there a way to include both the positive and negative signs (±) in the equation so that the final result reflects both possibilities?

>

(1)

>

Warning, A new binding for the name `gamma` has been created. The global instance of this name is still accessible using the :- prefix, :-`gamma`. See ?protect for details.

>

(2)

>

(3)

>

(4)

>

G(xi) = (1/2)*(-b+(-4*a*l+b^2)^(1/2))/l, G(xi) = -(1/2)*(b+(-4*a*l+b^2)^(1/2))/l, G(xi) = -4*a*exp(c__1*r*a^(1/2))/(exp(xi*r*a^(1/2))*(4*a*l-b^2+2*b*exp(c__1*r*a^(1/2))/exp(xi*r*a^(1/2))-(exp(c__1*r*a^(1/2)))^2/(exp(xi*r*a^(1/2)))^2)), G(xi) = -4*a*exp(xi*r*a^(1/2))/(exp(c__1*r*a^(1/2))*(4*a*l-b^2+2*b*exp(xi*r*a^(1/2))/exp(c__1*r*a^(1/2))-(exp(xi*r*a^(1/2)))^2/(exp(c__1*r*a^(1/2)))^2))

(5)

>

G(xi) = -4*a*exp(c__1*r*a^(1/2))/(exp(xi*r*a^(1/2))*(4*a*l-b^2+2*b*exp(c__1*r*a^(1/2))/exp(xi*r*a^(1/2))-(exp(c__1*r*a^(1/2)))^2/(exp(xi*r*a^(1/2)))^2))

(6)

>

G(xi) = 4*a*exp(c__1*r*a^(1/2))*exp(xi*r*a^(1/2))/(-4*a*l*(exp(xi*r*a^(1/2)))^2+b^2*(exp(xi*r*a^(1/2)))^2-2*b*exp(c__1*r*a^(1/2))*exp(xi*r*a^(1/2))+(exp(c__1*r*a^(1/2)))^2)

(7)

>

G(xi) = -4*a*exp(a^(1/2)*r*(c__1+xi))/(4*a*l*exp(2*xi*r*a^(1/2))-b^2*exp(2*xi*r*a^(1/2))+2*b*exp(a^(1/2)*r*(c__1+xi))-exp(2*c__1*r*a^(1/2)))

(8)

>

G(xi) = -4*a*(cosh(a^(1/2)*r*(c__1+xi))+sinh(a^(1/2)*r*(c__1+xi)))/(4*a*l*(cosh(2*xi*r*a^(1/2))+sinh(2*xi*r*a^(1/2)))-b^2*(cosh(2*xi*r*a^(1/2))+sinh(2*xi*r*a^(1/2)))+2*b*(cosh(a^(1/2)*r*(c__1+xi))+sinh(a^(1/2)*r*(c__1+xi)))-cosh(2*c__1*r*a^(1/2))-sinh(2*c__1*r*a^(1/2)))

(9)

>

G(xi) = -4*a*(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2)))/((cosh(xi*r*a^(1/2))+sinh(xi*r*a^(1/2)))*(4*a*l-b^2+2*b*(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2)))/(cosh(xi*r*a^(1/2))+sinh(xi*r*a^(1/2)))-(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2)))^2/(cosh(xi*r*a^(1/2))+sinh(xi*r*a^(1/2)))^2))

(10)

>

G(xi) = -4*a*(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2)))*(cosh(xi*r*a^(1/2))+sinh(xi*r*a^(1/2)))/((4*a*l-b^2)*cosh(xi*r*a^(1/2))^2+((8*a*l-2*b^2)*sinh(xi*r*a^(1/2))+2*b*(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2))))*cosh(xi*r*a^(1/2))+(4*a*l-b^2)*sinh(xi*r*a^(1/2))^2+2*b*(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2)))*sinh(xi*r*a^(1/2))-(cosh(c__1*r*a^(1/2))+sinh(c__1*r*a^(1/2)))^2)

(11)

Download tt.mw

Slider Plot and Value Substitution Help...

Asked by: Andiguys 65

April 29 2025

0 4

I need to create a slider plot for A10, A11, and A12 by varying the parameters theta, Pu, and a.
I have a syntax ready — could you suggest modifications to make it work correctly and generate the plot?

Additionally, is it possible to compute the values of A13 and A14 by substituting the obtained A10, A11, and A12 values for each combination of theta, Pu, and a from the slider plot?

Sheet attached: Slider_Q.mw

bug report. dsolve sometimes gives series solution...

Asked by: nm 11353

April 28 2025

1 3

The series to ode using 'series' option (if it exists) should always be series(...), i.e. with big O at end. but sometimes Maple forgets to add this. Here is an example

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

>	ode:=diff(y(x),x)+y(x)=1+x; IC:=y(0)=0; sol:=dsolve([ode,IC],y(x),'series')

>	lprint(sol); # notice solution is not series, it should be

y(x) = x

>	#above solution should be y(x) = series(x+O(x^6),x,6)

>	#this example below is correct ode:=diff(y(x),x)+y(x)=1+x; IC:=y(0)=1; sol:=dsolve([ode,IC],y(x),'series')

>	lprint(sol); #solution is series

y(x) = series(1+1/2*x^2-1/6*x^3+1/24*x^4-1/120*x^5+O(x^6),x,6)

Download bug_report_dsolve_series_april_28_2025.mw

Series of Complex Analysis Worksheets

by: Paras31 230 Product: Maple 2025

April 28 2025

1 0

Hello everyone,

I am creating this post to begin a thread where I will share a series of worksheets on important topics in Complex Analysis, written as part of my notes for my classes. Complex_Analysis_Notes.pdf

The planned sections include:

Section 1: Infinite Series
Section 2: Power Series
Section 3: The Radius of Convergence of a Power Series
Section 4: The Riemann Zeta Function and the Riemann Hypothesis
Section 5: The Prime Number Theorem

Each worksheet will include calculations, plots, and examples using Maple to illustrate key ideas.

✅ I plan to upload one worksheet every week to keep a steady pace and allow time for discussion and feedback between posts.

I hope this thread will be helpful both for learning and for deeper exploration.
Feel free to comment, suggest improvements, or ask questions as I post the materials.

Thank you!

>

restart; interface(imaginaryunit = 'I'); z := I*(1/3); S_N := proc (n) options operator, arrow; sum(z^k, k = 0 .. n) end proc; limit(S_N(n), n = infinity); S_N(10); S_N(100); S_N(1000); with(plots); points := [seq([Re(evalf(S_N(n))), Im(evalf(S_N(n)))], n = 0 .. 50)]; pointplot(points, connect = true, symbol = solidcircle, symbolsize = 10, color = blue, labels = ["Re", "Im"])

proc (n) options operator, arrow; sum(z^k, k = 0 .. n) end proc

>

restart; interface(imaginaryunit = 'I'); z := I*(1/2); S_N := proc (n) options operator, arrow; sum(z^k, k = 0 .. n) end proc; limit(S_N(n), n = infinity); S_N(10); S_N(100); S_N(1000); with(plots); points := [seq([Re(evalf(S_N(n))), Im(evalf(S_N(n)))], n = 0 .. 50)]; pointplot(points, connect = true, symbol = solidcircle, symbolsize = 10, color = blue, labels = ["Re", "Im"])

proc (n) options operator, arrow; sum(z^k, k = 0 .. n) end proc

>

restart; with(plots); interface(imaginaryunit = 'I'); H := proc (N) local n; sum(1/n, n = 1 .. N) end proc; limit(H(n), n = infinity); limit(Re(H(n)), n = infinity); limit(Im(H(n)), n = infinity); harmonic_pts := [seq([n, H(n)], n = 1 .. 100)]; harmonic_plot := pointplot(harmonic_pts, symbol = solidcircle, style = pointline, color = blue, labels = ["n", "Partial Sum Value"], axes = boxed)

>

Download infinite_series.mw

Controlling Variable Distribution in Polynomial Fo...

Asked by: salim-barzani 1555

April 28 2025

2 5

In this work, I do not intend to expand all the variables across the monomials. Instead, I want to restrict the distribution to only the variables x,y,z,tx, y, z, tx,y,z,t, possibly raising them to appropriate powers as needed, until I obtain the desired solution and satisfy the conditions of my PDE tests. However, I am uncertain whether "monomial" is the correct term to use here.

S1.mw

trail-1.mw

why such different result from integrate when pull...

Asked by: nm 11353

April 27 2025

2 3

These two expressions are the same, just pulled minus sign out

But look what happens when integrating them. the anti derivative of one is much more complicated than the other and contains complex numbers and logs. And no matter what I tried, I could not convert the complicated one to look same as the simpler result. Also could not verify the complicated one by back differentiating.

>	integrand_1:=x^2(-arctan(x) + x)exp(-arctan(x) + x)/(x^2 + 1);

x^2*(-arctan(x)+x)*exp(-arctan(x)+x)/(x^2+1)

>	integrand_2:=evala(integrand_1);

-x^2*(arctan(x)-x)*exp(-arctan(x)+x)/(x^2+1)

>	simplify(integrand_1 - integrand_2)

>	anti_1:=int(integrand_1,x);

>	anti_2:=int(integrand_2,x);

-(1-x+((1/2)*I)*ln(1-I*x)-((1/2)*I)*ln(1+x*I))*(1-I*x)^(-(1/2)*I)*(1+x*I)^((1/2)*I)*exp(x)

>	simplify(diff(anti_1,x)-integrand_1);

>	simplify(diff(anti_2,x)-integrand_2);

Error, (in simpl/simpl/ReIm/sum) too many levels of recursion

>	simplify(anti_1 - anti_2)

Error, (in simpl/simpl/ReIm/sum) too many levels of recursion

>	simplify(anti_2);

(1/2)*(I*ln(1+x*I)-I*ln(1-I*x)+2*x-2)*(1-I*x)^(-(1/2)*I)*(1+x*I)^((1/2)*I)*exp(x)

>	simplify(anti_2,ln);

(1/2)*(I*ln(1+x*I)-I*ln(1-I*x)+2*x-2)*(1-I*x)^(-(1/2)*I)*(1+x*I)^((1/2)*I)*exp(x)

Download int_strange_result_april_27_2025.mw

I would have expected same anti derivative to show. To check, I used another software, and that one gave same anti-derivative for both integrands.

The questions I have: Why Maple gives such different result for same integrand? And how could one convert the one with the logs and complex numbers to the first one?

Maple 2025

Compute integral exactly ...

Asked by: Zeineb 540

April 27 2025

0 0

Dear all
I have a double integral, i want to compute this integral and verify if the pproposed solution verify the proposed equation or not.
I can modify the right hand side of my equation or the exact solution, so that my equation has an exact solution with simple form of right hand side.

exact_solution.mw

Thank you for your help

Plotting contour integral diagrams...

Asked by: Paras31 230

April 27 2025

1 7

I am writing notes on complex analysis, I need to use figures of contour paths to integrate on them, i want to create something like this

I tried to plot the contour for
\oint_{|z|=2} \frac{1}{z^2+1}\,dz
I need to have connecting lines all around because the poles can not be isolated

Download CIF.mw

about updating the Physics Package...

Asked by: Jean-Michel 40

April 27 2025

0 2

Hello all,

After updating the Physics package I have this error :

Physics:-Version();
The "Physics Updates" version "1862" is installed in the

directory C:\Users\jm\maple\toolbox\2025\Physics Updates but

is not active. The active version of Physics is within the

library C:/Users/jm/maple/toolbox/2025/Physics Updates/lib\Ph\

ysics Updates.maple.

What am I supposed to do next?

Thanks a lot and kind regards to all,

Jean-Michel

bug report. Maple 2025 server crash on call to ode...

Asked by: nm 11353

April 27 2025

2 1

FYI;

You might have to try the command more than one time to see the above crash. Here is the worksheet

>

restart;

>	interface(version);

>	Physics:-Version();

>	SupportTools:-Version();

>

restart;

>	ode:=x^2-2xy(x)+5y(x)^2 = (x^2+2xy(x)+y(x)^2)diff(y(x),x);

>

sol:=y(x) = (-1/2*exp(RootOf(-exp(_Z)^2*ln(x*(exp(_Z)-2))+2*_C7*exp(_Z)^2+_Z*exp(_Z)^2+4*exp(_Z)*ln(x*(exp(_Z)-2))-8*_C7*exp(_Z)-4*exp(_Z)*_Z-2*exp(_Z)-4*ln(x*(exp(_Z)-2))+8*_C7+4*_Z+6))^2+3*exp(RootOf(-exp(_Z)^2*ln(x*(exp(_Z)-2))+2*_C7*exp(_Z)^2+_Z*exp(_Z)^2+4*exp(_Z)*ln(x*(exp(_Z)-2))-8*_C7*exp(_Z)-4*exp(_Z)*_Z-2*exp(_Z)-4*ln(x*(exp(_Z)-2))+8*_C7+4*_Z+6))-6+2*(exp(RootOf(-exp(_Z)^2*ln(x*(exp(_Z)-2))+2*_C7*exp(_Z)^2+_Z*exp(_Z)^2+4*exp(_Z)*ln(x*(exp(_Z)-2))-8*_C7*exp(_Z)-4*exp(_Z)*_Z-2*exp(_Z)-4*ln(x*(exp(_Z)-2))+8*_C7+4*_Z+6))^2-6*exp(RootOf(-exp(_Z)^2*ln(x*(exp(_Z)-2))+2*_C7*exp(_Z)^2+_Z*exp(_Z)^2+4*exp(_Z)*ln(x*(exp(_Z)-2))-8*_C7*exp(_Z)-4*exp(_Z)*_Z-2*exp(_Z)-4*ln(x*(exp(_Z)-2))+8*_C7+4*_Z+6))+9)^(1/2))/(1/2*exp(RootOf(-exp(_Z)^2*ln(x*(exp(_Z)-2))+2*_C7*exp(_Z)^2+_Z*exp(_Z)^2+4*exp(_Z)*ln(x*(exp(_Z)-2))-8*_C7*exp(_Z)-4*exp(_Z)*_Z-2*exp(_Z)-4*ln(x*(exp(_Z)-2))+8*_C7+4*_Z+6))^2-3*exp(RootOf(-exp(_Z)^2*ln(x*(exp(_Z)-2))+2*_C7*exp(_Z)^2+_Z*exp(_Z)^2+4*exp(_Z)*ln(x*(exp(_Z)-2))-8*_C7*exp(_Z)-4*exp(_Z)*_Z-2*exp(_Z)-4*ln(x*(exp(_Z)-2))+8*_C7+4*_Z+6)))*x:

>	odetest(sol,ode);

>

Download crash_maple_2025_april_27_2025.mw

Hopefully a fix could be found for this.

Initialization or Library Archive...

Asked by: Christopher2222 6030

April 27 2025

0 5

When generating a file to update parts of maple (for example constants) is it best to put it in an initialization file or make a library archive .mla?

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

E-Mail Address:
Password:
Remember Me:	Automatically sign in on future visits

Ask a Question

Create a Post