I=eye(2,2) N=I-A [C,kerN]=linsolve(N,-B)