Résumé en anglais

The method of simultaneous inverse iterations is unstable when it is used to compute the invariant subspace of a defective eigenvalue. In this paper we propose a stable method, of equal complexity, which is based on a modification of Newton 's method.