function VTB2_5(m,m0,e,z,rmin,rmax)

%VTB2_5 Rotating unbalance.
% VTB2_5(m,m0,e,zeta,rmin,rmax) plots the displacement and phase of 
% a system with rotating unbalance between the frequency 
% ratios rmin and rmax.  The mass of the system is m and the damping
% ratio is zeta.  The paramaters of the unbalance are the mass m0 and
% eccentricity e.

r=rmin:(rmax-rmin)/1000:rmax;

Xn=sqrt((r.^4)./((1-r.^2).^2+(2*z*r).^2));%(2.51)

X=Xn.*m0*e/m;%(2.51)

phi=atan2(2*z*r,1-r.^2);%(2.52)

plot(r,Xn)

xlabel('Frequency Ratio')

ylabel('Normalized Displacement Magnitude')

title('Normalized Displacement Magnitude versus Frequency Ratio')

text(0,0,'Press Return To Continue','sc')

pause

plot(r,X)

xlabel('Frequency Ratio')

ylabel('Displacement Magnitude')

title('Displacement Magnitude versus Frequency Ratio')

text(0,0,'Press Return To Continue','sc')

pause

plot(r,phi)

xlabel('Frequency Ratio')

ylabel('Phase')

title('Phase versus Frequency Ratio')

text(0,0,'Press Return To Continue','sc')