4x1+x2=1 x1+4x2+x3=2 x2+4x3+x4=3 ............ x998+4x999+x1000=999 x999+4x1000=1000求解该方程组的代码

要求解该线性方程组，可以使用NumPy库中的线性方程组求解函数`numpy.linalg.solve()`。下面是使用NumPy的示例代码：

import numpy as np

# 构造系数矩阵 A
A = np.zeros((1000, 1000))
for i in range(1000):
    A[i, i] = 4
    if i < 999:
        A[i, i+1] = 1

# 构造常数向量 B
B = np.zeros(1000)
B[0] = 1
B[1] = 2
B[2:-2] = np.arange(3, 999)
B[-2] = 999
B[-1] = 1000

# 求解线性方程组
X = np.linalg.solve(A, B)

print(X)

运行上述代码，将得到线性方程组的解向量 X。

希望这可以帮助到您！如果还有其他问题，请随时提问。