x = var('x')    
y = function('y',x)
DE= 4*(x^2)*diff(y,x,2)-(3*y)
desolve(DE,[y,x],ics=[1,3,2.5])

SageDays/Kreyszig-2.1-19 (last edited 2010-08-10 13:20:14 by manojmanu)