/* 

MLE for dirichlet distribution.


to compile : 


mex -g -output dirichlet_mle.dll dirichlet_mle.c


mex -f mexopts_intel10amd.bat -output dirichlet_mle.dll dirichlet_mle.c



Example 


d         = 4;
N         = 1000;
delta     = ceil(10*rand(d , 1));
p         = dirirnd(delta(: , ones(1 , N)));
delta_est = dirichlet_mle(p);

[delta , delta_est]



  Author : S閎astien PARIS