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

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

desolve(DE, [y,x])

SageDays/Kreyszig-1.3-9 (last edited 2010-08-11 11:54:17 by bhanukiran)