Can somebody help me what is the error in this code?

$\%$Define the symbolic Variable x

syms x;

$\%$Enter the upper and lower functions as y$1($x$)$ and y$2($x$),$ respectively.

k$=$randi$(20,1)$;

y$1($x$)=$ k$-$x$^2$;

y$2($x$)=$ x$^2+$x$-$k;

$\%$Solve the point of intersections. Save the array of solutions as roots.

roots $=$ solve$($y$1==$ y$2,$ x$)$;

$\%$Save the lower value as x$1$ and the higher limit as x$2.$ Use min and max to compare the roots for lower and upper limit$.$

$\%$Convert the roots from symbolic constants to double.

x$1=$ min$($roots$)$;

x$2=$ max$($roots$)$;

x$1=$ double$($x$1)$;

x$2=$ double$($x$2)$;

$\%$Find the points of intersection between the twop curves. Set it as an array P$1=[$x$,$ y$]$

P$1=[$x$1,$ y$1($x$1)]$;

P$2=[$x$2,$ y$2($x$2)]$;

$\%$Set the difference of the upper minus the lower function as I. Use absolute value function to ensure positive area

I $=$ abs$($y$1-$ y$2)$;

$\%$Use the integration to the find the area bounded by the two curves at the computed boundary.

Area $=$ int$($I$,$ x$,$ x$1,$ x$2)$;

$\%$Determine the moment along y axis.

My $=$ int$($x $*$ I, x$,$ x$1,$ x$2)$;

$\%$Determine the centroidal element ybar by finding the midoint of the lower and upper curves.

ybar $=$ double$($My $/$ Area$)$;

$\%$Determine the moment along y axis.

Mx $=$int$($ybar$*$I$,$ x$,$ x$1,$ x$2)$;

$\%$Find the centroidal elements xc$,$yc by dividing the Moments by the Area. Convert the answer as Doubles

xc $=$ double$($Mx $/$ Area$)$;

yc $=$ double$($ybar$)$;

$\%$Combine the centroidal elements as Centroid as an array of element xc and yc

Centroid $=[$xc$,$ yc$]$;

$\%$Set the graphing margin allowance to $1$

margin $=1$;

$\%$Graphing of the regions

ezplot$($y$1,[$x$1-$margin, x$2+$margin$])$;

hold on;

ezplot$($y$2,[$x$1-$margin,x$2+$margin$])$;

grid on;

k$=$ linspace$($x$1,$x$2)$;

k$2=[$k$,$fliplr$($k$)]$;

inBetween $=[$y$1($k$),$ fliplr$($y$2($k$))]$;

fill$($k$2,$ inBetween, $\text{'}$y$\text{'})$;

plot$($x$1,$y$1($x$1),"$r$*")$;

plot$($x$2,$y$1($x$2),"$r$*")$;

title $("$Centroid of Plane Region"$)$

plot$($xc$,$yc$,"$bo$")$