nm

8552 Reputation

19 Badges

13 years, 25 days

MaplePrimes Activity


These are questions asked by nm

Is this a  false positive, where Maple is solving an ODE which is supposed to be unsolvable?

Accoding to http://www.maplesoft.com/compare/mathematica_analysis/Comparison_Maple_Mathmatica_DEs_Kamke.pdf and considering ODE 13

Maple 18.01 does give an answer for the above ODE. I verfied the ODE from the book as well. The answer returned by Maple is very large, but it does solve it in 195 CPU seconds. Therefore the question is: Is this a false result? Or is the above document have an error in it and ODE 13 is actually solvable?

restart;
ode:=diff(y(x), x$2)-(a*y(x)^2+b*x*y(x)+c*x^2+alpha*y(x)+beta*x+gamma)^(-3/2);
sol:=dsolve(ode,y(x)) assuming a::NonZero; #I get an answer with or without this assumption. The book has the assumption
odetest(sol,ode);

btw,

 odetest(sol,ode)

gives an error as well. May be this is related to the issue or not. Not sure now.

For some reason, Maple does not use \frac{}{} for things like 1/2 that can occur in an expression

restart;
expr:=1/2;
latex(expr);

            1/2

This make the generated code not good looking. For example using this in Maple:

latex(diff(y(x), x)+y(x)*cos(x)-(1/2)*sin(2*x) = 0);

produces this latex when proccessed:

When the latex is put inside a display in the document.  This is terrible. Is there a way to tell Latex to use \frac{}{} for fractions?

 

I have few ode's which are solved by dsolve, but I am not able to get a zero from odetest(sol,ode). I tried the implicit option on the one which returns implicit solution, but I still do not get zero. I tried useInt as well.

Is there something else to do to verify the solution? My understanding is that if Maple returns a solution from dsolve, there there should be a way to get odetest() to verify the solution, but I could be wrong. Here are few examples, I have more if needed. This is Maple 18.01 on windows 7

restart;
MathematicalFunctions:-Version();
#       "C:\Program Files\Maple 18\lib\DEsAndMathematicalFunctions18.mla", `2014, July 25, 21:22 hours`
unassign(`print/ODESolStruc`):
ode1:=diff(y(x), x)+2*tan(y(x))*tan(x)-1:
ode2:=2*(diff(y(x), x))-3*y(x)^2-4*a*y(x)-b-c*exp(-2*a*x):
ode3:=(x^2+1)*(diff(y(x), x))+(y(x)^2+1)*(2*x*y(x)-1):
ode4:=x^7*(diff(y(x), x))+(2*(x^2+1))*y(x)^3+5*x^3*y(x)^2:
ode5:=(y(x)-x)*sqrt(x^2+1)*(diff(y(x), x))-a*sqrt((y(x)^2+1)^3):

sol1:=dsolve(ode1,y(x)):
sol2:=dsolve(ode2,y(x)):
sol3:=dsolve(ode3,y(x)):
sol4:=dsolve(ode4,y(x)):
sol5:=dsolve(ode5,y(x)):

odetest(sol1,ode1,implicit);  #not zero
odetest(sol2,ode2);             #not zero
odetest(sol3,ode3,implicit);  #not zero
odetest(sol4,ode4,implicit);  #not zero
odetest(sol5,ode5,implicit);  #not zero

I am going by the assumption that when Maple returns ODESolStruct as solution, then it means it could not solve the ODE. (example below)

My only complaint is that the syntax it uses for saying that the solution is ODESolStruct is not clear. I guess one has to look for & in the solution to know the result is ODEStruct.

http://www.maplesoft.com/support/help/Maple/view.aspx?path=dsolve%2fODESolStruc

Only when I convert the solution to string, then I can see the word "ODESolStruct" displayed.

My question is, how can I make maple display on the screen the word "ODESolStruct" in the solution, instead of using those "&" As that will make it more clear.

I am using worksheet on maple 18. Not document style. Here is an example:

restart;
ode:=diff(y(x),x$2)+a*exp(x)*sqrt(y(x));
sol:=dsolve(ode,y(x));


 The above was using 2D math display as default. If I use Maple notation as output I get:

----------------------------------------

restart:
ode:=diff(y(x),x$2)+a*exp(x)*sqrt(y(x)):
sol:=dsolve(ode,y(x));
sol := y(x) = `&where`(_a/exp(-2*(Int(_b(_a), _a))-2*_C1), [{diff(_b(_a), _a).......

-------------------------------------------

But now

convert(sol,string);
"y(x) = ODESolStruc(_a/exp(-2*Int(_b(_a),_a)-2*_C1),[{diff(_b(_a\ .............."

You can see now that the solution is ODESolStruct, but it is much more clear than the default solution above. But only when looking at the solution as string do I get it to show the word "ODESolStruct". 

Since odetest does not return zero, then maple did not solve it:

odetest(sol,ode);

btw, Compare the above to when Maple returns "DESol" structutre. In this case, it does now display on the screen the word "DESol":

restart;
ode:=diff(y(x), x, x)-y(x)*(a^2*x^(2*n)-1);
dsolve(ode,y(x));

Again, my question is:  Could I configure Maple to display in worksheet the solution using explicit ODESolStruct words instead of using "&" there to indicate more clearly the solution.

 

Maple 18.01, windows

restart;
ode:=2*a^2*y(x)-2*y(x)^3+3*a*(diff(y(x), x))+diff(y(x), x$2)=0;
dsolve(ode,y(x));

           returns y(x)=0

 

So does

ode:=2*a^2*y(x)-20*y(x)^3+3*a*(diff(y(x), x))+diff(y(x), x$2)=0;

ode:=2*a^2*y(x)-200*y(x)^3+3*a*(diff(y(x), x))+diff(y(x), x$2)=0;

ode:=2*a^2*y(x)-2000*y(x)^3+3*a*(diff(y(x), x))+diff(y(x), x$2)=0;

etc...

Is this a bug?

 

 

 

First 153 154 155 156 157 158 159 Last Page 155 of 164