Advanced Search

Maple 2024 Questions and Posts

These are Posts and Questions associated with the product, Maple 2024
View unanswered posts
Maple 2024 Questions and Posts Feed

How determine the ode type ? my answer and real an...

Asked by: Arya-S-AA 175
ode differential-equations duplicate-account
April 08 2025
1  2

my answer is so different but  and i want remove this lambert in my test? can we do something for giving question exactly solve the question by model i have a book of ode which a lot of time i do some trail but have a problem? how i can remove this issue specially for bernoli and  other type?

> 

restart
> 

with(DEtools, odeadvisor)
> 

ode := diff(y(x), x)+x/y(x)+2 = 0

diff(y(x), x)+x/y(x)+2 = 0

 (1)
> 

Student:-ODEs:-Type(ode)

{}

 (2)
> 

odeadvisor(ode)

[[_homogeneous, `class A`], _rational, [_Abel, `2nd type`, `class A`]]

 (3)
> 

Student:-ODEs:-Type(ode); W := dsolve(ode); odetest(W, ode)

{}

  

y(x) = -x*(LambertW(-c__1*x)+1)/LambertW(-c__1*x)

  

0

 (4)

Download test.mw

Why the code is not apply when i try to finding p...

Asked by: salim-barzani 1605
solve parameters
April 05 2025
3  0

A lot of time i finded but i have a dubt about this why this is happen each time number of equation for finding parameter a_12 is 4 but this time is 28 which i thoght some thing must be mistake also the author of paper use  u=2(ln(f))_xx which is wronge and not satisfy but i try to find R which is strange again is not number contain parameter but is satisfy also in equation 14 i don't know each i is 2 or 1 or it can be i remain itself?

thanks for any help ?

t1.mw

How get piecewise function by substitution? ...

Asked by: salim-barzani 1605
substitution piecewise limit
April 02 2025
3  2

In this example by applying the substitution i can get half of paicewise function but how get another  half ? i am looking for B_rs as Piecewise function ?

> 

restart
> 

eij := ((-3*k[i]*(k[i]-k[j])*l[j]+beta)*l[i]^2-(2*(-3*k[j]*(k[i]-k[j])*l[j]*(1/2)+beta))*l[j]*l[i]+beta*l[j]^2)/((-3*k[i]*(k[i]+k[j])*l[j]+beta)*l[i]^2-(2*(3*k[j]*(k[i]+k[j])*l[j]*(1/2)+beta))*l[j]*l[i]+beta*l[j]^2)

((-3*k[i]*(k[i]-k[j])*l[j]+beta)*l[i]^2-2*(-(3/2)*k[j]*(k[i]-k[j])*l[j]+beta)*l[j]*l[i]+beta*l[j]^2)/((-3*k[i]*(k[i]+k[j])*l[j]+beta)*l[i]^2-2*((3/2)*k[j]*(k[i]+k[j])*l[j]+beta)*l[j]*l[i]+beta*l[j]^2)

 (1)
> 

eval(eij, k[j] = b*k[i]); series(%, k[i], 3); convert(%, polynom); eval(%, b = k[j]/k[i]); Bij := (%-1)/(k[i]*k[j])

((-3*k[i]*(-b*k[i]+k[i])*l[j]+beta)*l[i]^2-2*(-(3/2)*b*k[i]*(-b*k[i]+k[i])*l[j]+beta)*l[j]*l[i]+beta*l[j]^2)/((-3*k[i]*(b*k[i]+k[i])*l[j]+beta)*l[i]^2-2*((3/2)*b*k[i]*(b*k[i]+k[i])*l[j]+beta)*l[j]*l[i]+beta*l[j]^2)

  

series(1+((-3*(-b+1)*l[j]*l[i]^2+3*b*(-b+1)*l[j]^2*l[i]+3*(b+1)*l[j]*l[i]^2+3*b*(b+1)*l[j]^2*l[i])/(beta*l[i]^2-2*beta*l[i]*l[j]+beta*l[j]^2))*k[i]^2+O(k[i]^4),k[i],4)

  

