Numerical solution of unsteady incompressible Navier-Stokes equations using high-order method of lines