clear all;
clc;

s= tf('s')  % Setting 's' as a frequency domain variable

D = 1;

G = (s+D)/((s^2+D*s+9.8))

num = [1 D];
den = [1 D 9.8];

Gs = tf(num,den)

[r,p,k] = residue(num,den);

step(G)

hold on;
impulse(G)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% hold on
% 
% D = 2;
% 
% G = (s+D)/((s^2+D*s+9.8))
% 
% step(G)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% hold on
% 
% D = 6;
% 
% G = (s+D)/((s^2+D*s+9.8))
% 
% step(G)