A new algorithm, which is exceptionally fast for certain choices of numerical parameters, is described for the study of nonlinear, incompressible flow between two rotating disks. Typical examples for Reynolds number R in the range 10 ≤ R ≤ 2000 are described and discussed. Comparisons are made with the limited available results generated by other methods.