Added the dstrauss function for approximate density.