• API
• FAQ
• Tools
• Archive
SHARE
TWEET

# Untitled

a guest Jun 19th, 2019 61 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. NN = 200; R = 0.05;
2. xlist = Table[x, {x, -0.2 [Pi], 0.2 [Pi], 0.01}];
3. modl[n_] := 2*^-3 (Quotient[n, 2] - NN/2);
4. t1 = -1 + Cos[x] - I Sin[x] + I R; t1p = -1 + Cos[x] + I Sin[x] + I R;
5. t2a[n_] := -1 - modl[n]; t2b[n_] := -1 + modl[n];
6. mat[x_] =
7.   DiagonalMatrix[
8.     Table[If[EvenQ[n], t1, t2a[n]], {n, 0, 2 NN - 1 - 1}], 1] +
9.    DiagonalMatrix[
10.     Table[If[EvenQ[n], t1p, t2a[n]], {n, 0, 2 NN - 1 - 1}], -1] +
11.    DiagonalMatrix[
12.     Table[If[EvenQ[n], t2b[n], 0], {n, 0, 2 NN - 1 - 3}], 3] +
13.    DiagonalMatrix[
14.     Table[If[EvenQ[n], t2b[n], 0], {n, 0, 2 NN - 1 - 3}], -3];
15. list0 = Sort@Re@Eigenvalues[mat[xlist[[3]]]];
16. list0p = Table[list0[[i]] + list0[[2 NN - i + 1]], {i, NN}];
17. ListPlot[Tooltip@list0p, PlotRange -> All]
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.

Top