$\%$ Define the natural frequency $($wn$)$ and the range of damping factors $($zeta$)$wn $=1$; $\%$ Normalized natural frequencyzeta$\_$values $=[0\mathrm{.}1,0\mathrm{.}2,0\mathrm{.}3,0\mathrm{.}5,0\mathrm{.}7,1,1\mathrm{.}5]$; $\%$ Different damping factors$\%$ Frequency range for the Bode plotw $=$ logspace$(-1,1,500)$; $\%$ Frequency range from $0\mathrm{.}1$ to $10$ in log scale$\%$ Loop over different zeta values and plot Bode plotfigure;hold on;for zeta $=$ zeta$\_$values $\%$ Define the transfer function for each zeta num $=[100]$; $\%$ Numerator s$^2$ term den $=[12*$zeta$*$wn wn$^2]$; $\%$ Denominator s$^2+2*$zeta$*$wn$*$s $+$ wn$^2$ sys $=$ tf$($num$,$ den$)$; $\%$ Create transfer function $[$mag$,$ phase$]=$ bode$($sys$,$ w$)$; $\%$ Compute magnitude and phase $\%$ Convert magnitude to dB mag$\_$dB $=20*$log$10($squeeze$($mag$))$; phase$\_$deg $=$ squeeze$($phase$)$; $\%$ Convert phase to degrees $\%$ Plot magnitude and phase subplot$(2,1,1)$; semilogx$($w$/$wn$,$ mag$\_$dB$,$ 'DisplayName', $[\text{'}\backslash$zeta $=\text{'}$ num$2$str$($zeta$)])$; xlabel$(\text{'}\backslash$omega $/\backslash$omega$\_$n$\text{'})$; ylabel$(\text{'}20$ log $|$G$($j$\backslash$omega$)|\text{'})$; title$(\text{'}$Magnitude Response'$)$; grid on; hold on; subplot$(2,1,2)$; semilogx$($w$/$wn$,$ phase$\_$deg, 'DisplayName', $[\text{'}\backslash$zeta $=\text{'}$ num$2$str$($zeta$)])$; xlabel$(\text{'}\backslash$omega $/\backslash$omega$\_$n$\text{'})$; ylabel$(\text{'}$Phase $($degrees$)\text{'})$; title$(\text{'}$Phase Response'$)$; grid on; hold on;end$\%$ Add legendssubplot$(2,1,1)$;legend$(\text{'}$show$\text{'})$;subplot$(2,1,2)$;legend$(\text{'}$show$\text{'})$;

Make a graph in code above symmetrical to the y$-$axis