nm

11538 Reputation

20 Badges

13 years, 115 days

MaplePrimes Activity


These are replies submitted by nm

I do not understand the result  2*z*z'+y'

This is not valid Maple expression. Why not simply use DifferentialGeometry for total derivative?


 

restart;

f:=z^2+y;

z^2+y

DifferentialGeometry:-DGsetup([x], [z, y]):

DifferentialGeometry:-ExteriorDerivative(f);

_DG([["form", "Euc", 1], [[[2], 2*z], [[3], 1]]])

 


 

Download total_derivative.mw

looks like the search there can be improved. Try using google search. Type this exactly as is

"Created In: Maple 10" site:maplesoft.com/applications

in your browser search.  This returns many pages of hits now

 

@Preben Alsholm 

I really do not remember. Too many bugs and can't remember eveything I find :)

If you do find similar one, then feel free to delete this.  I switched to using solve instead, so I am OK now.

@C_R 

as far as I know there is nothing wrong with  applying exp() to both sides of an equation. This is always valid. So if we have an equation that looks like    ln(x)=y, we can always write this as  x=exp(y) since   exp(ln(x)) = x always.

Maple does this automatically, this mathematically is always valid
 

eq:=ln(diff(y(x),x)) = x+y(x);

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

PDEtools:-Solve(eq,diff(y(x),x))

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

expand(%);

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

 


And now it is separable ode which is solved by integration.

Download exp_op.mw

@Carl Love 

fyi, Actually it works if I do not pass the assumptions to odetest.

But first obtain the result of odetest, and then in separate step, call simplify with the same assumptions. Now it works and gives zero.

So no need to change the assumptions the way you had them. I just need to replace x0,y0 with _C1 and it still works but not with odetest, but with simplify.

interface(version);

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

restart;

ode:= y(x) = x + 3*ln(diff(y(x), x)):
sol_separable:= dsolve(ode, [separable]);
assmp:= sol-> ((x,_C1)::~real, op(indets(sol, specfunc(ln))[]) >~ 0):

#do not add assuming here.
residual:=odetest(sol_separable, ode):

#add assuming here.
simplify( residual ) assuming assmp(sol_separable);

y(x) = -3*ln(exp(-(1/3)*x)-(1/3)*c__1)

0

 

 

Download V3_IC.mw

@Carl Love 

