Problem with for loops

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 all
clear
clc

[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);

end

f5(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);

figure
plot(f5)
covariance5=cov(f5');
figure
surf(X,Y,covariance5)
figure
surf(X,Y,Z)

Answer 1

[Geoff Hayes]

Jaime - please don't dynamically create variables. See Stephen's post at https://www.mathworks.com/matlabcentral/answers/304528-tutorial-why-variables-should-not-be-named-dynamically-eval.

In your code, you are creating ten arrays of long (should rename this variable to make it more clear what it represents). This is then equivalent to a single 10 x m array (where m is the length of your X) which should be just

 v_a = normrnd(0,1,10,length(X));

Then, in your for loops, you can just reference the appropriate column of v_a instead of trying to figure out which array variable you will need.

[Jaime de la Mota]

This helps me. Thanks.

[Torsten]

And the reason that results are not reproducable is because you use normally distributed random numbers which will change from one execution of the code to the next.