#solve and permutation matrix from #factorize_lu

#factorize_lu would return a permutation array in a very strange manner, which left me confused until https://github.com/andrewcsmith clarified it. So I've put a provision to return the permutation matrix along with the result.

Linear algebra equations can now be solved LAPACK free with the new #solve method.