3 レポート課題13

締め切りは 7月20日(金) 18:00。Oh-o! Meiji レポートシステムを使って、 kadai13.pdf を提出する。

neumann.m を参考に、 行列の指数関数

$$\exp A=\sum_{n=0}^\infty \frac{1}{n!}A^n$$

を計算する関数 exponential() を作って下さい。 例えば

$$A=t J,\quad J:=\begin{pmatrix} 0 & 1\\ -1 & 0 \end{pmatrix}$$

のとき、 $\exp A$ が何になるか調べて (対馬先生のテキストの付録 http://nalab.mind.meiji.ac.jp/~tsushima/senkei.html の第9章の章末問 題2参照)、それを MATLAB の計算結果と比較しなさい ($t$ は適当に 1 つ 2 つ選べば良い。$t=1$ とか $t=\pi/4$ とか。)。