8.15-1(1)

fun(x, y) = sin(x) * sqrt(y) + x
q, e = integral2(fun, 0, 2, 0, 1)

8.15-1(2)

fun(z, y, x) = y * sin(x) + z * cos(x)
q, e = integral3(fun, 0, pi, 0, 1, -1, 1)

8.15-2

fun(y, x) = 1
ymin(x) = 0
ymax(x) = sqrt(4 - x^2)
q, e = integral2(fun, 0, 2, ymin, ymax)
P = q / 4