A numerical solution of axisymmetric cavity flows