function y = interpolc(x, n, xt, yt)
%
% Interpolación lineal.
%
% Suponemos que xt(i) < xt(i+1), i = 1, ..., n-1
%
%

global NroInvocaciones NroIteraciones;

NroInvocaciones = NroInvocaciones + 1;

  if ((x >= xt(1)) && (x <= xt(n)))
    jplus1 = 2;
    while (x > xt(jplus1))
      NroIteraciones = NroIteraciones + 1;
      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
%
