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: