fplot(@(x)1/12*x.^3-13/12*x.^2+3.5*x,[0,7])