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

k = var('k')

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

desolve(DE, [y,x])

SageDays/Kreyszig-1.3-6 (last edited 2010-08-11 11:49:15 by bhanukiran)