You can do this in Python with 1 function call:
https://docs.scipy.org/doc/scipy/reference/generated/scipy.n...
https://docs.scipy.org/doc/scipy/reference/generated/scipy.n...
If you don't have that knowledge, I think you can get quite far by just starting with a unit matrix and changing the coefficients iteratively until you get something that looks right.