这段代码最后有点问题，输出不对

非常抱歉，上一段代码确实有问题。以下是修改后的代码：

import numpy as np
import sympy

# 定义下三角矩阵L和列矩阵b
L = np.array([[1, 0, 0], [2, 3, 0], [4, 5, 6]])
a1, a2, a3 = sympy.symbols('a1 a2 a3')
b = sympy.Matrix([[a1], [a2], [a3]])

# 求解Ly=b中的y
n = L.shape[0]  # 矩阵L的行数
y = sympy.zeros(n, 1)  # 初始化y
for i in range(n):
    y[i] = b[i]
    for j in range(i):
        y[i] -= L[i, j] * y[j]
    y[i] /= L[i, i]

# 定义上三角矩阵U
U = np.triu(L)

# 求解Ux=y中的x
n = U.shape[0]  # 矩阵U的行数
x = sympy.zeros(n, 1)  # 初始化x
for i in range(n-1, -1, -1):
    x[i] = y[i]
    for j in range(i+1, n):
        x[i] -= U[i, j] * x[j]
    x[i] /= U[i, i]

# 将y和x中的数值替换为符号变量
y = y.subs({a1: 1, a2: 2, a3: 3})
x = x.subs({a1: 1, a2: 2, a3: 3})

# 输出结果
print("y = \n", y)
print("x = \n", x)

在修改后的代码中，我们将y和x中的数值直接替换为了1、2和3，这是因为这个例子中的符号变量a1、a2和a3对应的实际值就是1、2和3，这样就可以验证我们的代码是否正确了。在输出结果时，我们直接输出了y和x的值，而不是使用sympy库中的evalf()函数进行数值运算，这样就可以得到带有符号的解。

如果需要对其他的符号变量进行求解，只需要将a1、a2和a3替换成对应的符号变量即可。