private static double computeJSD(double[] px,double[] py){
  double[] m=computeMean(px,py);
  return 0.5 * (computeKLD(px,m) + computeKLD(py,m));
}

private static double computeJSD ( double [ ] px , double [ ] py ) { double [ ] m = computeMean ( px , py ) ; return 0.5 * ( computeKLD ( px , m ) + computeKLD ( py , m ) ) ; }
compute the jensen - shannon divergence
