Quasi-Newton method