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

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

function z = transc(y)
         z=y.^3-exp(y);
