1/(sqrt(pi*2* [1.00] ^2)) * integrate( e^(-(x- [0.00] )^2/(2* [1.00] ^2)) * 0.5*(1+erf(x/sqrt(2))), {x,-inf,inf})