clear;

a = 1;
b = 3;

// *** 関数の定義 ***
// 関数f(x)
deff('y = f(x)','y = exp(- x .^ 2)');
// 原始関数F(x)
deff('y = F(x)','y = sqrt(%pi) .* erf(x) ./ 2');

X = linspace(0,5);

// *** グラフの描画 ***
plot(X,f(X));
xlabel("x");
ylabel("f(x)");

// *** 積分の計算 ***
// 解析的な積分
F(b) - F(a)
// 数値積分
integrate('f(x)','x',a,b)
// モンテカルロ積分
n = 1000;
mean((b - a) * f( a + (b - a) * rand(1,n)))