clear;

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

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