creating multidiagonal matrix from a basic square matrix

For a finite difference code I need to create a multidiagonal matrix from a basic square matrix A:

How can I do this without using a for loop?

Answer 1

[Andrei Bobrov]

A = magic(3);
n = 8;

B = repmat({A},1,n);
out = blkdiag(B{:});

[Saeid]

Thanks Andrei. That was great help!