Added a struct that holds the LU decomposition and has a 'solve' member function. Fixed some typos and added some clarifications to the README file.