clear;

r = 1;

// *** 積分する関数の定義 ***
function z = func(x,y)
  z = (x .^ 4) * sin(y)
endfunction
// 積分範囲
x0 = 0; x1 = r;
y0 = 0; y1 = %pi;

// *** 二重積分 ***
X = [x0, x0;
     x1, x1;
     x1, x0];
Y = [y0, y0;
     y0, y1;
     y1, y1];
// 数値解
I = int2d(X, Y, func)
// 解析解
A = 2*(r^5)/5