Advertisement
Guest User

Untitled

a guest
Feb 22nd, 2017
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.38 KB | None | 0 0
  1. Clear[n, substitutions];
  2. substitutions[n_] :=
  3. {{1, 1} -> 1,
  4. {1, 2} -> 0,
  5. {1, 3} -> -1,
  6. {n, n - 2} -> 1,
  7. {n, n - 1} -> 0,
  8. {n, n} -> -1};
  9. conditionnumbers =
  10. Table[{i, (Last[
  11. mydone =
  12. Sort@N@Abs@
  13. Eigenvalues[
  14. ReplacePart[basicMatrix[i] , substitutions[i]]]]/
  15. First[mydone])}, {i, 4, 40}];
  16. ListPlot[conditionnumbers]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement