function zerofinder()
clear; close all
a=10; b = 3;
x = linspace(-20,20,100);
y = func(x,a,b);
zz = zeros(100)
plot(x,y,x,zz)

%[z,val,iter] = newton(@func,1,.0001,40,a,b);
%z = newton(@func,1,.001,40,a,b)
[zero, res, niter] = bisect(func, -5, 5, .0001, 40, a,b)

function z = func(y,a,b)
z=y.^2+a*y+b;