x = var('x') 
y = function('y', x)
deq1 = diff(y , x) + 2*x*y == 0
r=  desolve(deq1, [y,x], ics = (0,1)) 
show(r)

Kreyszig-1.3-5 (last edited 2010-12-18 06:11:35 by sai kumar)