I tried your assumptions on general solution but it does not verify it. Would you know why? I replaced your sol-> ((x,x0,y0)::~real  with sol-> ((x,_C1)::~real since now there is no x0,y0 and only _C1,x are there.

restart;

ode:= y(x) = x + 3*ln(diff(y(x), x)):
sol_separable:= dsolve(ode, [separable]);
assmp:= sol-> ((x,_C1)::~real, op(indets(sol, specfunc(ln))[]) >~ 0):
odetest(sol_separable, ode) assuming assmp(sol_separable);

y(x) = -3*ln(exp(-(1/3)*x)-(1/3)*c__1)

-3*ln(3*exp(-(1/3)*x)-c__1)-x-3*ln(-1/(c__1*exp((1/3)*x)-3))

 

 

Download V2_IC.mw

works OK for me in Maple 2024.1


 

interface(version);

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

restart;

assume(g > 0):
dsolve({diff(r(t), t$2) + g * t * r(t), r(0) = r0, D(r)(0) = r1},  r(t));

r(t) = (1/2)*((g^(1/3)*GAMMA(2/3)^2*3^(2/3)*r0+(2/3)*3^(5/6)*Pi*r1)*AiryAi(-g^(1/3)*t)+AiryBi(-g^(1/3)*t)*(r0*GAMMA(2/3)^2*g^(1/3)*3^(1/6)-(2/3)*r1*Pi*3^(1/3)))/(g^(1/3)*GAMMA(2/3))

 


 

Download works_for_me.mw

@C_R 

Both solutions are correct, but the separable one requires assumptions to obtain zero result from odetest while the dAlembert one does not. In this sense, the dAlembert solution is preferable.

ode:=y(x)-x=3*ln(diff(y(x),x));
sol_dalembert:=dsolve(ode,[dAlembert]);

y(x)-x = 3*ln(diff(y(x), x))

y(x) = x, y(x) = x+3*ln(exp(-(1/3)*x)*c__1/(-1+exp(-(1/3)*x)*c__1))

odetest([sol_dalembert][2],ode)

0

sol_separable:=dsolve(ode,[separable]);

y(x) = -3*ln(exp(-(1/3)*x)-(1/3)*c__1)

residual:=odetest(sol_separable,ode);

-3*ln(3*exp(-(1/3)*x)-c__1)-x-3*ln(exp(-(1/3)*x)/(3*exp(-(1/3)*x)-c__1))

simplify(residual) assuming exp(-x/3) - c__1/3>0, x>0

0

 


 

Download difference_in_separable_and_dalembert.mw

I guess the rule of thumb here is this:

If we can find a method to solve non-linear ode in y', as in this example,  without first having to solve for y' , then use that method.

The dAlembert method does that, it does no require isolating y' first.

While in the case of solving as separable, we have to first solve for y' to isolate it, and then solve the resulting ode as seprable. This is the main difference I see.

I solved this ode by hand using both separable and dAlembert methods. I see nothing wrong with the separable approach, this is how we solve this at school. The only difference as I said, is we have to isolate y' first. And this seems to be the cause of the trouble even though I can't still see why and what is wrong with doing that. May be someone can see the subtle mistake if any in the separable solution method.

 

@Kitonum 

Thanks, I think the Maple 2024.1 gives the correct solution, not Maple 2018, because this Clairaut  ode can only have one general solution and possibly a number of singular solutions (if they are valid ofcourse).

In the above result from Maple 2018 you show, only the first general solution must be the correct one and all the others are not. 

Singular solutions do not have a constant of integration in them. So this looks like something that was fixed in later version of Maple after 2018 or in Maple 2024.

But the strange thing I was asking about are those warning messages.

I do not yet know what causes them and why they go away when trying the command one more time and do not think I saw them before from calling dsolve. They seem to be new.  I do not have Maple 2024.0 to check.

 

@dharr 

Oh, yes ofcourse. Phase plot does not work for non-autonomous system of odes. I completely forgot that when I was looking at this. It will be nice if Maple issues a WARNING message to alert the user so they know :)

 

it looks like Maple does not have the cos(2*x) trig relations (double angle trig) in its lookup tables.You could help it by using simplify with side-relation

e:=1/2 + cos(2*x)/2;
simplify(e,{cos(2*x)=1-2*sin(x)^2}):
simplify(%);

btw, the intutive thing for the user is just to do simplify() and nothing else. This is what another software does

In Maple one has to try 10 different things to get the answer. Nothing new here. See my recent question on Limit where the same thing there. One has to try different things to find which one works.

 

 

@ecterrab 

good news. After rebooting my windows 10 PC, now the error installing V 1780 is gone! And now dsolve works

No error like before. I have no idea why restarting the PC was needed. There is no protection issues on the folder that Physics uses to install to. Never touched that. There was something causing the update to fail and rebooting fixed it. 


 

interface(version);

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

Physics:-Version(latest);

`You already have the latest "Physics Updates" installed`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1780 and is the same as the version installed in this computer, created 2024, August 4, 11:27 hours Pacific Time.`

restart;

libname;

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

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

diff(diff(y(x), x), x) = (diff(y(x), x))^3-(diff(y(x), x))^2

y(0) = 3, (D(y))(0) = 1

dsolve([ode,IC])

y(x) = ln(LambertW(-RootOf(LambertW(-_Z*exp(-1)))*exp(x-1)))+4-ln(-RootOf(LambertW(-_Z*exp(-1))))

 


 

Download divide_by_zero_maple_FIXED_2024_1_august_6_2024.mw

Next time I get this error updating Physics, I will reboot the PC and try again.  

Thanks for the help.

 

 

@C_R 

How did you manage to show the JDK version column there in task manager? I can't see how to show this on my task manager;

I do not think I have any other windows application that uses JDK, at least for GUI that I know above. Matlab does not, Mathematica does not. Scientific work place does not. notepad++ that I use does not. Word does not.

It is only Maple that I know uses Java for GUI. 

This is a brand new windows 10 installation, so I think everything is up-to-date. I took the PC to best buy and they installed windows 10 fresh since my C drive died and everything is new. But I will double check again that graphics drivers are update to date or not.

 

 

@Ronan 

Displays show the right monitor as primary (the one I get shaking on) because it says "1" on it. the  left monitor says "2". Is this right? This makes the right one the primary? And when I scroll down it has the checkmark set on it. So the right monitor is the primary. I also checked by clicking "DETECT" and it is correct.  

 

I will try swapping them by making the left monitor primary instead as you suggested and reboot windows and see if this fixes it.  

 

First 13 14 15 16 17 18 19 Last Page 15 of 90