new_contents: import hpat
  from numba import prange
  import numpy as np
  import h5py
  import argparse
  import time
  
  @hpat.jit(locals={'X':hpat.float64[:]})
  def kde():
      f = h5py.File("lr.hdf5", "r")
      X = f['responses'][:]
      b = 0.5
      points = np.array([-1.0, 2.0, 5.0])
      N = points.shape[0]
      n = X.shape[0]
      exps = 0
      t1 = time.time()
      for i in prange(n):
          p = X[i]
          d = (-(p-points)**2)/(2*b**2)
          m = np.min(d)
          exps += m-np.log(b*N)+np.log(np.sum(np.exp(d-m)))
      t = time.time()-t1
      print("Execution time:", t,"\nresult:", exps)
      return exps
  
  def main():
      parser = argparse.ArgumentParser(description='Kernel-Density')
      parser.add_argument('--size', dest='size', type=int, default=1000000)
      args = parser.parse_args()
      size = args.size
  
      res = kde()
  
  if __name__ == '__main__':
      main()
  
subject: add kde example