len_mono = 0
Xfit_mono = zeros(Xfit.size)
Xfit_mono_ind = zeros(Xfit.size)
Xfit_mono[len_mono] = Xfit[0]
for iX, x in enumerate(Xfit):
    if iX > 0:
        if x > Xfit_mono[len_mono]:
            len_mono = len_mono + 1
Xfit_mono[len_mono] = x
Xfit_mono_ind[len_mono] = iX
Xfit_mono_ind = Xfit_mono_ind[:len_mono + 1]
Xfit_mono = Xfit_mono[:len_mono + 1]
print(Xfit_mono_ind)
print(Xfit_mono)
Fit a monotonic curve (preferably in python)