function [xn,n] = sekantenverfahren(b, x, x0, x1) 
fneu = auswertung_newton(x,b,x1);
falt = auswertung_newton(x,b,x0); 
xn(1)=x0;
xn(2)=x1;
e=1e-10;                            %Schwelle für Abbruchkriterium
n=2;
nmax=100;

while (abs(xn(n)-xn(n-1)) > e) && n<=nmax
        xn(n+1)= xn(n) - ((xn(n)-xn(n-1)) / (fneu-falt))* fneu;
        falt = fneu;
        fneu = auswertung_newton(x,b,xn(n+1));
        n=n+1;
end
xn=xn(n);