function x = nrsecsolve(F, a, b)
%
% Resolución de ecuación algebraica por el método Newton-Raphson
% secante.
%
% Constantes
epsilon = 1e-12;
delta   = 1e-6;
%
x  = a;
xp = b;
%
fxp = F(xp);
%
while (abs(x - xp) >= epsilon)
%
  fx = F(x);
  Df = fx - fxp;
  if (abs(Df) > delta)
    DfDx = Df / (x - xp);
  end
  %fprintf('x=%16.14f F(x)=%16.14f DFDx=%16.14f\n', x, fx, DfDx);
  %
  xp  = x;
  fxp = fx;
  x   = x - fx / DfDx;
end
%
end

