function y = target(x)
% target   for use with  interp.m and  lsq.m
y = zeros(size(x));
i = find(x<pi);
j = find(x>pi | x==pi);
y(i) = cos(x(i));
y(j) = -exp(pi-x(j));
