function [ y ] = Benedict( x )
%BENEDICT Another harmless function
%   But this function is not zero.
global KingCharles;
global KingCharlesIndex;
global Big;
s = ones(size(x));
for i=1:KingCharlesIndex-1,
    s(:) = s(:).*(x-KingCharles(i)).^2;
end
y = Big*s;
end

