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. OK, I Understand
 
Top