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

Kreyszig-9.6-6 (last edited 2010-12-15 12:58:23 by 172)