MaplePrimes Questions

Hi,

I’m having fun animating a beautiful geometric shape starting from a few trigonometric functions. I’m wondering if there’s a way to link each curve in the animation to its name.

Animation_Trigo.mw

How to get amultiple line plots in one graph for differnt values of E1 := .1,0.2,0.3,and 0.4 with differnt color line red,blue,black and green and how to get numerical values for the all E1 value in one matrix form 

at E1=0.1 value of diff(g(x), x, x), diff(f(x), x, x),diff(f(x), x) and diff(g(x),  x).

AGM_method_single_line_plot.mw

I recently upgraded from Maple 23 to Maple 24. While many display issues have been resolved, I’ve encountered a new real problem: when entering operations or a factorial in a denominator or exponent, the cursor unexpectedly jumps to the inline position. This forces me to manually reposition the cursor using backspace plus left-arrow keys, and forgetting to do so can lead to errors. I have a perpetual Maple license through my university and haven’t purchased maintenance this time. Is there any way to fix or work around this cursor-jumping issue in Maple 24 without purchasing a new license?

The_Bohrs_Model_-_MaplePrimes.mw

Look at the equation (11) in the Maple's document. I would like to force Maple to let the variable "r" inside the squared root so to get the equation (12). Any idea of doing that?  Thank you in advance for your help.

Already  by Help of my favorite Dr.David he did find the thus three step for non schrodinger equation but in here i got some issue of coding which is so different from before, is about transformation of pdes to two parts od real and imaginary part and then substitution our function the functions is clear but combining them and findinf leading exponent and resonance point and finding function in step 3 in different and jsut the function is different with eperate the real and imaginary part for finding step one ...

note:=q=u*exp(#) then |q|=u

schrodinger-test.mw

paper-1

paper-2

(2500iw/(1+5iw) )+(200iw/1-10iw)+5 rationalize and simplify

Hello All

I have the system of 4 differential equations with retarted delay (t-tau) and advanced argument (t+tau). The problem is also considered as a two-value boundary problem since some variables have initial conditions and some have terminal ones.

Here is a model:

Data:

A := 10;
B := 5;
C := 2;
tau = 0.5;
T := 10;
 z := 5;
 N := 20;
 

System:

sys := diff(K(t), t) = A*Q(t - tau) + K(t), diff(Q(t), t) = l_Q(t)*z - N, diff(l_K(t), t) = B*Q(t - tau) + l_K(t), diff(l_Q(t), t) = l_K(t + tau)*K(t + tau) + C;

Initial and terminal conditions:
 
inc := K(0) = 50, Q(0) = 200, l_K(T) = -20, l_Q(T) = 0;
   

I will appreciate any direction/suggestion on how this system can be solved. 

Dmitry

I was surprised that Maple can't solve this first order ode which is exact ode.

I solved by hand and Maple says my solution is correct.

Any one can find why Maple failed to solve this and if older versions can solve it? Also tried implicit option, but that did not help.

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

restart;

ode:=diff(y(x),x) = (2*sin(2*x)-tan(y(x)))/x/sec(y(x))^2;

diff(y(x), x) = (2*sin(2*x)-tan(y(x)))/(x*sec(y(x))^2)

sol:=dsolve(ode);

mysol:=cos(2*x)+x*tan(y(x))=c__1;

cos(2*x)+x*tan(y(x)) = c__1

odetest(mysol,ode);

0

 

 

Download maple_solving_exact_ode_august_25_2025.mw

In the attached file, the trigonometric term (2, term) is transformed into a term (3, term1) consisting of radicals. Is there a Maple procedure that can be used to reverse this process? Given an algebraic term (e.g., consisting of radicals, powers, etc.), under what conditions can it be transformed into a trigonometric form (not a Fourier series) in the sense of (3) according to (2)?test.mw

 interface(version);

`Standard Worksheet Interface, Maple 2024.2, Windows 11, October 29 2024 Build ID 1872373`

(1)

restart

term := 2*cos(5*arcsin((1/2)*x))

2*cos(5*arcsin((1/2)*x))

(2)

term1 := expand(term)

-3*(-x^2+4)^(1/2)*x^2+(-x^2+4)^(1/2)+(-x^2+4)^(1/2)*x^4

(3)

convert(term1, trig)

-3*(-x^2+4)^(1/2)*x^2+(-x^2+4)^(1/2)+(-x^2+4)^(1/2)*x^4

(4)

simplify(term1, trig)

(-x^2+4)^(1/2)*(x^4-3*x^2+1)

(5)

solve(term1 = sqrt(2), x)

(1/2)*(8-2*(10+2*5^(1/2))^(1/2))^(1/2), (1/2)*(8+2*(10-2*5^(1/2))^(1/2))^(1/2), (1/2)*(8+2*(10+2*5^(1/2))^(1/2))^(1/2), -(1/2)*(8-2*(10+2*5^(1/2))^(1/2))^(1/2), -(1/2)*(8+2*(10-2*5^(1/2))^(1/2))^(1/2), -(1/2)*(8+2*(10+2*5^(1/2))^(1/2))^(1/2)

(6)

evalf(solve(term1 = sqrt(2), x))

.3128689302, 1.782013048, 1.975376681, -.3128689302, -1.782013048, -1.975376681

(7)

plot(term, x = -2.5 .. 2.5)

 

plot(term1, x = -2.5 .. 2.5)

 

NULL

Download test.mw

As we can see, RealDomain:-solve gives an incorrect solution to the following system: 

restart;

sys := `~`[diff](sqrt(2*a^2-8*a+10)+sqrt(b^2-6*b+10)+sqrt(2*a^2-2*a*b+b^2), [a, b]):

RealDomain:-solve(`~`[`=`](sys, 0), {a, b})

{a = 5/3, b = 5/2}, {a = a, b = 2*a/(a-1)}

(1)

plot(eval(sys, {max(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2)), min(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2))}[-1]), a = -infinity .. infinity)

 

