Maple 2024 Questions and Posts

These are Posts and Questions associated with the product, Maple 2024

I sort of remember there is a special syntax for setting initial condition for an ode derivative as    y'(a)=b, where and are symbols. I forgot what it is, everything I try gives error. (I thought eval was used to work for this, but can't get it to work now).

Any one knows how to set this IC?   For an example, given this ode y''(x)+y'(x)+y(x)=0 I want to solve it with the IC as   y'(a)=b where has no numerical value. Just a symbol.

For an example, using another software, it is done as follows

ClearAll[x,y,a,b];
ode=y''[x]+y'[x]+y[x]==0;
DSolve[{ode,y'[a]==b},y[x],x]

Below is my attempts in Maple. I tried eval, subs, and the normal D(y)(a)=b but none of these works when is symbol. 

I looked at help and all examples I saw use constants, as in D(y)(0)=value. 

Iam sure this is possible to do in Maple, but I forgot how. Below is worksheet.
 

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

restart;

ode:=diff(y(x),x$2)+diff(y(x),x)+y(x)=0;
IC1:=eval(diff(y(x),x),x=a)=b;
IC2:=D(y)(a)=b;
IC3:=subs(x=a,D(y)(x)=b);

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

diff(y(a), a) = b

(D(y))(a) = b

(D(y))(a) = b

dsolve([ode,IC1]);

Error, (in dsolve) found differentiated functions with same name but depending on different arguments in the given DE system: {y(a), y(x)}

dsolve([ode,IC2]);

Error, (in dsolve) found differentiated functions with same name but depending on different arguments in the given DE system: {y(a), y(x)}

dsolve([ode,IC3]);

Error, (in dsolve) found differentiated functions with same name but depending on different arguments in the given DE system: {y(a), y(x)}

#eval works if the point is actual value as in
IC1:=eval(diff(y(x),x),x=0)=b;

eval(diff(y(x), x), {x = 0}) = b

#but I want x=a in the above.


Download how_to_set_IC_point_to_symbol.mw

There is no issue for dirichlet initial condition, where  y(a)=b works. It is the  neumann one which I can't figure its syntax. So this works OK

ode:=diff(y(x),x$2)+diff(y(x),x)+y(x)=0;
IC1:=y(a)=b;
dsolve([ode,IC1]);

Maple help pages keep getting worst with each release.

I want all input to be displayed using Maple 1D notation, so I can copy the example to my worksheet since I only use worksheet and not document (2D) mode.

So even though the first thing I do when I open help it to turn off the 

                     view->Display examples with 2D

So it is no longer checked, I still see many pages using 2D math for input. 

Here is one example ?D  page

If I copy one such input to my worksheet now it looks like this

eveything in 2D becomes ?? when I copy it.

So one can only look but not copy?

Is there any other option to make sure, really make sure, all examples have 1D as input?

The problem is that it is all not consistent. Some examples have a mix of 2D and 1D as the above page. Some are in 2D and some are in 1D.

And this is all on the same help page!

Does no one inside Maplesoft even look at their own help pages?
 

At school the teacher always said that if we have second order ode and only one initial conditions (say y'(0)=0 or y(0)=0) then the solution should have one constant of integration in it.

And if we have no initial conditions, then the solution should have 2 constants of integrations in it.

And if we have two initial conditions, then the solution should have zero constants of integrations in it.

In this example, Maple is given second order ode with one IC. But the solution it gives when asked to solve it explicit, has no constant of integration in it at all. 

When asked to solve it using implicit, then the constant of integration shows up. 

Both solutions actually verify to be fully correct using odetest. So it looks like the solution as explicit is particular solution and not a general solution.

Why is that? Why it did not give general solution when asked to solve the ode as explicit?

restart;

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1793 and is the same as the version installed in this computer, created 2024, August 25, 9:6 hours Pacific Time.`

libname;

"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"

restart;

ode:=y(x)*diff(y(x),x$2)+diff(y(x),x)^2+1=0;
IC:=D(y)(0)=1;

y(x)*(diff(diff(y(x), x), x))+(diff(y(x), x))^2+1 = 0

(D(y))(0) = 1

sol1:=dsolve([ode,IC],explicit);
 

y(x) = (-x^2+x*2^(1/2)+1/2)^(1/2)

sol2:=dsolve([ode,IC],implicit);

-(1/2)*y(x)^2+x*y(0)-(1/2)*x^2+c__2 = 0

odetest(sol1,[ode,IC]);
odetest(sol2,[ode,IC]);

[0, 0]

[0, 0]

 

 

Download why_missing_constant_of_integration_august_25_2024.mw

The second result is definitely wrong, but why is this? 

MmaTranslator:-FromMma("Subscript[x, 0] -> 2 (y - z) (y - z)");
                       MmaTranslator:-FromMma(

MmaTranslator:-FromMma("Subscript[x, 0] -> 2 (y - z)^2"); # should return: x[0] = 2*(y - z)^2
                      MmaTranslator:-FromMma(

The result is still wrong even if I type an explicit multiplication sign (for instance, "Subscript[x, 0] -> 2*(y - z)^2"). 

On a Windows 10 installation with

I get after system restart and Maple launch

followed by

What could be the cause that the update is not downloading and how to fix it?

I am trying to understand under which conditions submenus disappear.

Now I have seen that the context panel is empty when submenues disappear.

Has anybody noticed the same and are there other instances when the context panel is emtpy?

Both effects could have the same rootcause.

Here is an example:
A -> no submenu

B -> Empty context panel

C -> Menu appears inside workspace (first time I see that)

 

Update

C seems to happen when the Maple window is not maximized to the full screen. -> then the context pannel becomes empty (for some reason) -> then if the Maple window is maximized to full screen we see C. Opeing new documents do not show this when toggling between and

Hi everyone,

I’m encountering an issue with Maple 2024 on my MacBook Air M1 (2020) running macOS Sonoma 14.6.1. Every time I try to log in to Maple Cloud from within the application, the entire program freezes and becomes unresponsive. I’ve tried restarting the application and my MacBook, but the problem persists.

Is anyone else experiencing similar issues with Maple 2024 on an M1 Mac running macOS Sonoma? If so, have you found any solutions or workarounds? I’ve checked for updates and even reinstalled the software, but the problem still exists.

Any help or advice would be greatly appreciated!

Thanks in advance!

This problem began about 3 nights ago.

When I load the subpackage UHG to RationalTrigonometry no warning.

When I rerun the worksheet warnings appear. If I put in arestart and reload the packages no warning. The problem only occurs with the sub packageUHG. The RationalTrigonometry package also has a few aliases.

Screen shot of one rerun  then added restart.

I've modified my call to coulditbe() to always use assuming real to prevent false positives.

Now I find it gives internal error Error, (in type/complex) too many levels of recursion which can not even be trapped when the thing I am checking happened to have complex I in it. (This is done in code).

So I modied the code to only do this check if there is no complex in the expression being checked.

But I think Maple should not generate such error in first place. This indicates some internal problem in Maple, correct?

Here is worksheet. (ps. updated, wrong worksheet for somereason was uploaded).

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1792 and is the same as the version installed in this computer, created 2024, August 22, 12:6 hours Pacific Time.`

libname;

"C:\Users\Owner\maple\toolbox\2024\Physics Updates\lib", "C:\Program Files\Maple 2024\lib"

the_residual:=-1/4*I*x^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)*6^(1/2)-1/4*x^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)*2^(1/2)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)-x^(1/2);
try
    coulditbe(the_residual=0) assuming real;
