def solve_least_squares(self, rhs, method='CH'):
        if method == 'CH':
            return self.cholesky_solve(rhs)
        t = self.T
        return (t*self).inv(method=method)*t*rhs
Return the least-square fit to the data.
dataset/ETHPy150Open sympy/sympy/sympy/matrices/matrices.py/MatrixBase.solve_least_squares