sympy symbolic matrix racine carrée

Question

Est-il possible d'utiliser la fonction sympy sqrtm pour calculer la racine carrée d'une matrice symbolique? Par exemple, si j'ai la matrice [[1, a], [a, 1]] avec un> 0, puis-je utiliser sympy pour calculer la racine carrée?

Answer 1

from sympy import * 
a=symbols('a') 
m=Matrix([[1, a], [a, 1]]) 
m**(1/2) 

donne

Matrix([ 
[ (-a + 1)**0.5/2 + (a + 1)**0.5/2, -(-a + 1)**0.5/2 + (a + 1)**0.5/2], 
[-(-a + 1)**0.5/2 + (a + 1)**0.5/2, (-a + 1)**0.5/2 + (a + 1)**0.5/2]]) 

Pour comparer les résultats, dans Mathematica