# solution 3 : evaluation in double integral var('u') var('v') uu = integral(u*u+v*v, u, 0, 2) vv = integral(uu, v, 0, 2) result = vv/4 print 'result of double integral : ', result