powered by
Computes phi based on gamma
McCullagh_compute_phi(gamma, j)
phi[j]
vector of gamma parameters
index of phi to compute