%Hope this could help you!
% Define the transfer function
num = [10 120 660 2280 4730 4600 1600];
den = [1 14 113 628 2379 6350 11115 9800 3200];
sys = tf(num, den);
% Plot the unit step response
step(sys);
% Add title and axis labels
title('Unit Step Response');
xlabel('Time (seconds)');
ylabel('Output');
% Customize the axis
xlim([0 10]); % set x-axis limits
ylim([0 1.2]); % set x-axis limits
grid on; % add grid lines
