x=var('x')
y=function('y',x)
deq1=diff(y,x)+(y/3)-((1/3)*(1-2*x)*y^4)==0
r=desolve(deq1,[y,x],y)
t=r.simplify_full()
t.simplify_log()