Reproduce two of the curves from fig 3b
Use np.power