Numerical solution of parameter-dependent linear systems