coeff_mat = scipy.sparse.coo_matrix((data, (row, col))).tocsc()
return scipy.sparse.linalg.spsolve(coeff_mat, np.ones(2 * (n - 1)) * n)
Использование scipy разреженных матриц для решения системы уравнений