function vtb5_1(m,c,k)

%VTB5_1 Transmissibility ratio of a SDOF system.
%  VTB5_1(m,c,k) plots the tranmissibility ratio for the system
%  described by the mass m, damping coefficient c, and stiffness k.
%  VTB5_1(zeta,w) plots the transmissibility ratio for the system
%  described by the damping ratio zeta and undamped natural
%  frequency w.



if nargin==2

   z=m;

   w=c;

else

   w=sqrt(k/m);

   z=c/(2*m*w);

end

r=[0:.01:2]';



%Calculate various parts of equation 5.7 for TR



TR1=ones(201,1)+(2*z.*r).^2;



TR2=(ones(201,1)-r.^2).^2+(2*z.*r).^2;



TR=sqrt(TR1./TR2);



plot(r,TR)



title(['Transmissibility plot for zeta = ',num2str(z),' w = ',num2str(w),' rad/s'])





grid

ylabel('Transmissibility Ratio')

xlabel('Dimensionless Frequency')