How can i solve this equation to n ?
조회 수: 1 (최근 30일)
이전 댓글 표시
I try to solve this equation 1-(γ(a,x)/(n-1)!)=R to n but i can't find how, where R>=0 and γ(a,x) is the lower incomplete gamma function as i found. Does anyone have any idea?
댓글 수: 0
답변 (4개)
Walter Roberson
2012년 2월 3일
If 1-(γ(a,x)/(n-1)!)=R then 1 + R = y(a,x) / (n-1)! and so
(n-1)! = y(a,x) / (1 + R);
(n-1)! is gamma(n), so you want to solve gamma(n) = y(a,x) / (1+R)
c = gammainc(a,x) / (1+R);
fzero( @(n) gamma(n) - c, 5 )
댓글 수: 1
Sean de Wolski
2012년 2월 3일
Definitely more optimized than mine. One comment: gammainc() expects (x,a) not(a,x).
Sean de Wolski
2012년 2월 3일
fzero(@(n)1-(gammainc(x,a)/gamma(n))-R,5)
Maybe? You'll hit overflow if you start with a big n.
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Gamma Functions에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!