A numerical technique has been developed to solve the time-dependent equation of continuity of plasma for the ionosphere. The solution of the equation gives the behaviour of electron density at any time of the day of the year. A computer program for this has been developed which has been used to simulate different physical situations of the ionosphere.