> Let B=[ [A I];
>         [0 I] ]
> 
> B^(k+1) = [ [A^k   I+A+...+A^k];
>             [0          I     ] ]

thanks