In this study a robust iterative technique is developed for solving the one-dimensional human arterial blood flow problem by adopting Locally Conservative Galerkin method (LCG). Using Newton method with two different linear solvers (i.e. Gauss elimination and Jacobi methodologies), the non-linear governing equations are solved. Such strategies result in rapid convergence and fast solution without excessive memory cost for semi and full implicit LCG discretizations. In the proposed methods, the numerical strategies require computing a 4 x 4 matrix per element to determine blood flow characteristics. The novel methods developed are employed to study the blood flow through the major vessels of a complex human systemic circulation network. The validity and stability of the present methods are investigated by comparing the results against the available data in the literature.