/*fabden called by dlaw command, returns density for any density law */
double fabden(double radius, 
  double depth);