def phi1(x, alpha):
    return ((np.log((tau(x, alpha) + 1)) + (alpha * np.log(((tau(x, alpha) / alpha) + 1)))) - tau(x, alpha))
def phi1(x, alpha):
    return ((np.log((tau(x, alpha) + 1)) + (alpha * np.log(((tau(x, alpha) / alpha) + 1)))) - tau(x, alpha))