vastus225=1/60*2.25^3-1/4*2.25^2+37/30*2.25
fplot(@(x)[sqrt(x),1/60*x.^3-1/4*x.^2+37/30*x],[0,10])
legend('sqrt(x)','interpoletsioonipolünoom')