Trench Inverse Update - testing

Initialization

<< FGN.m

Testing

rA = Table[FGNAcf[k, 0.7], {k, 0, 6}];
r = Take[rA, 5];
GI = Inverse[G = toe[r]];
rnew = 
        rA[[6]];
MatrixForm[Inverse[toe[Take[rA, 6]]]]

TrenchInverseUpdate[GI, r, rnew] // MatrixForm

Flatten[
        Inverse[toe[Take[rA, 6]]] - TrenchInverseUpdate[GI, r, rnew] // Abs] // 
      Total