Advertisement
Guest User

Untitled

a guest
Jul 25th, 2017
100
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.30 KB | None | 0 0
  1. n = 10000;
  2. M = N[SparseArray[
  3. Join[Table[{i, i} -> 1.0, {i, 1, n}],
  4. Table[{i, i - 1} -> 2.0, {i, 2, n}],
  5. Table[{i, i + 1} -> 2.0, {i, 1, n - 1}]]]];
  6. val = Eigenvalues[M, -1][[1]]
  7.  
  8. A=coo_matrix((data, (row, col)), shape=(nn, nn)).toarray()
  9. eigs(A, 1, sigma=0)[0][0]
  10.  
  11. (-0.00013788630102868301+0j)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement