Hello everyone. I have coded a function, however, I cannot reproduce the same result using for loops. Could someone please tell me how to do it?
Thanks.
Here it is the code in question:
close allclearclc[X,Y] = meshgrid(-1:0.5:1,-1:0.5:1);Z=exp((-1)*abs(X-Y));long=length(X);[eigenvectorsZ, eigenvalues_matZ]=eig(Z);eigenvaluesZ=eig(Z);for i=1:long v_a(i)=normrnd(0,1); v_a2(i)=normrnd(0,1); v_a3(i)=normrnd(0,1); v_a4(i)=normrnd(0,1); v_a5(i)=normrnd(0,1); v_a6(i)=normrnd(0,1); v_a7(i)=normrnd(0,1); v_a8(i)=normrnd(0,1); v_a9(i)=normrnd(0,1); v_a10(i)=normrnd(0,1);endf5(1,1)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a(1);f5(2,1)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a(1);f5(3,1)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a(1);f5(4,1)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a(1);f5(5,1)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a(1); f5(1,2)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a2(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a2(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a2(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a2(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a2(1);f5(2,2)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a2(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a2(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a2(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a2(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a2(1);f5(3,2)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a2(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a2(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a2(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a2(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a2(1);f5(4,2)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a2(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a2(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a2(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a2(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a2(1);f5(5,2)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a2(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a2(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a2(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a2(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a2(1);f5(1,3)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a3(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a3(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a3(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a3(1);f5(2,3)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a3(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a3(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a3(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a3(1);f5(3,3)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a3(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a3(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a3(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a3(1);f5(4,3)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a3(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a3(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a3(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a3(1);f5(5,3)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a3(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a3(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a3(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a3(1);f5(1,4)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a4(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a4(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a4(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a4(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a4(1);f5(2,4)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a4(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a4(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a4(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a4(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a4(1);f5(3,4)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a4(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a4(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a4(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a4(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a4(1);f5(4,4)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a4(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a4(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a4(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a4(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a4(1);f5(5,4)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a4(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a4(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a4(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a4(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a4(1);f5(1,5)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a5(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a5(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a5(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a5(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a5(1);f5(2,5)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a5(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a5(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a5(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a5(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a5(1);f5(3,5)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a5(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a5(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a5(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a5(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a5(1);f5(4,5)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a5(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a5(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a5(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a5(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a5(1);f5(5,5)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a5(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a5(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a5(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a5(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a5(1);f5(1,6)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a6(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a6(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a6(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a6(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a6(1);f5(2,6)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a6(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a6(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a6(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a6(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a6(1);f5(3,6)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a6(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a6(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a6(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a6(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a6(1);f5(4,6)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a6(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a6(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a6(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a6(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a6(1);f5(5,6)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a6(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a6(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a6(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a6(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a6(1);f5(1,7)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a7(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a7(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a7(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a7(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a7(1);f5(2,7)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a7(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a7(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a7(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a7(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a7(1);f5(3,7)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a7(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a7(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a7(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a7(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a7(1);f5(4,7)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a7(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a7(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a7(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a7(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a7(1);f5(5,7)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a7(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a7(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a7(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a7(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a7(1);f5(1,8)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a8(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a8(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a8(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a8(1);f5(2,8)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a8(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a8(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a8(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a8(1);f5(3,8)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a8(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a8(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a8(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a8(1);f5(4,8)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a8(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a8(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a8(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a8(1);f5(5,8)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a8(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a8(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a3(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a8(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a8(1);f5(1,9)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a9(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a9(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a9(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a9(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a9(1);f5(2,9)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a9(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a9(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a9(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a9(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a9(1);f5(3,9)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a9(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a9(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a9(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a9(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a9(1);f5(4,9)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a9(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a9(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a9(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a9(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a9(1);f5(5,9)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a9(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a9(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a9(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a9(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a9(1);f5(1,10)=eigenvectorsZ(1,5)*sqrt(eigenvaluesZ(5))*v_a10(5)+eigenvectorsZ(1,4)*sqrt(eigenvaluesZ(4))*v_a10(4)+eigenvectorsZ(1,3)*sqrt(eigenvaluesZ(3))*v_a10(3)+eigenvectorsZ(1,2)*sqrt(eigenvaluesZ(2))*v_a10(2)+eigenvectorsZ(1,1)*sqrt(eigenvaluesZ(1))*v_a10(1);f5(2,10)=eigenvectorsZ(2,5)*sqrt(eigenvaluesZ(5))*v_a10(5)+eigenvectorsZ(2,4)*sqrt(eigenvaluesZ(4))*v_a10(4)+eigenvectorsZ(2,3)*sqrt(eigenvaluesZ(3))*v_a10(3)+eigenvectorsZ(2,2)*sqrt(eigenvaluesZ(2))*v_a10(2)+eigenvectorsZ(2,1)*sqrt(eigenvaluesZ(1))*v_a10(1);f5(3,10)=eigenvectorsZ(3,5)*sqrt(eigenvaluesZ(5))*v_a10(5)+eigenvectorsZ(3,4)*sqrt(eigenvaluesZ(4))*v_a10(4)+eigenvectorsZ(3,3)*sqrt(eigenvaluesZ(3))*v_a10(3)+eigenvectorsZ(3,2)*sqrt(eigenvaluesZ(2))*v_a10(2)+eigenvectorsZ(3,1)*sqrt(eigenvaluesZ(1))*v_a10(1);f5(4,10)=eigenvectorsZ(4,5)*sqrt(eigenvaluesZ(5))*v_a10(5)+eigenvectorsZ(4,4)*sqrt(eigenvaluesZ(4))*v_a10(4)+eigenvectorsZ(4,3)*sqrt(eigenvaluesZ(3))*v_a10(3)+eigenvectorsZ(4,2)*sqrt(eigenvaluesZ(2))*v_a10(2)+eigenvectorsZ(4,1)*sqrt(eigenvaluesZ(1))*v_a10(1);f5(5,10)=eigenvectorsZ(5,5)*sqrt(eigenvaluesZ(5))*v_a10(5)+eigenvectorsZ(5,4)*sqrt(eigenvaluesZ(4))*v_a10(4)+eigenvectorsZ(5,3)*sqrt(eigenvaluesZ(3))*v_a10(3)+eigenvectorsZ(5,2)*sqrt(eigenvaluesZ(2))*v_a10(2)+eigenvectorsZ(5,1)*sqrt(eigenvaluesZ(1))*v_a10(1);figureplot(f5)covariance5=cov(f5');figuresurf(X,Y,covariance5)figuresurf(X,Y,Z)
Best Answer