change double density_negbinomial(double x,double mu,double tau) to double negbinomial_density(double x,double mu,double tau) in cnegbin.cpp