In [1]:

    

import numpy as np
import matplotlib.pyplot as plt
import scipy.stats


    

In [142]:

    

a = 3
x = np.linspace(0, 5, 1e7)
y = (a ** 3/2) * x ** 2 * np.exp(-2 * a * x)


    

In [143]:

    

##plt.plot(x, y, '.')
plt.show()


    

In [144]:

    

ymedian = np.mean(y)
idx = np.nanargmin(np.abs(y-ymedian))


    

In [145]:

    

5 * idx / 1e7


    

    
Out[145]:





1.0716665000000001

In [146]:

    

ymedian = np.max(y)
idx = np.nanargmin(np.abs(y-ymedian))


    

In [147]:

    

idx / 1e7 * 5


    

    
Out[147]:





0.33333349999999995

In [ ]: