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

x=[-5:0.001:5];        %przedział

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

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

while(stop==1)
    i=i+1;
    x0=(x1+x2)/2;
    y0=8+x0*2+x0^2+4*x0^3;
    y1=8+x1*2+x1^2+4*x1^3;
    
        if((y1*y0)<0)
            x2=x0;
        else
            x1=x0;
        end
    plot(x0,y0,'ys');
    delta=(x2-x1)/x2;
    wzgl=abs(delta);
        if(wzgl<eps)
            stop=0;
        else
            stop=1;
        end
        
    
end
fprintf('Miejsce zerowe funkcji to: %d', x0, s);