nm

11458 Reputation

20 Badges

13 years, 76 days

MaplePrimes Activity


These are questions asked by nm

I solved this linear constant coefficients ode by hand, and wanted to use Maple to check my solution.

dsolve works with no problem. I wanted to check the particular solution on its own, so called DEtools:-particularsol(ode,y(x)); but was surprised it hangs.

Why would DEtools:-particularsol(ode,y(x)); hang, when dsolve gives the complete solution (including particular solutin ofcourse) right away?  It is 7th order ODE. But nothing special. Linear and constant coefficients. Since dsolve solves it right away, I expected DEtools:-particularsol to have no problem as well.

restart;
ode:=diff(y(x),x$7)-2*diff(y(x),x$6)+9*diff(y(x),x$5)-16*diff(y(x),x$4)+24*diff(y(x),x$3)-32*diff(y(x),x$2)+16*diff(y(x),x)=exp(2*x)+x*sin(x)+x^2;

dsolve(ode);  #no problem.

#this hangs. Why?
DEtools:-particularsol(ode,y(x));


 

interface(version);

`Standard Worksheet Interface, Maple 2020.1, Windows 10, July 30 2020 Build ID 1482634`

Physics:-Version()

`The "Physics Updates" version in the MapleCloud is 851. The version installed in this computer is 847 created 2020, October 17, 17:3 hours Pacific Time, found in the directory C:\Users\me\maple\toolbox\2020\Physics Updates\lib\`

restart;
ode:=diff(y(x),x$7)-2*diff(y(x),x$6)+9*diff(y(x),x$5)-16*diff(y(x),x$4)+24*diff(y(x),x$3)-32*diff(y(x),x$2)+16*diff(y(x),x)=exp(2*x)+x*sin(x)+x^2;

diff(diff(diff(diff(diff(diff(diff(y(x), x), x), x), x), x), x), x)-2*(diff(diff(diff(diff(diff(diff(y(x), x), x), x), x), x), x))+9*(diff(diff(diff(diff(diff(y(x), x), x), x), x), x))-16*(diff(diff(diff(diff(y(x), x), x), x), x))+24*(diff(diff(diff(y(x), x), x), x))-32*(diff(diff(y(x), x), x))+16*(diff(y(x), x)) = exp(2*x)+x*sin(x)+x^2

#this works
sol:=dsolve(ode)

y(x) = (5/16)*x+(1/128)*exp(2*x)+(1/8)*x^2-(1/4)*exp(-(2*I)*x)*_C3-(1/4)*exp((2*I)*x)*_C3-exp(x)*_C4+(1/8)*exp(-(2*I)*x)*_C5+(1/8)*exp((2*I)*x)*_C5+((1/36)*I)*exp(-I*x)-((1/36)*I)*exp(I*x)-((67/57600)*I)*exp(-(2*I)*x)+((67/57600)*I)*exp((2*I)*x)+((1/4)*I)*exp(-(2*I)*x)*_C5*x-((1/4)*I)*exp((2*I)*x)*_C5*x+(67/28800)*x*exp((2*I)*x)+(67/28800)*x*exp(-(2*I)*x)-(1/4)*exp(-(2*I)*x)*_C6*x-(1/4)*exp((2*I)*x)*_C6*x-((1/36)*I)*exp(I*x)*x+((1/36)*I)*exp(-I*x)*x+((67/21600)*I)*exp(-(2*I)*x)*x+((1/4)*I)*exp(-(2*I)*x)*_C2+((1/8)*I)*exp(-(2*I)*x)*_C6-((67/21600)*I)*exp((2*I)*x)*x-((1/4)*I)*exp((2*I)*x)*_C2-((1/8)*I)*exp((2*I)*x)*_C6+(67/43200)*exp((2*I)*x)+(67/43200)*exp(-(2*I)*x)+(1/54)*exp(-I*x)+(1/54)*exp(I*x)+_C4*exp(x)*x+(1/48)*x^3+_C7+_C1*exp(x)

#this hangs. Why?
DEtools:-particularsol(ode,y(x));

 


 

Download particular_sol.mw

I was saying to myself today, that one nice thing in Maple, is that it has good type system so one can check type of arguments aginst wrong types being used in the call.

So I was surprised when I called ArrayTools:-AllNonZero  using list as argument, where this function is supposed to accept only Matrix,Vector or Array. And it worked.  But it gave wrong answer at the same time.

Am I doing something wrong here? Should not have this call failed to go through?

restart;
interface(warnlevel=4);
kernelopts('assertlevel'=2):

A:=[0,0,0];
ArrayTools:-AllNonZero(A)

Maple replied 

                       true

But 

A:=Array([0,0,0]);
ArrayTools:-AllNonZero(A)

Now Maple gives the expected result  false

The question is why Maple did not detect the wrong type? Should it have detected wrong type?

And why it even gave wrong answer (but this is not as important, since the call should not be allowed in first place).

btw, this applied to all the other API's listed

 

 

I found a new strange thing with odetest that I have not seen before.

Would you say that these two solutions are the same or not?

They look the same to me.  One just have all the terms on one side, that is all.

Why would Maple odetest verify the first one but not the second? I did not know that all terms have to be on one side before. Is this documented somewhere?

Please see worksheet below.

restart;

interface(version);

`Standard Worksheet Interface, Maple 2020.1, Windows 10, July 30 2020 Build ID 1482634`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 847 and is the same as the version installed in this computer, created 2020, October 17, 17:3 hours Pacific Time.`

ode:=diff(y(x), x) - 2*y(x) = 2*sqrt(y(x));
ic:=[y(0)=1];
maple_sol:=dsolve([ode,op(ic)],implicit);
odetest(maple_sol,[ode,op(ic)]) assuming x>0

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

[y(0) = 1]

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

[0, 0]

#now move some terms to one side, and odetest no longer verifies it
my_sol:=y(x)^(1/2) = -1+2*exp(x);
odetest(my_sol,[ode,op(ic)]) assuming x>0

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

[diff(y(x), x)-2*y(x)+2-4*exp(x), 0]

#the solution is to put all terms on one side
my_sol:=y(x)^(1/2) = -1+2*exp(x);
odetest(rhs(my_sol)-lhs(my_sol)=0,[ode,op(ic)]) assuming x>0

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

[0, 0]

 


 

Download odetest_issue_oct_18_2020.mw

 

Using solve on this example:

restart;
eq:=x^3 - 3*x^2 + 3*x - 1=0;
solve(eq,x);

gives solution x=1, with multiplicty 3

              1, 1, 1

When using PDEtools:-Solve

restart;
eq:=x^3 - 3*x^2 + 3*x - 1=0;
solve(eq,x);
PDEtools:-Solve(eq,x);

it gives

         x = 1

How can one make it show x=1, x=1, x=1 ?

Found that only when using PDEtools:-Solve(eq,x,numeric) it gives 

           x = 1., x = 1., x = 1.

How to make it do the same without using numeric?

Maple 2020.1.1

 

What is the reason that error accepts its arguments without being inside paranthesis? This is different from all other Maple functions I've seen. For example, both these below work the same

f := proc (x) 
     if x<0 then 
        error "invalid x: %1", x; 
     else 
        x^(1/2); 
     end if 
end proc:

And

g := proc (x) 
     if x<0 then 
        error("invalid x: %1", x); 
     else 
        x^(1/2)
     end if 
end proc:

One can't do for other Maple functions. For example sin x  gives an error. It must be sin(x)

Looked at help page for error and did not see something about this difference.

First 111 112 113 114 115 116 117 Last Page 113 of 201