数値計算ライブラリのnumpyはほんとに便利。 import numpy as np def jacobi(A,N,check): B = np.fabs(A - np.diag(list(np.diag(A)))) nondiagmax = np.max(B) while nondiagmax > check: k = int(np.argmax(B) / 3) m = np.argmax(B) % 3 cos2phi = np.fab…