在Octave中实现错误功能

Question

我是编程的完全初学者，我在执行error function作业时遇到了麻烦。我的erf(x)适用于x的值，从0到2.0，但不适用于更高的值。适当值的表格为here。

function y = erf (z) 
    konst = 2/(pi^0.5);
    vsota=0;

    n=0;
    while n <= 500;
        vsota = vsota + ((-1)^n)*(z^(2*n+1))/(factorial(n)*(2*n+1));
        n=n+1;
    end
    y=konst*vsota;
end

请帮我stackoverflow，你是我唯一的希望！

Answer 1

也许octave在评估大数字factorial()时遇到问题。我在Python中使用n <= 170得到了很好的结果。