McCullagh_compute_condition function

Compute the linear constraint on psi elements for identifiablity.