InvP:

--------------------------------------------------------------------------
   The product form of the inverse:

               [ 1     v1]-1
   inv(Bnew) = [    1  v2]  inv(Bold)
               [       v3] 

   v = inv(Bold)y where y is the new column k of B

   Since version 1.
--------------------------------------------------------------------------
   Form:
   bInv = InvP( bInvOld, v, k )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   bInvOld       (n,n)    Old inverse
   v             (n,1)    Vector = inv(Bold)*y where y is a new column of B
   k                      Column index of v

   -------
   Outputs
   -------
   bInv          (n,n)    New inverse

--------------------------------------------------------------------------