function y=bisection(a,b)
fprintf('
for i=1:5
if (f(a)*f(b))<0
    mid=(a+b)/2
    if f(a)*f(mid)<0
        b=mid;
    else
        a=mid;
    end
else 
    disp ('Choose another interval')
end
end
end

function ret=f(x)
ret=x.^2-3;
end