clc
clear
x1 = 2.;
x2 = 3.;
xsoln = bisect(@fn,x1,x2);
disp(xsoln)
