Create functions to solve and cache matrix inversions
Add documentation