#Solution of Problem Set 2.6 Problem no 15

reset()
#resets the initial values
x = var('x')
y = function ('y',x)
ics=[1,2,-6]
eq=4*x^2*(diff(y,x,2 ))+24*x*diff(y,x)+25*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 "is"
print
p=desolve(eq,y,ics)
print "y =" ,
show(p)
print
plot(p,(x,0.2,1.5))

SageDays/Kreyszig-2.6-15 (last edited 2010-08-12 09:46:45 by RahulPalkar)