function [V]=bisect(a,b,m)
% This assumes that the function is assumed to 
%  be set to my_function

u = my_function(a); v = my_function(b); 

[a b], [u v] 

V=[];e=b-a;

for n=1:m+1;
    c = (a+b)/2; e=e/2;
    w = my_function(c);
    V=[V; [n c e w]];
    if w*v == 0 
       return
    elseif (w*v < 0) 
        a=c;u=w;
     else
        b=c;v=w;
    end
end
end
