A numerical model based on auxiliary difference equation (ADE) is derived to simulate the plasma in the time domain and the calculating results are compared with that using Z transform.