clear;

// *** 関数の定義 ***
// 関数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));

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