function phiIcross = cros(phiI)
phiIcross = inv(phiI'*phiI)*phiI';