Mathematica 计算格林函数矩阵元：n+1 格点链的 G(n+1, 1)

以下是使用 Mathematica 编写的代码，计算 n+1 个格点的链的格林函数矩阵元 G(n+1, 1)：

ClearAll["Global`*"]

(*定义哈密顿量*)
H[n_] := SparseArray[{
   Band[{1, 2}, {n, n}] -> -1,
   Band[{2, 1}, {n, n}] -> -1,
   Band[{1, 1}] -> {0, 1}
  }]

(*定义格林函数*)
G[n_] := Inverse[{{ω + I ε, H[n]}, {-H[n], ω + I ε}}]

(*计算G(n+1, 1)*)
n = 10; (*取n=10为例*)
ω = 0; (*取ω=0*)
ε = 0.1; (*取ε=0.1*)
G[n + 1][[1, 1]]

执行以上代码，将得到 G(n+1, 1) 的值。可以根据需要修改 n、ω、ε 的值来计算不同情况下的格林函数矩阵元。