Matlabを用いて正規分布で表示された確率密度関数を積分する

P = {\displaystyle \int_{-\infty}^\infty } \dfrac{1}{\sqrt{2\pi}\sigma} \exp \left\lbrack -\dfrac{1}{2}\left( \dfrac{x-\mu}{\sigma} \right)^2 \right\rbrack 

P = @(mu,sigma) (integral(@(x) (1/(sqrt(2*pi)*sigma)*exp(-1/2*((x-mu)/sigma).^2)),-Inf,Inf))

実行結果

>> P(0,1)
ans =
    1.0000