
%figure 1;




%figure 2;

step = .01;
f = @(K)  bisectCool(1,5,.0000000001,K);
x = 1:step:5;
y = arrayfun(f,x);

plot(x,y);


%for n = 1:(5-1)/step
%    y(n) = f(x(n));
%end