catch:
   print("error");
end try;

-((1/4)*I)*x^(1/2)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)*6^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)-(1/4)*x^(1/2)*2^(1/2)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)-x^(1/2)

Error, (in type/complex) too many levels of recursion

the_residual:=-1/4*I*x^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)*6^(1/2)-1/4*x^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)*2^(1/2)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)-x^(1/2);
try
    if not has(the_residual,I) then
       coulditbe(the_residual=0) assuming real;
    else
       print("by passing, since residual is complex");
    fi;
catch:
   print("error");
end try;

-((1/4)*I)*x^(1/2)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)*6^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)-(1/4)*x^(1/2)*2^(1/2)*(-(8*x^(3/2)+12*c__1)^(2/3)-I*(8*x^(3/2)+12*c__1)^(2/3)*3^(1/2))^(1/2)/(8*x^(3/2)+12*c__1)^(1/3)-x^(1/2)

"by passing, since residual is complex"

 


 

Download internal_error_too_many_Levels_from_coulditbe_august_23_2024.mw

Hello,

Are there any tools in Maple that allow you to find real roots of complex integral equations?

I have attached examples of such equations.

Unfortunately, I can't solve them.

Best regards.

Test.mw

Studying the attached documentDuffing Equation from the University of Colorado, I tried to plot the Duffing equation in three different ways. My worksheet is attached.
But I didn't get the desired result as you can see below. Any suggestion on how to fix my problem?

 

Download Example.mw

 

This ode has solution when solving for the IC by hand. But Maple gives new error I did not see before when asking it to solve the ode with IC. 

No error if asked to solve the ODE without the IC.

Is this new error in dsolve? I have no earlier Maple versions to check. Below is worksheet. Solving for the IC by hand gives the solution which odetest verified.

restart;

interface(version);

