function y = FcnEval(x)
y = zeros(size(x));
for ii = 1:numel(x)
if x(ii) < -2;
y(ii) = 1./((x(ii).^2) + 1);
elseif -2 <= x(ii) && x(ii) < -1;
y(ii) = sin(x(ii));
elseif -1 <= x(ii) && x(ii) < 1;
y(ii) = x(ii)^2;
elseif 1 <= x(ii) && x(ii) < 2;
y(ii) = 2;
else
y(ii) = 1./(x(ii)+1);
end
end
or
y = zeros(size(x));
t1 = x < -2;
y(t1) = 1./(x(t1).^2 + 1);
t2 = -2 <= x & x < -1;
y(t2) = sin(x(t2));
t3 = -1 <= x & x < 1;
y(t3) = x(t3).^2;
y(1 <= x & x < 2) = 2;
t4 = x >= 2;
y(t4) = 1./(x(t4)+1);
or
f = {@(x)1./(x.^2 + 1);@sin;@(x)x.^2;@(x)2;@(x)1./(x+1)};
[b,b]=histc(x,[-inf;-2;-1;1;2;inf]);
y = arrayfun(@(ii,jj)f{ii}(jj),b,x);
