ddir function

Dirichlet density