added LU Decomposition::inverse() method. Also changed solve() method to use regular arrays instead of vectors as the argument and return value, to reduce the number of vcector instances needed.