Mathematica 计算 3 格点链格林函数矩阵元教程

以下是使用 Mathematica 计算 3 格点链的格林函数矩阵元的示例代码：

首先，我们定义 3 格点链的哈密顿量：

h = {{0, 1, 0}, {1, 0, 1}, {0, 1, 0}};

接下来，我们定义能量：

e = 1;

然后，我们定义单位矩阵：

id = IdentityMatrix[3];

接下来，我们定义格林函数矩阵元的表达式：

g[i_, j_] := Inverse[e id - h][[i, j]];

最后，我们计算格林函数矩阵元：

gmatrix = Table[g[i, j], {i, 1, 3}, {j, 1, 3}];

运行以上代码，就可以得到 3 格点链的格林函数矩阵元。