y=0;
stop=1;
eps=0.001;

x=[-5:0.001:5];      

y=8+x.*2+x.^2+4*x.^3;   % fjunckja
plot(x,y,'-r'); hold on;

x1=-5;
x2=5;
i=0;

while(stop==1)
    i=i+1;
    x0=(x1+x2)/2;           % ustalanie polowy przedzialu
    y0=8+x0*2+x0^2+4*x0^3;  %
    y1=8+x1*2+x1^2+4*x1^3;  % pierwsza pochodna
    
        if((y1*y0)<0)       %warunek wybrania przedziału
            x2=x0;
        else
            x1=x0;
        end
    plot(x0,y0,'o');
    delta=(x2-x1)/x2;
    wzgl=abs(delta); % moduł
        if(wzgl<eps)
            stop=0;
        else
            stop=1;
        end
        
    
end
fprintf('miejsce zerowe: %d', x0, i);