code: private double anomalyScore(double x, double mean, double sd) {
          x = (x - mean) / sd;
          x = normPDF(x);
          x = (max_norm_prob - x)/max_norm_prob;
          x *= 100;
          return x;
      }
docs: Output: anomaly score between 0 and 100