I think that the docstring from the abstract class should be inherit to the child classes

needs further testing

numpy transpose is not working on 1D vectors [it returns the same vector]

We should save the inverse of the approximated hessian right ?

since it is dynamically allocated len(points) = steps or maximum_iterations

If we pass them error rises when we choose other line search i.e. inexact