What is this?
how do i go about finding all the polynomials of degree 5 or less that satisfy this equation?
i found the two obvisous families of the solution to be y=a and y=bx, where a and b are constants.....
polynomial solution of degree 5 or less:
d:=diff(y(x),x)*diff(y(x),x$3)-diff(y(x),x$2); yp:=add(a||i*x^i,i=0..5); eval(d,y(x)=yp): collect(%,x): {coeffs(%,x)}; solve(%,{a0,a1,a2,a3,a4,a5}); {a0 = a0, a1 = a1, a2 = 0, a3 = 0, a4 = 0, a5 = 0}
A first integral is y''=ln(C1*y').
y=a+bx are the only polynomial solutions.
Look at it that way. If y is a polynomial of degree <3, then y'''=0, so the lhs is 0. So the r.h.s. should be also 0. That gives the solution y=a+bx.
If y is a polynomial of degree n>=3, the l.h.s is a polynomial of degree
(n-1)+(n-3) = 2n-4,
and the r.h.s. is a polynomial of degree n-2. That gives an equation
2n-4 = n-2,
from which n=2 that contradicts to n being >=3.
Alec
solve for y'y'''=y'
If that's really what the equation is, it's trivial: the difference of the two sides of the equation factors as
y' (y''' - 1)
so there are two cases: y' = 0 and y''' = 1.
Maple has a command for this. It gives the single solution that Alec found above
d:=diff(y(x),x)*diff(y(x),x$3)-diff(y(x),x$2):
DEtools[polysols](d,y(x));
gives
[y(x) = _C3+_C4*x]
no other
polynomial solution of degree 5 or less:
d:=diff(y(x),x)*diff(y(x),x$3)-diff(y(x),x$2); yp:=add(a||i*x^i,i=0..5); eval(d,y(x)=yp): collect(%,x): {coeffs(%,x)}; solve(%,{a0,a1,a2,a3,a4,a5}); {a0 = a0, a1 = a1, a2 = 0, a3 = 0, a4 = 0, a5 = 0}A first integral is y''=ln(C1*y').
Also
y=a+bx are the only polynomial solutions.
Look at it that way. If y is a polynomial of degree <3, then y'''=0, so the lhs is 0. So the r.h.s. should be also 0. That gives the solution y=a+bx.
If y is a polynomial of degree n>=3, the l.h.s is a polynomial of degree
(n-1)+(n-3) = 2n-4,
and the r.h.s. is a polynomial of degree n-2. That gives an equation
2n-4 = n-2,
from which n=2 that contradicts to n being >=3.
Alec
solve for y'y'''=y'
solve for y'y'''=y'
Really?
If that's really what the equation is, it's trivial: the difference of the two sides of the equation factors as
y' (y''' - 1)
so there are two cases: y' = 0 and y''' = 1.
polysols
Maple has a command for this. It gives the single solution that Alec found above
d:=diff(y(x),x)*diff(y(x),x$3)-diff(y(x),x$2):
DEtools[polysols](d,y(x));
gives
[y(x) = _C3+_C4*x]