1+(-3*(-b+1)*l[j]*l[i]^2+3*b*(-b+1)*l[j]^2*l[i]+3*(b+1)*l[j]*l[i]^2+3*b*(b+1)*l[j]^2*l[i])*k[i]^2/(beta*l[i]^2-2*beta*l[i]*l[j]+beta*l[j]^2)

  

1+(-3*(-k[j]/k[i]+1)*l[j]*l[i]^2+3*k[j]*(-k[j]/k[i]+1)*l[j]^2*l[i]/k[i]+3*(k[j]/k[i]+1)*l[j]*l[i]^2+3*k[j]*(k[j]/k[i]+1)*l[j]^2*l[i]/k[i])*k[i]^2/(beta*l[i]^2-2*beta*l[i]*l[j]+beta*l[j]^2)

  

(-3*(-k[j]/k[i]+1)*l[j]*l[i]^2+3*k[j]*(-k[j]/k[i]+1)*l[j]^2*l[i]/k[i]+3*(k[j]/k[i]+1)*l[j]*l[i]^2+3*k[j]*(k[j]/k[i]+1)*l[j]^2*l[i]/k[i])*k[i]/((beta*l[i]^2-2*beta*l[i]*l[j]+beta*l[j]^2)*k[j])

 (2)
> 

simplify((-3*(-k[j]/k[i]+1)*l[j]*l[i]^2+3*k[j]*(-k[j]/k[i]+1)*l[j]^2*l[i]/k[i]+3*(k[j]/k[i]+1)*l[j]*l[i]^2+3*k[j]*(k[j]/k[i]+1)*l[j]^2*l[i]/k[i])*k[i]/((beta*l[i]^2-2*beta*l[i]*l[j]+beta*l[j]^2)*k[j]))

6*l[j]*l[i]*(l[i]+l[j])/((l[i]-l[j])^2*beta)

 (3)


Download Lim.mw

How plot complex function including two complex f...

Asked by: salim-barzani 1605
contourplot substitution plotting
March 29 2025
4  5

why i get error in end and how i can fix this error?

> 

restart
> 

with(PDEtools)
> 

undeclare(prime, quiet); declare(u(x, y, t), quiet); declare(f(x, y, t), quiet)
> 

theta := i -> t*w[i]+y*l[i]+x:
> 

eqf := f(x, y, t) = theta(1)*theta(2)+Bij(1, 2):
> 

eqfcomplex := eval(eqf, l[2] = conjugate(l[1])):
> 

eq17 := u(x, y, t) =2*diff(f(x, y, t), x)/f(x, y, t):
> 

equ := eval(eq17, eqfcomplex):
> 

sys := map(normal, {diff(rhs(equ), x), diff(rhs(equ), y)}):
> 

nsys  := map(numer, sys):
nroot := solve(nsys, {x, y}, explicit):
> 

dsys  := map(denom, sys):
droot := solve(dsys, {x, y}, explicit):
> 

{nroot} intersect {droot}

{}

 (1)
> 

compact_ans1 := nroot[1]:
> 

__w := seq(w[i] = (-beta*l[i]^2 - b*l[i] - a), i=1..2):
> 

__Bij := (i,j) -> 12*alpha/(beta*(l[i] - l[j])^2):
> 

eval(eval(compact_ans1, {__w, Bij(1, 2) = __Bij(1, 2)}), l[1]=lambda[1]+I*lambda[2])
assuming lambda[1]::real, lambda[2]::real:
 
> 

ans1 := map(simplify, %, size): # it's up to you to use another simplification strategy
> 

eqp1 := eval(eval(ans1, l[2] = conjugate(l[1])), l[1] = lambda[1]+I*lambda[2])
> 

NULL
> 

# Do the same for the other nroot solutions
> 

eqp := {x = xp+((1/2)*beta*lambda[2]^3+I*(-beta*lambda[1]-b)*lambda[2]^2*(1/2)-((1/2)*beta*conjugate(lambda[1]+I*lambda[2])^2-(1/2)*beta*lambda[1]^2+(1/2)*b*conjugate(lambda[1]+I*lambda[2])+a)*lambda[2]+I*lambda[1]*(conjugate(lambda[1]+I*lambda[2])-lambda[1])*(beta*conjugate(lambda[1]+I*lambda[2])+beta*lambda[1]+b)*(1/2))*t/lambda[2], y = yp-(I*beta*lambda[2]^2+(2*beta*lambda[1]+b)*lambda[2]+I*((conjugate(lambda[1]+I*lambda[2])+lambda[1])*beta+b)*(conjugate(lambda[1]+I*lambda[2])-lambda[1]))*t/(2*lambda[2])}
> 

