MATLAB scribt showing error

$\%$ Parameters

m $=167\mathrm{.}67$; $\%$ Mass of motorcycle and rider $($kg$)$

x$0=0\mathrm{.}30$; $\%$ Initial displacement $($m$)$

v$0=0\mathrm{.}1$; $\%$ Initial velocity $($m$/$s$)$

$\%$ Mono Shock Suspension

k$1=32\mathrm{.}37$e$3$; $\%$ Stiffness $($N$/$m$)$

zeta$1=0\mathrm{.}23$; $\%$ Damping ratio

$\%$ Gas$-$Filled Suspension

k$2=24\mathrm{.}16$e$3$; $\%$ Stiffness $($N$/$m$)$

zeta$2=0\mathrm{.}17$; $\%$ Damping ratio

$\%$ Calculate natural frequencies

omega$\_$n$1=$ sqrt$($k$1/$ m$)$; $\%$ Mono shock

omega$\_$n$2=$ sqrt$($k$2/$ m$)$; $\%$ Gas filled

$\%$ Calculate damping coefficients

c$1=2*$ zeta$1*$ omega$\_$n$1*$ m;

c$2=2*$ zeta$2*$ omega$\_$n$2*$ m;

$\%$ Define time vector

t $=$ linspace$(0,2,1000)$; $\%$ Simulation from $0$ to $2$ seconds

$\%$ Define transfer functions for each suspension system

sys$1=$ tf$([0,0,$ k$1],[$m$,$ c$1,$ k$1])$;

sys$2=$ tf$([0,0,$ k$2],[$m$,$ c$2,$ k$2])$;

$\%$ Simulate the system response with initial conditions

$[$response$1,$ t$1]=$ initial$($sys$1,[$x$0$; v$0],$ t$)$;

$[$response$2,$ t$2]=$ initial$($sys$2,[$x$0$; v$0],$ t$)$;

$\%$ Plot responses

figure;

plot$($t$1,$ response$1,\text{'}$b$\text{'},$ 'LineWidth', $1\mathrm{.}5)$; hold on;

plot$($t$2,$ response$2,\text{'}$r$\text{'},$ 'LineWidth', $1\mathrm{.}5)$;

title$(\text{'}$Displacement Response of Mass for Different Suspensions'$)$;

xlabel$(\text{'}$Time $($s$)\text{'})$;

ylabel$(\text{'}$Displacement $($m$)\text{'})$;

legend$(\text{'}$Mono Shock Suspension', 'Gas$-$Filled Suspension'$)$;

grid on;

$\%$ Decision Analysis

disp$(\text{'}$Analysis based on maximum displacement observed:$\text{'})$;

if max$($abs$($response$1))$ max$($abs$($response$2))$

disp$(\text{'}$Mono Shock Suspension provides better control over displacement.$\text{'})$;

else

disp$(\text{'}$Gas$-$Filled Suspension provides better control over displacement.$\text{'})$;

end