function z = damp(xi)
% damp  Evaluate f(xi) function for optimum damping ratio.  Problem 6-21
z = cos( 4*xi.*sqrt(1-xi.^2) ) + 1 - 8*xi.^2 + 8*xi.^4;