build_phi_matrix function

Build intermediate C(phi,n) matrix in inversion of Exchangeable variance matrix

build_phi_matrix(n, phi, directed = T, sta_flag = F)