Update the derivation of Crank-Nicolson model.

Includes changes to indexing notation to avoid
conflict with "i" as used in stability analysis
for complex numbers. The non-linear source term
is handled better.