Have you tried adding another forloop for the columns?
for j = 1:10000
I dont know what corr() returns, so you might have to adjust your z() output.
Also, your row indices are going to go out of bound. for example, at max i = 982,
you will be looking at rows 982:982+59 for your A and b matrices (which do not exist)