For sure. I just wanted to point out the approach that arguably required "the least math knowledge".

Also, OP was solving a degree four problem and wanted a closed form.

For better stability just use `x = np.linalg.solve(M, y)`, or whichever linear solver you expect to be fast.