x = var('x')
y = function('y',x)

DE  = diff(y,x) + 3*(x^2)*(y^2)

desolve(DE, [y,x])

SageDays/Kreyszig-1.3-4 (last edited 2010-08-11 11:44:58 by bhanukiran)