function y = interpol(x, n, xt, yt)
%
% Interpolación lineal.
%
% Suponemos que xt(i) < xt(i+1), i = 1, ..., n-1
%
%
  if ((x >= xt(1)) && (x <= xt(n)))
    jplus1 = 2;
    while (x > xt(jplus1))
      jplus1 = jplus1 + 1;
    end
    j = jplus1 - 1;
    y = yt(j) * (x - xt(jplus1)) / (xt(j) - xt(jplus1)) + ...
        yt(jplus1) * (x - xt(j)) / (xt(jplus1) - xt(j));
  else
    y = 0;
  end
end
%
%       fprintf('It: j+1=%d, x=%g, x(j+1)=%g\n', jplus1, x, xt(jplus1));
