u = var('u')
v = var('v')
0<=u<=2*pi
-pi/2<=v<=pi/2
a = var('a')
in1 = integral(a**2*cos(v), u, 0, 2*pi)
print "First half",in1
fp1 = integral(in1,v,-pi/2,pi/2)
print "Second half",fp1
print "Final Result: ",show(expand(fp1))

Kreyszig-9.6-5 (last edited 2010-12-15 12:54:44 by Enhance_vasu.naresh)