2.7（2）

• 代码

点击查看代码

import numpy as np

n=1000
A=np.zeros((n,n))
b=np.arange(1,n+1)

np.fill_diagonal(A[1:,:-1],1)
np.fill_diagonal(A[:-1,1:],1)
np.fill_diagonal(A[:-1,1:],1)

x=np.linalg.solve(A,b)
print("解为：")
print(x[:5])
print(x[-5:])
print("学号：3015")