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

f = my_function(a);
df = der_my_function(a); 

V=[0 a f df];

for n=1:m;
    f  = my_function(a);
    df = der_my_function(a); 
    a  = a-f/df;
    V  = [V; [n a f df]];
  end
end
