Here is another ode solution by Maple, I am not able to get zero from odetest.
If someone could come up with a trick or method to verify this, that will be great. I tried all sort of assumptions and not able to get zero. Even coulditbe() could not give zero.
It is a Chini ode and also homogeneous, `class G`.
Solving using either method, odetest do not returns zero. There are no initial conditions. I am sure the solutions are correct, but odetest need some help or may be this requires different method to verify not using odetest?
Worksheet below.
|
| > |
ode:=diff(y(x),x) = a*x^(n/(1-n))+b*y(x)^n;
DEtools:-odeadvisor(ode,y(x));
|
![[[_homogeneous, `class G`], _Chini]](/view.aspx?sf=238925_question/c1ad2cf9ac0a85cf2d58fedc79934fc6.gif)
| > |
maple_sol:=dsolve(ode,y(x),[Chini],useInt,implicit); #explicit makes it take very long time
|
| > |
residual:=odetest(maple_sol,ode); #how to show this is zero?
|
| > |
coulditbe(residual=0) assuming real;
|
| > |
#now solving as homog. (which gives simpler looking solution
|
| > |
maple_sol:=dsolve(ode,y(x),[homogeneous]);
|
| > |
residual:=odetest(maple_sol,ode); #how to show this is zero?
|
| > |
coulditbe(residual=0) assuming real;
|
|
Download challenge_odetest_sept_1_2024.mw
