Remark 1.

To run implicit time stepping the value of $\epsilon$ is set to about the size of the time step, $dt$ divided by ten (or smaller).

To run explicit time stepping the value of $\epsilon$ is set to 1.