bisekcja

function [wyn] = bisekcja(f,a,b)
epsil = 0.001;
fa = f(a);
fb = f(b);
c = (a+b)/2;
fc = f(c);
wyn = [0 0];
if(fa*fb > 0)
    disp('Funkcja nie spelnia zalozen')
else
while(abs(fc)>epsil)
    c = (a + b)/2;
    fc = f(c);
    if(abs(fc)<epsil)
        wyn(1,1) = a;
        wyn(1,2) = b;
        break
    end
    if(fa*fc<0)
        b = c;
    else
        a = c;
        fa = fc;
    end
end
end
end
%funkcja anonimowa i wywołanie
%funk1=@(x)x^3+x^2-1
%bisekcja(funk1,0,2)