function indicator = re16point2(r,y,x)
% Indicator function for Eq. 16.2
indicator = -sin(y.^2) + ...
   (1-x.^r)./(1 + x.^r)*cos(r.*sqrt(x)./y);