`Standard Worksheet Interface, Maple 2024.1, Windows 10, June 25 2024 Build ID 1835466`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1789 and is the same as the version installed in this computer, created 2024, August 10, 8:50 hours Pacific Time.`

ode := diff(y(x), x) = sqrt(2)*sqrt(-(y(x) - 6)*(2*y(x) - 3))/(-y(x) + 6);
IC:=y(0)=3;

diff(y(x), x) = 2^(1/2)*(-(y(x)-6)*(2*y(x)-3))^(1/2)/(-y(x)+6)

y(0) = 3

sol:=dsolve(ode);

x-(1/4)*(-4*y(x)^2+30*y(x)-36)^(1/2)-(9/8)*arcsin((4/9)*y(x)-5/3)+c__1 = 0

sol_with_IC:=dsolve([ode,IC]); #why this error??

Error, (in dsolve) when calling 'series/RootOf'. Received: 'unable to compute series'

#lets solve for the IC by hand:
eq:=eval(sol,[y(x)=3,x=0])

-(1/4)*18^(1/2)+(9/8)*arcsin(1/3)+c__1 = 0

C_sol:=PDEtools:-Solve(eq,_C1);

c__1 = (3/4)*2^(1/2)-(9/8)*arcsin(1/3)

my_new_sol:=eval(sol,C_sol)

x-(1/4)*(-4*y(x)^2+30*y(x)-36)^(1/2)-(9/8)*arcsin((4/9)*y(x)-5/3)+(3/4)*2^(1/2)-(9/8)*arcsin(1/3) = 0

odetest(my_new_sol,[ode,IC])

[0, 0]

 


 

Download internal_error_when_unable_to_find_solution_august_22_2024.mw

Dear Expert Users,

I am still struggling with my transition from Mathcad 15 to Maple. Would someone be so kind to look at the attached worksheet and how I can obtain the desired result? Thanks in advance. The data used in the example come from TestFleck.xlsx

QuestionPrimes.mw TestFleck.xlsx

*****************************

EDIT: Thanks for the help i used the restore backup feature. I made my anti virus not react to maple since i thought that might work, so far no more problems. Else this is a fix if anybody else has this problem.

******************************

Hi, so my maple freezes during simple tasks, such as trying to type something, or somtimes copying and inerting the results from a previous equation, the feature where you mark something and drag the answer to somewhere else also does this. When this happens, maple completely freezes and i'm forced to close the app with taksmanager, where i lose everything i have worked on, sometimes for the past hour because i didn't save every minute which is basically not something you can do.

Is there any fix for this or anything i can do to prevent this, cause it really sucks. 

I'm on a windows 11 with maple 2024

I am preparing some examples to start the year in the differential equations lesson. And I was wondering if there is any other way to represent the equilibrium points than the following which I found and then plotted using pointplot.

restart

NULL

alpha := 0.4e-1; beta := 0.8e-2; sys := {diff(g(t), t) = -.25*g(t)+beta*p(t)*g(t), diff(p(t), t) = .7*p(t)*(1-(1/100)*p(t))-alpha*p(t)*g(t)}; ics := {g(0) = 5, p(0) = 20}
 

sol := dsolve(`union`(sys, ics), {g(t), p(t)}, numeric, range = 0 .. 50); with(plots); plots:-display(plots:-odeplot(sol, [t, p(t)], 0 .. 50, color = blue, legend = ["p(t)"]), plots:-odeplot(sol, [t, g(t)], 0 .. 50, color = red, legend = ["gt)"]), labels = ["t", "Population"], title = "Population Dynamics")

 

eqns := {.7*p(t)*(1-(1/100)*p(t))-alpha*p(t)*g(t) = 0, -.25*g(t)+beta*p(t)*g(t) = 0}; eq_points := solve(eqns, {g(t), p(t)}); eq_points

{g(t) = 0., p(t) = 0.}, {g(t) = 0., p(t) = 100.}, {g(t) = 12.03125000, p(t) = 31.25000000}

(1)

eqpoints := pointplot([[0, 0], [100, 0], [31.25000000, 12.03125000]], color = [red], symbol = diamond, symbolsize = 15)

``
phaseplot := odeplot(sol, [p(t), g(t)], 0 .. 35, color = green, thickness = 2)

with(DEtools); vectorfield_plot := dfieldplot([diff(p(t), t) = .7*p(t)*(1-(1/100)*p(t))-alpha*p(t)*g(t), diff(g(t), t) = -.25*g(t)+beta*p(t)*g(t)], [p(t), g(t)], t = 0 .. 35, p = 0 .. 150, g = 0 .. 30, arrows = small, color = blue, axes = boxed)

 

 

display([vectorfield_plot, phaseplot, eqpoints])

 

NULL


 

Download aquarium.mw

First 26 27 28 29 30 31 32 Last Page 28 of 43