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()

Kreyszig-1.6-33-U (last edited 2010-12-18 06:16:43 by 10)