xc1 = np.power(x, -c)