In this paper, we propose a numerical algorithm to solve a class of dynamic time scale equation which is called the q-difference equation. First, we apply the method for solving initial value problems (IVPs) which contain the first and second order delta derivatives. Illustrative examples show the usefulness of the method. Then we present applications of the method for solving the strongly non-linear damped q-difference equation. The results show that our method is more accurate than the other existing method.