行列 A のルート A^{1/2} の計算について
A^{1/2} = {I + (A-I)}^{1/2}
= I + (1/2)*(A-I) - (1/8)*(A-I)^2 + (1/16)*(A-I)^3 - ・・・
自分はこれで収束する場合に近似(数値)計算するくらいしか考えつかないのですが、
Mathematicaだと MatrixPower[A,1/2] で、なにやら厳密解らしき 値がでてきます。
一体どういう計算を行っているのでしょうか？ Mathematicaで採用されている方法でなくてもいいので教えてください。