clear;

// *** 積分する関数の定義 ***
function z = func(x,y)
  region = y >= x .^ 2
  z = - 1 ./ ((2 * x + y + 1) .^ 2) .* region
endfunction

// *** 二重積分 ***
X = [0;
     1;
     1];
Y = [0;
     0;
     1];
// 数値解
I = int2d(X, Y, func)
// 解析解
A = (1/3) * log(4) - 1/2