x = var('x')   
y = function('y',x)   
sol=desolve (x*diff(y,x) + 4*y - 8*x^4, [y,x], ics=(1,2))
sol

SageDays/Kreyszig-1.6-22 (last edited 2010-08-11 10:52:20 by SWATI AHUJA)