x = var('x')    # define a variable x
y = function('y', x)   # define y to be a function of that variable
DE = diff(y, x) + (3*y^2)/x
desolve(DE, [y,x],ics=(1,1/2))

Arnab Basu, Student, IIT Roorkee

Saurabh Agrawal, Student, IIT Roorkee

Rachit Jha, Student, IIT Roorkee

SageDays/Kreyszig-1.5-13 (last edited 2010-08-09 14:06:52 by saurabhagrawal)