x = var('x')   
y = function('y',x)   
sol=desolve (diff(y,x) - 1- y^2, [y,x], ics=(0,0))
sol

SageDays/Kreyszig-1.6-21 (last edited 2010-08-11 10:53:54 by SWATI AHUJA)