https://quadexcel.com/wp/useful-matrix-functions-in-r-programming-tutorial-13/