Mathematica 9 includes algorithms for efficiently evaluating matrix functions. For a function with a power series for , the matrix function for an × matrix is an × matrix effectively based on the power series approximation using matrix powers of , but can be computed more accurately and efficiently using advanced algorithms.
In this example, the fundamental trig identities and will be tested for random matrices .
Compute the matrix sine and cosine.
Show the magnitudes of and .
Theoretically, the plot on the left should just have the values along the diagonal. The noise that appears comes from numerical error and is actually quite small, but shows up because the logarithmic default color function used by MatrixPlot shows values with vastly different orders of magnitude.
Testing the error in the matrix norm shows the numerical error is quite small.