NULL
> 

vx, vy := diff(eval(x, eqp), t), diff(eval(y, eqp), t); dydx := simplify(vy/vx)
> 

eqfp := dchange(eqp, eqfcomplex, [xp, yp], params = [a, b, alpha, beta, `λ__1`, `λ__2`], simplify); eq17p := dchange(eqp, eq17, [xp, yp], params = [a, b, alpha, beta, `λ__1`, `λ__2`], simplify); eqt := simplify(eval(eq17p, eqfp))
> 

eqt1 := eval(subs({xp = x, yp = y}, eqt), l[1] = lambda[1]+I*lambda[2])
> 

with(plots); lambda[1] := .14; lambda[2] := .68; alpha := -.46; beta := 1.83; a := 1.56; b := -.19; eq := y = (-beta*conjugate(lambda[1]+I*lambda[2])^2-b*conjugate(lambda[1]+I*lambda[2])-beta*lambda[2]^2+I*(2*beta*lambda[1]+b)*lambda[2]+lambda[1]*(beta*lambda[1]+b))*(x+(2*I)*sqrt(3)*lambda[1]*sqrt(alpha/(beta*(lambda[1]+I*lambda[2]-conjugate(lambda[1]+I*lambda[2]))^2))/lambda[2])/((lambda[1]+I*lambda[2])*beta*conjugate(lambda[1]+I*lambda[2])^2+(lambda[1]+I*lambda[2])*b*conjugate(lambda[1]+I*lambda[2])-I*beta*lambda[2]^3+(-beta*lambda[1]-b)*lambda[2]^2+I*(-beta*lambda[1]^2+2*a)*lambda[2]-beta*lambda[1]^3-b*lambda[1]^2); U := proc (x, y, a, b, alpha, beta, `λ__1`, `λ__2`) options operator, arrow; rhs(eqt1) end proc; contour1 := contourplot(eval(U(x, y, a, b, alpha, beta, `λ__1`, `λ__2`), t = -50), x = -200 .. 200, y = -100 .. 100, contours = 30, color = red, grid = [100, 100], transparency = .1); contour2 := contourplot(eval(U(x, y, a, b, alpha, beta, `λ__1`, `λ__2`), t = 0), x = -200 .. 200, y = -100 .. 100, contours = 30, color = green, grid = [100, 100], transparency = .1); contour3 := contourplot(eval(U(x, y, a, b, alpha, beta, `λ__1`, `λ__2`), t = 50), x = -200 .. 200, y = -100 .. 100, contours = 30, color = blue, grid = [100, 100], transparency = .1); trajectory_plot := implicitplot(eq, x = -200 .. 200, y = -200 .. 200, color = black, thickness = 1); T := textplot([[100, 45, "t=50", color = blue], [45, -10, "t=0", color = green], [-100, -45, "t=-50", color = red]], font = [Times, Roman, 16]); display(contour1, contour2, contour3, trajectory_plot, T, labels = ["x", "y"], scaling = constrained, size = [1200, 800])

.14

  

.68

  

-.46

  

1.83

  

1.56

  

-.19

  

y = (.4755583090+0.*I)*(x+(-0.+.1517971372*I)*3^(1/2))

  

proc (x, y, a, b, alpha, beta, lambda__1, lambda__2) options operator, arrow; rhs(eqt1) end proc

  

Error, (in plot/iplot2d) invalid input: Plot:-ColorBar expects its 2nd argument, ymin, to be of type numeric, but received infinity

  

Error, (in plot/iplot2d) invalid input: Plot:-ColorBar expects its 2nd argument, ymin, to be of type numeric, but received infinity

  

Error, (in plot/iplot2d) invalid input: Plot:-ColorBar expects its 2nd argument, ymin, to be of type numeric, but received infinity

  

Error, (in plots:-display) expecting plot structure but received: contour1

  
8

Download line-plot.mw

How do I change the legend position so it doesn't ...

