var('x')
y=function('y',x)
print desolve(diff(y,x) + y/3 - 1/3*(1-2*x)*y^4,[y,x])