У меня та же проблема, что и здесь https://comp.soft -sys.matlab.narkive.com / HV95Az6b / как вычислить нормальную кумулятивную вероятность с чрезвычайно высокой точностью
Я загружаю пакет 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 не работает, возможно, мне также следует отредактировать исходный кодэрф (или написать свой эрф)?
Заранее спасибо!