plotPolynomialFit.m

function output = plotPolynomialFit(ds_x, ds_y, N,confidence, resolution)
        [output, error] = polyfit(ds_x,ds_y,N);
        if confidence
            range = linspace(0, max(ds_x), resolution)
            [outputline, delta] = polyconf(output, range, error, 'predopt', 'curve')
            plot(range, outputline, "-");
            hold on
            plot(range,outputline+delta,"r--",range,outputline-delta,"r--");
            hold off
        else
            [outputline, delta] = polyval(output,ds_x, s1_e);
            plot(range,outputline,"--");
        end
end