M = var('M')
a = var('a')
u = var('u')
v = var('v')
0<=u<=2*pi
-pi/2<=v<=pi/2
temp = M/(4*pi*(a**2))
show(expand(temp))
in1 = integral((temp*(a**4)*(cos(v)**3)),u,0,2*pi)
show(expand(in1))
ans = integral(in1,v,-pi/2,pi/2)
show(expand(ans))

Kreyszig-9.6-7 (last edited 2010-12-15 13:00:04 by 172)