オクターブでのエラー機能の実装

Question

私はプログラミングの初心者です。私は宿題のためにerror functionを実装することに問題があります。私のerf(x)は、0から2.0までのxの値に対して正常に動作しますが、それ以上の値では機能しません。適切な値の表は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 

私は唯一の希望です！

Answer 1

多分、オクターブは大きい数字のfactorial()を評価するのに問題があります。私はn <= 170を使ってPythonで良い結果を得ました。