fplot(@(x)[sqrt(x),-1/6*x^2+7/6*x],[0,5])
legend('sqrt(x)','-1/6*x^2+7/6*x')
grid on