Metodo de la Secante

%Metodo de la Secante
function y=Secante(FuncionS,Xi,Xj,errorp)
c=1;
fXi=feval(FuncionS,Xi);
fXj=feval(FuncionS,Xj);
fprintf(' i Xi f(Xi) |ep| \n')
fprintf('---------------------------------------------------\n')
%imprimir primera iteracion ep=0
ep=0;
fprintf('%5d %10.8f %10.8f %10.8f \n',c,Xi,fXi,ep);
%imprimir segunda iteracion ep=100
c=c+1;
ep=100;
fprintf('%5d %10.8f %10.8f %10.8f \n',c,Xj,fXj,ep);
while(ep>errorp)
c=c+1;
%guardamos en segundo valor
aux=Xj;
%Calculamos el sgte valor
Xj=Xj-((fXj*(Xi-Xj))/(fXi-fXj));
Xi=aux;
%calcular el error porcentual
ep=abs((Xj-Xi)/Xj)*100;
%calcular f(Xj) y f(Xi)
fXi=fXj;
fXj=feval(FuncionS,Xj);
fprintf('%5d %10.8f %10.8f %10.8f \n',c,Xj,fXj,ep);
end