Numerical Stability of Finite Difference Algorithms for Electrochemical Kinetic Simulations. Matrix Stability Analysis of the Classic Explicit, Fully Implicit and Crank-Nicolson Methods, Extended to the 3- and 4-point Gradient Approximation at the Electrodes.