We consider the differential-algebraic system for the blood flow and pressure in the systemic arteries. By the operator splitting method, we transform the system into the hyperbolic one, introduce the bicharacteristics, and perform the time-space nonuniform discretization, obtaining the innovative difference scheme. Our results are illustrated with numerical experiments.