The Kelvin-Helmholtz mode instability driven by an equilibrium radial electric field in a cold plasma is studied numerically in cylindrical geometry. The numerical analysis is applicable to both toroidal and mirror plasmas. The growth rate and a threshold in the radial gradient scale length and also the magnitude of the E × B drift frequency for instability are obtained.