% function whose zero is to be found by bisection
% function val = func(x)
function val = func(x)
% coef contains coefficients of the polynomial (x-2)^13, multiplied out, i.e.
% (x-2)^13 = x^13 - 26*x^12 + 312*x^11 - 2288*x^10 + ... - 8192
coef=[1 -26 312 -2288 11440 -41184 109824 -219648 329472 -366080 292864 ...
     -159744 53248 -8192];
% Evaluate polynomial (x-2)^13 using coef
val = polyval(coef,x);
%
% Evaluate polynomial (x-2)^13 using factored form
% val = (x-2).^13;
%
% val = sin(x);