# Solution 2 : page 483 # double integral in polar coordintes var('r') sin2 = sin(x)*sin(x) rr = integral(r**3, r, 0, 1) ss = integral(sin2, x, 0, pi/2) Ix = ss*rr print 'Ix : ', Ix