def JCB(A,b,N=25,x=None):
if x is None:
x = zeros(len(A[0]))
D = diag(A)
R = A - diagflat(D)
for i in range(N):
x = (b - dot(R,x))/D
pprint(x)
return x
微信
公众号
扫码关注公众号