Added polynomial solvers from 1 to 4 degrees and fixed
matrix_solve.