Application of sparse matrix solvers and Newton's method to fluid flow problems