function v=wigglyfn(x,y)
% y may be a vector, but x is assumed scalar.

v=zeros(1,length(y));
if abs(x)<10
w=(abs(y)<10);
v(w)=(2+cos(x-y(w)))./(1+.1*(y(w)).^2);
end
    