%Lähteandmed
x=[1 2 3 4];
y=[0 8 5 6];
plot(x,y,'r*') %tähistame algandmedmustade tärnidega
hold on
clear
fplot(@(x)5/2*x^3-41/2*x^2+52*x-34,[0,5],'g') %piirkonna valime ise sobivalt
legend('lähteandmed','interpolatsioonipolünoom')