Numerical solution of a human circulation model