I need help with my MATLAB code. I want to create a function of the following code, so if I have multiple S and z values.

S $=$ B$\_$mat $+$ B$\_$mat';

kappa $=$ trace$($adjoint$($S$))$;

$\%$ Characteristic Equation $($Eqn$.3\mathrm{.}73)$

f $=$ @$($x$)($x$^2-$ trace$($B$\_$mat$)^2+$ kappa$)*($x$^2-$ trace$($B$\_$mat$)^2-...$

dot$($z$,$ z$))-($x $-$ trace$($B$\_$mat$))*($z$\text{'}*$S$*$z $+$ det$($S$))-...$

z$\text{'}*($S$\text{'}*$S$)*$z;

$\%$ Derivative of the characteristic equation

df $=$@$($x$)2*$x$*($x$^2-$ trace$($B$\_$mat$)^2-$ dot$($z $,$ z$))+...$

$2*$x$*($x$^2-$ trace$($B$\_$mat$)^2+$ kappa$)+($z$\text{'}*$S$*$z $+$ det$($S$))$ ;

$\%$ Calculate the eigenvalue Newton$-$Raphson iteration

lambdas $=$ zeros$(4,1)$;

lambdas$(1)=10$;

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

quat $=$ zeros $(4,3)$;

DCM $=$ zeros $(3,3,3)$;

for i $=1$:$3$

Xn $=$ f$($lambda$\_$init$)/$ df$($lambda$\_$init$)$;

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

$\%$ get the corresponding attitude estimates

rho $=$ lambda$\_$init $+$ trace$($B$\_$mat$)$;

rho$\_$mat $=$ rho $*$ eye $(3)-$ S;

q$\_$unit $=[$adjoint$($rho$\_$mat$)*$ z ; det$($rho$\_$mat$)]$;

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

$\%$ store the step size and attitude estimate

step$($i$)=$ Xn;

lambdas$($i $+1)=$ lambda$\_$init;

quat$($:$,$ i$)=$ q$\_$unit;

DCM$($:$,$ :$,$ i$)=$ q$\_$DCM$($q$\_$unit$)$;

$\%$ Capture lambda$\_$init and DCM after first and third iterations

if i $==1$

lambda$\_$after$\_$first $=$ lambda$\_$init;

DCM$\_$after$\_$first $=$ DCM$($:$,$ :$,$ i$)$;

elseif i $==3$

lambda$\_$after$\_$third $=$ lambda$\_$init;

DCM$\_$after$\_$third $=$ DCM$($:$,$ :$,$ i$)$;

end

end

$\%$ Compare lambda$\_$init after the first and third iterations

lambda$\_$difference $=$ abs$($lambda$\_$after$\_$third $-$ lambda$\_$after$\_$first$)$;

$\%$ Compare DCM after the first and third iterations

DCM$\_$difference $=$ abs$($DCM$\_$after$\_$third $-$ DCM$\_$after$\_$first$)$;