Computing the maximum robust controlled invariant subspace