SAS IML求矩阵相关系数程序有问题,

SAS IML求矩阵相关系数程序有问题,
proc iml ;
A={0.483 13.2682 0 4.3646 5.1070,
0.4035 13.4909 39.0131 3.6151 5.5005,
0.8979 25.7776 9.0513 4.8920 7.5342,
0.5927 16.0245 13.2935 4.4529 6.5913
};
r=nrow(A);
c=ncol(A);
B=A;
do i=1 to c;
do p=1 to r;
if i=3
then B[p,i]=(A[,i]-A[p,i])/(A[,i]-A[>