numerical solution of linear systems