close all
clear
clc

% deklaracja funkcji
f = @(x) 0.00015.*x.^5 - 0.001.*(x - 3).^4 + 0.1.*x.^3 - 5.*(x - 3) - 17;
X = linspace(-10, 10, 1000);
szum = randn(1, 1000);
hold on
plot(X, f(X), 'r');
plot(X, f(X) + szum);
for i=2:6
    p = polyfit(X, f(X) + szum, i);
    plot(X, polyval(p, X));
end
hold off
legend({'f(x)', ...
    'f(x) + szum', ...
    'wielomian stopnia: 2', ...
    'wielomian stopnia: 3', ...
    'wielomian stopnia: 4', ...
    'wielomian stopnia: 5', ...
    'wielomian stopnia: 6', ...
    }, 'Location', 'northwest')
xlabel('x')
ylabel('y')