fullrfk - full rank factorization of A^k
This function computes the full rank factorization of A^k i.e. Bk*Ck=A^k where Bk is full column rank and Ck full row rank. One has range( Bk )=range( A^k ) and ker( Ck )=ker( A^k ).
For k=1 , fullrfk is equivalent to fullrf .
A=rand(5,2)*rand(2,5);[Bk,Ck]=fullrfk(A,3); norm(Bk*Ck-A^3,1)
F.D (1990);