Asked by: nmacsai 260
label axis legend
March 29 2025
2  1

How do I change the legend position so it doesn't cover the plot label?

Before entertaining non-standard solutions to this problem such as building my own from scratch, is the a standard method for changing the position of the default legend generated along with a standard plot? I couldn't find anything about custom legend positions in the help sheets.

I can't help but wonder why this observed behaviour occurs by default? In what kind of plot would I prefer the plot labels to be obscured by the plot legend?



legend_covers_plot_labels.mw

Can I open Maple 2025 files in Maple 2024, if the ...

Asked by: nmacsai 260
March 27 2025
0  4

Can I open Maple 2025 files in Maple 2024? Further, what if the files don't use Maple 2025 features/packages? Does that change the outcome?

Inttrans invfourier(exp(-I*T*w),w,t) not calculat...

Asked by: Bachatero 45
March 27 2025
1  2

Using inttrans package fourier (Dirac(t-T),t,w) gives the correct answer

exp(-I*T*w).

Taking immediatel the inverse transform I get only

invfourier(exp(-I*T*w), w, t),

but the the expression remains unevaluated and cannot be brought to evaluate by any means that I know. Funny thing: Taking just plain Fourier Transform and Inverse in the usual integral form "int" it works flawlessly. Apparently Maple knows how to deal with distributions in this context. What am I doing wrong?

How determine interval of singularity in complex f...

Asked by: salim-barzani 1605
explore singularity plotting
March 27 2025
4  5

in a lot of my function i have a interval which is make my function singular and i don't know how remove this singularity even when i am change a lot of parameter with explore which explore option for plot is a little bit heavy for more than  7 or 8 parameter for running , and i know the shape of the graph is 2-soliton and 1-breather(zig-zag) but i have to see the shape and make my plot have a best shape  there is any idea for fixing this issue?

singular-interval.mw

what is wrong with this command?...

Asked by: awass 301
worksheet editing
March 27 2025
0  12

I tried the following procedure in a worksheet; Maple did not like it and claimed there was an error. However, I cannot even copy this to a Maple prompt; it jumps to another type of region. Any ideas? If I retype the command there is no problem with an error.

It reminds me of Maple 2 and the letter t which sometimes had to be retyped to get Maple to respond-a very strange bug which was eliminated years ago.

Can I forbid SemiAlgebraic to use RootOf and deman...

Asked by: Saalehorizontale 55
March 27 2025
0  0

Hey guys, 

 

