# 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

Kreyszig-9.3-2 (last edited 2010-12-15 13:01:02 by 172)