By iteration, we further get that u t ≤ u 0 · (E + IP) t , where E is a unit matrix.