x=[0 1 4 9]
y=[0 1 2 3]
plot(x,y,'r*')
hold on
fplot(@(x)1/60*x^3-1/4*x^2+37/30*x,[0,9])
fplot(@(x)sqrt(x),[0,9])
legend('algandmed','interpolant','tegelik funktsioon')