（2）使用克莱姆法则求解二元一次方程组。 [25.6x+3y=28 [8x+13.4y=65给出Python代码

克莱姆法则（Cramer's Rule）用于求解线性方程组时，当系数行列式不为零时，可以直接计算出每个变量的值。对于给定的二元一次方程组：

\[ \begin{align*}
a_1 x + b_1 y &= c_1 \\
a_2 x + b_2 y &= c_2
\end{align*} \]

其中 \( a_1, b_1, a_2, b_2 \) 是方程中的系数，\( c_1, c_2 \) 是常数项，克莱姆法则表示解为：

\[ x = \frac{\text{Determinant of the coefficient matrix with column 1 replaced by } [c_1]_{(1)} }{\text{Determinant of the coefficient matrix}} \]
\[ y = \frac{\text{Determinant of the coefficient matrix with column 2 replaced by } [c_2]_{(2)}}{\text{Determinant of the coefficient matrix}} \]

这里，`[c_1]_{(1)}` 和 `[c_2]_{(2)}` 表示将原方程组的第一列和第二列分别替换为常数列。

下面是用 Python 实现的代码片段：

import numpy as np

# 定义方程组的系数矩阵和常数项
coeff_matrix = np.array([[25.6, 3], [8, 13.4]])
constant_terms = np.array([28, 65])

# 计算系数矩阵的行列式
determinant_coeff = np.linalg.det(coeff_matrix)

# 计算替换后的行列式，用于计算 x 和 y
det_column_1_replaced = np.linalg.det(np.column_stack((coeff_matrix[:, 0], constant_terms)))
det_column_2_replaced = np.linalg.det(np.column_stack((coeff_matrix[:, 1], constant_terms)))

# 检查是否能应用克莱姆法则，即行列式是否非零
if determinant_coeff != 0:
    x = det_column_1_replaced / determinant_coeff
    y = det_column_2_replaced / determinant_coeff
else:
    print("方程组无唯一解或无穷多个解，因为系数行列式的值为零。")

# 输出结果
print(f"x = {x}, y = {y}")

执行这段代码后，如果系数行列式不为零，你会得到 \( x \) 和 \( y \) 的值。如果行列式为零，则说明方程组有无穷多解或者没有唯一解。