#Solution of Problem Set 2.6 Problem no 16

reset()
#resets the initial values
x = var('x')
y = function ('y',x)
ics=[1,2,0]
eq=x^2*(diff(y,x,2 ))+x*diff(y,x)+9*y==0
show(eq)

print " The General solution of the above equation is "
p=desolve(eq,y,show_method=True)
show(p[0])
print " The method of solution is",p[1]
print
print "The Particular solution for the above equation when"
print
print "y(",ics[0],") =",ics[1],"and Dy(",ics[0],") =", ics[2] ,
print
p=desolve(eq,y,ics)
print "y =" ,
show(p)
print
plot(p,(x,0.01,1))

SageDays/Kreyszig-2.6-16 (last edited 2010-08-12 09:50:12 by RahulPalkar)