Numerical approach to the Stokes problem with high contrasts in viscosity

An algorithm based on Sherman–Morrison–Woodbury formula to solve numerically the Stokes problem is described. The algorithm allows to obtain the solution for viscosity contrasts up to ten orders of magnitude, moreover solution speed does not depend on viscosity contrast. Tests of accuracy of the algorithm are provided.