from sympy import * import sympy var('ax ay bx by') a = sympy.Matrix([ax, ay]) b = sympy.Matrix([bx, by]) S = sympy.sqrt(a.norm()**2 * b.norm()**2 - a.dot(b)**2) S = sympy.simplify(S) myValue={ax:1, ay:1,bx:1-sympy.sqrt(3), by:1+sympy.sqrt(3)} print(f'# 問題2の解答 面積は {S}') print(f'# 問題2の解答 面積は {S.subs(myValue).simplify()}') # 問題2の解答 面積は sqrt(-(ax*bx + ay*by)**2 + (Abs(ax)**2 + Abs(ay)**2)*(Abs(bx)**2