Блог

У меня та же проблема, что и здесь https://comp.soft-sys.matlab.narkive.com/HV95Az6b/how-to-compute-the-normal-cumulative-probability-at-extremely-precision

Я загружаю пакет https://www.advanpix.com/documentation/users-manual / и знать, как использовать его в простом выражении

 mp('pi/4')

ans = 

    0.78539816339744830961566.....
  

Но я не знаю, как использовать этот mp. в функции normcdf я пытаюсь отредактировать выражение последней строки в следующем коде, но пока не получается.

 function uu = normcdf(xx, mu, sigma)
%NORMCDF Gaussian CDF without using the stats toolbox
%     uu = normcdf(xx, mu, sigma)
% Iain Murray, July 2012
if ~exist('mu', 'var')
    mu = 0;
end
if ~exist('sigma', 'var')
    sigma = 1;
end
uu = 0.5*erf((xx-mu)./(sigma*sqrt(2)))   0.5;
  

Добавить mp() после 0.5 не работает, может быть, мне также следует отредактировать исходный код erf (или написать свой собственный erf)?

Заранее спасибо!