I try to solve big systems of polynomial equations and inequalities. Therefore I use the command SemiAlgebraic. In the moment I take those result and want to go on calculating with them. Sadly it turns out, that solve has some problems with RootOf expressions. It doesnt find a solution (althoug the graph shows that there is one) and gives the warning solution may have been lost. So now I though I might just aks SemiAlgebraic to give me solutions without RootOf expressions. For example you can write {x = RootOf(_Z^2 - y)+1, 1 < y, y < 2} as {x=t+1, y=t^2,1<t<2^0.5 . This might be easier to work with for solve. 

So my question is: Is there any way I can tell SemiAlgebraic precisely in what form the solution should be? 
Since the websites are down Im not able to do a first own research on this problem. So thank you in advance. 

Regards

Felix

Why i can't find x&y C-point ?...

Asked by: salim-barzani 1605
substitution duplicate-question
March 27 2025
4  2

in some equation i don't have problem but in a lot of them this problem is come up for me and i don't know how fix this issue?

> 

restart
> 

with(PDEtools)
> 

undeclare(prime, quiet); declare(u(x, y, t), quiet); declare(f(x, y, t), quiet)

``

 (1)
> 

thetai := t*w[i]+y*l[i]+x
> 

eqw := w[i] = (-1+sqrt(-4*b*beta*l[i]-4*a*beta+1))/(2*beta)
> 

Bij := proc (i, j) options operator, arrow; -24*alpha*beta/(sqrt(1+(-4*b*l[j]-4*a)*beta)*sqrt(1+(-4*b*l[i]-4*a)*beta)-1+((2*l[i]+2*l[j])*b+4*a)*beta) end proc
> 

NULL
> 

theta1 := normal(eval(eval(thetai, eqw), i = 1)); theta2 := normal(eval(eval(thetai, eqw), i = 2))
> 

eqf := f(x, y, t) = theta1*theta2+Bij(1, 2)
> 

eqfcomplex := eval(eval(eval(eqf, l[2] = conjugate(l[1])), l[1] = lambda[1]+I*lambda[2]))
> 

eq17 := u(x, y, t) = 2*(diff(f(x, y, t), x))/f(x, y, t); equ := simplify(eval(eq17, eqfcomplex))

u(x, y, t) = 8*(-(1/2)*(-4*b*beta*conjugate(lambda[1]+I*lambda[2])-4*a*beta+1)^(1/2)*(1+((-(4*I)*lambda[2]-4*lambda[1])*b-4*a)*beta)^(1/2)-b*beta*conjugate(lambda[1]+I*lambda[2])+1/2-(b*(lambda[1]+I*lambda[2])+2*a)*beta)*((1/2)*t*(1+((-(4*I)*lambda[2]-4*lambda[1])*b-4*a)*beta)^(1/2)+(1/2)*t*(-4*b*beta*conjugate(lambda[1]+I*lambda[2])-4*a*beta+1)^(1/2)+conjugate(lambda[1]+I*lambda[2])*y*beta+((lambda[1]+I*lambda[2])*y+2*x)*beta-t)/((1+((-(4*I)*lambda[2]-4*lambda[1])*b-4*a)*beta)^(1/2)*(-(-4*b*beta*conjugate(lambda[1]+I*lambda[2])-4*a*beta+1)^(1/2)*((2*y*((lambda[1]+I*lambda[2])*y+x)*beta+t*(b*t-y))*conjugate(lambda[1]+I*lambda[2])+2*x*((lambda[1]+I*lambda[2])*y+x)*beta+((b*(lambda[1]+I*lambda[2])+2*a)*t-(lambda[1]+I*lambda[2])*y-2*x)*t)+4*(I*lambda[2]-conjugate(lambda[1]+I*lambda[2])+lambda[1])*((1/2)*conjugate(lambda[1]+I*lambda[2])*b*y*beta+(a*y-(1/2)*b*x)*beta+(1/4)*b*t-(1/4)*y)*t)-4*t*(-4*b*beta*conjugate(lambda[1]+I*lambda[2])-4*a*beta+1)^(1/2)*(conjugate(lambda[1]+I*lambda[2])*(beta*(-y*((1/2)*b*(lambda[1]+I*lambda[2])+a)+(1/2)*b*x)-(1/4)*b*t+(1/4)*y)+(((I*lambda[1]*lambda[2]+(1/2)*lambda[1]^2-(1/2)*lambda[2]^2)*b+a*(lambda[1]+I*lambda[2]))*y-(1/2)*(lambda[1]+I*lambda[2])*b*x)*beta+(1/4)*(b*t-y)*(lambda[1]+I*lambda[2]))+4*y*beta*b*conjugate(lambda[1]+I*lambda[2])^2*(-((lambda[1]+I*lambda[2])*y+x)*beta+(1/2)*t)+conjugate(lambda[1]+I*lambda[2])*(-4*beta^2*(y^2*(b*(lambda[1]^2-lambda[2]^2+(2*I)*lambda[1]*lambda[2])+2*a*(lambda[1]+I*lambda[2]))+2*x*(b*(lambda[1]+I*lambda[2])+a)*y+b*x^2)+2*beta*(-4*b*(b*(lambda[1]+I*lambda[2])+a)*t^2+2*t*(y*(b*(lambda[1]+I*lambda[2])+a)+b*x)+y*((lambda[1]+I*lambda[2])*y+x))+b*t^2-t*y)+4*(-2*((I*lambda[1]*lambda[2]+(1/2)*lambda[1]^2-(1/2)*lambda[2]^2)*b+a*(lambda[1]+I*lambda[2]))*x*y-(lambda[1]+I*lambda[2])*x^2*b-2*a*x^2+12*alpha)*beta^2+2*beta*(-4*a*(b*(lambda[1]+I*lambda[2])+a)*t^2+t*(y*(b*(lambda[1]^2-lambda[2]^2+(2*I)*lambda[1]*lambda[2])+2*a*(lambda[1]+I*lambda[2]))+2*(b*(lambda[1]+I*lambda[2])+2*a)*x)+x*((lambda[1]+I*lambda[2])*y+x))+((b*(lambda[1]+I*lambda[2])+2*a)*t-(lambda[1]+I*lambda[2])*y-2*x)*t)

 (2)
> 

ans := solve({diff(rhs(equ), x), diff(rhs(equ), y)}, {x, y}, explicit)
 

``

Download critical-point.mw

Why do I get the warning solutions may have been l...

Asked by: Saalehorizontale 55
solve implicit rootof
March 26 2025
1  12

Hey guys, 

From a former calculation I got a set of points as a implicit RootOf function for an intervall. Now I want to check, if these points are in a certain area. So i thougt I take the RootOf function, the intervall and the inequalities (which describe the target area) and use the solve command. But then I get the warning, solutions may have been lost and no solution. When you draw the implicit function you can see thats in the right area (above y=1 and below y=x). So there should be a clear anwer, giving me back the whole RootOf function in the intervall.

Download QUESTI~2.MW

Since there was an error uploading the picture here the code 

restart;
Sol := {x = RootOf(_Z^2 - y, index = real[2]) + 1, 1 < y, y < 2};
area := {1 < y, y < x};
Sol_area := solve(Sol union area);
print(Sol_area);

So why do I get this warning, the calculation seems quite easy? And is there a workaround? Or a diffrent kind of solve function? SemiAlgebraic is as far as i know only for polynomials. So I got an error as well. Since the websites are down I could start an own reasearch before. So thank you in advance. 

Regards

Felix

Using Solve for special variable leeds to missing ...

Asked by: Saalehorizontale 55
solve inequality
March 26 2025
0  2

Hey guys, 

I have a problem with the solve command. And since the websites are down, I cant help myself. I have a function x(y) and an intervall for y. This function or the set of points described by the function should now be transformed to another area of the plane. So now I can get a(x,y)=1/y and b(x,y)=x/(x+y-1) with just plugging in, then everything depeends from the y Invtervall. But I want to get the form b(a) and an intervall for a. So that I can see the function directly. So i thought I just put everything into the solve command and than ask for a solution for {a,b} and expected to get what I want. (Uploading the script here gives an error)
 

restart;
Sol := solve({a = 1/y, b = x/(x + y - 1), x = (y - 1)^2, 1 < y, y < 5/4});
Sol_ab := solve({a = 1/y, b = x/(x + y - 1), x = (y - 1)^2, 1 < y, y < 5/4}, {a, b});

#expected (or wanted) solution
#with y = 1/a and the inequalities we get 4/5 < a and a < 1 with y
#with y = we get x(a)/a and y(a) and reach b = 1/a - 1
#all together: Sol={b=1/a - 1, 4/4<a, a<1};

THe problem is, that I get an empty set which is obviously wrong. So I somehow make an error when making the variables I want concrete with adding {a,b}. What do I do wrong? Or is there a better command for what I want to achieve? In this case I can solve the problem via hand, but I have more complex tranformations and mor ecomplex functions x(y), so thats why I ask for general help with my problem. 

Thank ypu in advance

Felix

Download QUESTI~1.MW

My interactive Planck Curve works in Maple 2023, b...

Asked by: erik10 791
March 25 2025
1  7

Dear Maple users

I have an Interactive Planck Curve working great in Maple 2023, but not in Maple 2024 and 2025. Can you explain why? When I drag in the Temperature slider it fails with the window: (in plottools:-getdata) range out of bounds. 

Is it a bug or has some command changed?

File attached.

Kind regards,

Erik

Planck_Curve.mw

Converting math to text but still executing...

Asked by: PerKirkegaard64 20
March 25 2025
0  2

I have a student who has a problem when closing and opening a Maple file.

It seems as if Maple turns math fields into text, but still execute when using ! or !!!

The dark red part is written in a text field, but Maple still executes

If I try to write in a math field and executes, closes Maple and opens again, this does not happen, so it is not the file that is the problem.  The student is running 2024.2 version.

Can anyone explain the problem and how to solve it.

 

I am unable to add a comment or the file. I have tried several times, without any luck

5 6 7 8 9 10 11 Last Page 7 of 44