gamma function for N>100
9 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Ladan Avazpour
am 8 Okt. 2018
Kommentiert: Ladan Avazpour
am 9 Okt. 2018
Hi, everyone.
I have written a simple code by using a gamma function. when I choose N>98, the answer is Nan. What should I do to fix this problem?
function sii = diagonals(n, al, alp, bet)
ssi = zeros(1,n+1);
for i = 1:n+1
sii(i) = (2^(al+1)*(gamma(i-1+alp+1))*(gamma(i-1+bet+1)))/((2*(i-1)+al+1)*(factorial(i-1))^2);
end
end
2 Kommentare
John D'Errico
am 8 Okt. 2018
Bearbeitet: John D'Errico
am 8 Okt. 2018
Is there a good reason why you have not just used gamma itself? Why would you write your own poor code to compute something that is computed well already for you? I'm not going to try to decrypt what you are doing based on undocumented, poorly written code. Perhaps you are trying to write your own incomplete gamma. Oops, also exists already, part of MATLAB.
Akzeptierte Antwort
Walter Roberson
am 9 Okt. 2018
If gammaln cannot be made to suit your purpose then you should use the symbolic toolbox or equivalent
Weitere Antworten (1)
John D'Errico
am 8 Okt. 2018
No. You will not make it work for large values, when they exceed the limits of gamma itself. Your code is not better than what you will find inside gamma itself. Sorry, but NOT.
Instead, consider using the function gammaln, the natural log of the gamma function. In essentially all cases when gamma produces an overflow, gammaln will be the solution, by re-writing what you are doing to accommodate logs. Only consider exponentiating at the end.
gammaln(10000)
ans =
82099.7174964424
Those big numbers are meaningless in a double precision context, because they would vastly overflow the dynamic range of a double precision number.
Siehe auch
Kategorien
Mehr zu Gamma Functions finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!