#!/usr/bin/python
import sys
import numpy as np
Nbins=int(sys.argv[2])
with open(sys.argv[1], 'r') as HistoFile:
    L = np.loadtxt(HistoFile,usecols=0,dtype=float)
    H=np.histogram(L,bins=Nbins); Bins=H[1][1] - H[1][0]
    B=np.delete(H[1], Nbins) + (Bins/2.)
    data=zip(B,H[0])
    np.savetxt(sys.argv[3], data)