load polinom.mat

figure

% 10 pontos ------------------------------------------
predikaltTartomany10 = linspace(-4,8, 10);

hatodfoku = polyfit(meresiPozicio, mertErtekek, 6);

predikaltErtekek10 = polyval(hatodfoku, predikaltTartomany10);
subplot(1,2,1)
plot(predikaltTartomany10, predikaltErtekek10, 'r')

hold on

% 100 pontos ------------------------------------------
predikaltTartomany100 = linspace(-4,8, 100);

hatodfoku = polyfit(meresiPozicio, mertErtekek, 6);

predikaltErtekek100 = polyval(hatodfoku, predikaltTartomany100);
plot(predikaltTartomany100, predikaltErtekek100, 'b')

% Derivált
dErtek10 = diff (predikaltErtekek10);

dErtek100 = diff (predikaltErtekek100);


subplot(1,2,2)
plot(predikaltTartomany10(1:length(dErtek10)), dErtek10,'r')
hold on
plot(predikaltTartomany100(1:length(dErtek100)), dErtek100,'b')