return c*pow(x, c-1)*np.exp(-pow(x, c))