I need help with my MATLAB code. If I have multiples files of data, for example vBody$1,$ vBody$2,$ vRef$1,$ vRef$2,$ I don't want to copy and paste this code again and again. Can you help me create a function of this code.

$\%$ Initializing arrays

z $=$ zeros$(3,1)$;

B$\_$mat $=$ zeros$(3)$;

no$\_$stars $=$ length$($vBody$)$;

$\%$ Creating for loop to find B and z of the Davenport's q$-$method

for i $=1$:no$\_$stars

$\%$ Extracting the vectors

b $=$ vBody$($:$,$ i$)$;

r $=$ vRef$($:$,$ i$)$;

$\%$ Creating the B matrix

B$\_$mat $=$ B$\_$mat $+$ b $*$ r$\text{'}$;

$\%$ Creating the z vector

z $=$ z $+$ cross$($b$,$ r$)$;

end

$\%$ Calculating the K matrix

K$11=$ B$\_$mat $+$ B$\_$mat' $-$ trace$($B$\_$mat$)*$ eye$(3)$;

K$\_$mat $=[$K$11$ z ; z$\text{'}$ trace$($B$\_$mat$)]$;

$\%$ Finding the eigenvalues and q vectors

$[$q$\_$vec, eigs$]=$ eig$($K$\_$mat $,$'vector'$)$;

$\%$ Finding the biggest eigenvalue and its q vector

$[$lambda $,$ i$]=$ max$($eigs$)$;

q $=$ q$\_$vec$($:$,$ i$)$;

q$\_$hat $=$ q $/$ norm$($q$)$;

$\%$ Calculating the DCM using the function

A$\_$davenport $=$ q$\_$DCM$($q$)$;

$\%$ Confirm that lambda$\_$dave is close to lambda$0$

lambda$\_$init $=$ no$\_$stars;

diff$\_$lambda $=$ lambda$\_$init $-$ lambda;

function skew$\_$mat $=$ Skew$($v$)$

$\%$ This function creates the skew symmetric matrix of a vector

skew$\_$mat $=[0-$v$(3)$ v$(2)$; v$(3)0-$v$(1)$; $-$v$(2)$ v$(1)0]$;

end

function A $=$ q$\_$DCM$($q$)$

$\%$ This function creates the rotation matrix by multiplying two matrices

$\%$ Calculating psi and $($don$\text{'}$t know what it is called$)$

psi $=[($q$(4)*$ eye$(3)-$ Skew$($q$(1$:$3)))$; $-$q$(1$:$3)\text{'}]$;

greek$\_$letter $=[($q$(4)*$ eye$(3)+$ Skew$($q$(1$:$3)))$; $-$q$(1$:$3)\text{'}]$;

$\%$ DCM

A $=$ greek$\_$letter' $*$ psi;

end