Implementation of Quasi Newton methods