%composite newtons 3/8 rule
function [gs] = newton(f,a,b,n)
    h = (b-a)/n;
    gs = cell(n,1);
    counter = 1;
    for i = a:h:(b - h)
        %i = a, i+h = b
        gs{counter} = @(y) ((h)/8 * (f(i, y) + 3 * f(((2*i)+(i+h))/3,y) + 3 * f((i+(2*(i+h)))/3,y) + f(i+h,y)));
        counter = counter + 1;
    end
end