for j=1:m-1 for k=j+1:n A(j,k)=A(j,k); A(j,j)=A(j,j); A(k,j)=A(k,j)/A(j,j); for p=j+1:m for q=k:n A(p,q)=A(p,q)-(A(q,j)*A(j,p)); end end end end display('aaaaaa'); A what is the problem with this,

2.0000         0    4.0000
     0    1.0000    2.0000
0.5000    4.0000   -8.0000