swh:1:snp:d91631a580497fff421f83ef3e706159a0d67602
Tip revision: ca0cbd288f03c7a29050b3d6698a96b45ccfa4b2 authored by Marko Laine on 15 August 2021, 17:59:37 UTC
cleanup
cleanup
Tip revision: ca0cbd2
betapf.m
function y=betapf(x,a,b)
%BETAPF Beta probability density function
% y = betapf(x,a,b) returns G(a+b)/G(a)/G(b)*x^(a-1)*(1-x)^(b-1)
% where G is the gamma function
% Marko Laine <marko.laine@fmi.fi>
% $Revision: 1.5 $ $Date: 2012/09/27 11:47:33 $
%y = gamma(a+b)/gamma(a)/gamma(b)*x.^(a-1).*(1-x).^(b-1);
y = x.^(a-1).*(1-x).^(b-1)./beta(a,b);