clc
clear
x1 = 1.;
x2 = 3.;
xsoln = bisect1(@fn,x1,x2,25);
disp(xsoln)
