#Solution of Problem Set 2.6 Problem no 19

reset()
r = var('r')
bc=[5,30,10,300]
print "Let 'r' denote the radius of the spheres"
v = function ('v',r) # 
print
print"Let 'v' be the electrostatic potential between the 2 concentric spheres"
print
print"""The physical information is v(r) is a solution of rv"+2v'=0"""
eq=r*diff(v,r,2)+ 2*diff(v,r)==0
show(eq)
print " The General Solution of the equation written above is"
desolve(eq,v)
print
print "The boundary condition of the distance and electrostatic voltages is"
print "at r =", bc[0],"v =",bc[1]
print "at r =", bc[2],"v =", bc[3]

p=desolve(eq,v,bc)
show(p)
plot(p,(r,5,10))

SageDays/Kreyszig-2.6-19 (last edited 2010-08-12 10:12:36 by RahulPalkar)