Advertisement
Guest User

Untitled

a guest
Oct 24th, 2014
131
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.36 KB | None | 0 0
  1. Big[A0_, m0_, a0_, b0_, c0_, d0_] :=
  2.  
  3. Module[{A = A0, m = m0, a = a0, b = b0, c = c0, d = d0},
  4.  
  5. For[i = 1, i <= 64, i = i + 4,
  6. For[j = 1, j <= 64, j = j + 4,
  7. A[[i ;; i + 3, j ;; j + 3]] =
  8. DiagonalMatrix[(m + 1)*
  9. Diagonal[A[[i ;; i + 3, j ;; j + 3]]] + {a, b, c, d}*
  10. RotateLeft[Diagonal[A[[i ;; i + 3, j ;; j + 3]]], 1]]
  11. ]
  12. ];
  13. A]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement