  The iterated Crank-Nicolson is a predictor-corrector algorithm commonly used
in numerical relativity for the solution of both hyperbolic and parabolic
partial differential equations. We here extend the recent work on the stability
of this scheme for hyperbolic equations by investigating the properties when
the average between the predicted and corrected values is made with unequal
weights and when the scheme is applied to a parabolic equation. We also propose
a variant of the scheme in which the coefficients in the averages are swapped
between two corrections leading to systematically larger amplification factors
and to a smaller numerical dispersion.
