Computer solution of large sparse positive definite systems