Crank-Nicolson applied to a plug profile

Crank-Nicolson never blows up, so any $F$ can be used (modulo loss of accuracy).

$N_x=50$. $F=5$ gives instabilities.

$N_x=50$. $F=0.5$ gives a smooth solution.