Numerical solution for unsteady Couette flow of viscoplastic fluids