clear
fplot('sqrt(x)',[0,17]);
hold on
grid on
fplot('-1/1008*x^4+77/2520*x^3-301/1008*x^2+1599/1260*x',[0,17])
x=[0 1 4 9 16]; %lähteandmed
y=[0 1 2 3 4];

%lähendamine kuupsplainiga
%teeme võrgu
xvork=0:0.0001:16; %NB!~semikoolon
%esimene element x-i esimene väärtus, viimane element on x-i viimane
%väärtus ja sammupikkuse valime ise
ysplain=interp1(x,y,xvork,'spline'); %NB!~semikoolon

plot(x,y,'r*',xvork,ysplain)
legend('sqrt(x)','interpolant','algandmed','kuupsplain')






