To solve this with Python, we can use the SciPy Scientific Python package to compute the Moore-Penrose pseudo-inverse \(\text{A}^+\)of the matrix \(\text{A}\), where \(\text{A}^+\) is defined the matrix that 'solves' (by least squares minimization) the matrix equation \(\text{Ax=y}\), with solution  \(\overline{x}=\text{A}^+y\).  Sample code is provided below: