Added a function to construct Lagrange interpolation functions
Fixed the return of poly, which returned arrays in reversed order
Added documentation for the lagrange function