jared = @(x) (exp(x))+5*x^3-5*x^2+6*x-10;
djared = @(x) (exp(x))+15*x^2-10*x+6;

i = input('Initial value: ?');
low = input('Lower bound: ?');
up = input('Upper bound: ?');
err = input('Error: ?');

bisection2(jared,low,up,err)
newton(i, err, jared, djared)
fzero(jared,i)