extrema(sqrt(2*a^2-8*a+10)+sqrt(b^2-6*b+10)+sqrt(2*a^2-2*a*b+b^2), {}, {a, b})

{max(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2)), min(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2))}

(2)

Download solve_returns_an_unsatisfiable_real_solution.mw

This appears to be a bug; is it possible to fix it? 
Text: 

sys := diff~(sqrt(2*a^2 - 8*a + 10) + sqrt(b^2 - 6*b + 10) + sqrt(2*a^2 - 2*a*b + b^2), [a, b]):
RealDomain:-solve(sys =~ 0, {a, b});

THis is problem from textbook. Maple do not give solution. 

But when asked for implicit solution, it gives one.  Should it not have done this automatically?

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

ode:=y(x)*diff(y(x),x) = a;
ic:=y(0) = b;
sol:=dsolve([ode,ic]);

y(x)*(diff(y(x), x)) = a

y(0) = b

sol:=dsolve([ode,ic],'implicit')

-2*a*x+y(x)^2-b^2 = 0

 

 

Download why_no_solution_maple_2025_1.mw

We see now there are two solutions for y(x), since quadratic.

So why dsolve do not solve this and at least give implicit solution automatically? Should this be reported as defect?

is been a while i work on a test still i am study and there is a lot paper remain and is so important in PDEs, a lot paper explain in 2003 untill know and there is other way to find it too but i choose a easy one but is 2025 paper  which is explanation is so beeter than other paper, also some people write a package for take out this test with a second but maybe is not work for all i  search for that  but i didn't find it i will ask the question how we can find thus as shown in graph i did my train but need a little help while i am collect more information and style of solving 

Download p1.mw

Hello everyone,
I hope this message finds you well. I am trying to plot a function f(x, y) and overlay its contour on a quarter ellipse using Maple 2015. However, I’ve encountered some difficulties and have not been successful so far. I would greatly appreciate any assistance in resolving this issue. Thank you!

Plotting in 2D

restart:with(plots):  aa := 4: bb := 2:  
f := -((x^(2))/(aa^(2))+(y^(2))/(bb^(2))-1)*((aa^(2)*bb^(2))/(aa^(2)+bb^(2))):  
plot3d(f,x = 0 .. aa/(2),y = 0 .. bb/(2),region = (x, y) -> ((2 x)/(aa))^(2) + ((2 y)/(bb))^(2)<= 1,axes = boxed,style = patchcontour, grid = [50, 50],orientation = [-120, 45],shading = zhue,title = "f(x,y) over quarter ellipse domain");

Contour plotting
xrange := 0 .. aa/(2): yrange := 0 .. bb/(2):  
nx := 100:   ny := 100:  
dx := (rhs(xrange) - lhs(xrange))/(nx-1):dy := (rhs(yrange) - lhs(yrange))/(ny-1):  
Z := Matrix(nx, ny, (i, j) -> local x, y, inside;x := lhs(xrange) + (i-1)*dx;y := lhs(yrange) + (j-1)*dy;inside := (((2 x)/aa)^2 + ((2 y)/bb)^2 <= 1);if inside then f(x, y) else NULL end if):  
contourplot(Z, xrange, yrange,contours = 15, filled = true, coloring = [blue, green, yellow, red], axes = boxed, title = "Contour plot over quarter ellipse", grid = [nx, ny]);  

Hi,
I would like to solve a system of 2 PDE with this code, but it does not work.
lin_eqs := {diff(h(x, t), t) + H*diff(u(x, t), x) = 0, diff(u(x, t), t) + g*diff(h(x, t), x) = 0};
sol := pdsolve(lin_eqs);
print('sol', sol);

Can you help me?Thanks
Best regards

Any idea why Maple simplifies 1+sin(x)^2 to 2-cos(x)^2?  Leaf count is larger also.

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

e1:=1+sin(x)^2;

1+sin(x)^2

e2:=simplify(e1)

-cos(x)^2+2

MmaTranslator:-Mma:-LeafCount(e1)

6

MmaTranslator:-Mma:-LeafCount(e2)

8

 

 

Download strange_simplification_august_20_2025.mw

1 2 3 4 5 6 7 Last Page 1 of 2428