function Mar27()
% L(x) = (1/2) + (2/3)*(x-.25) is the interpolant
% e(x) = L(x) - sqrt(x). You can find the max err by calculus or
x = linspace(.25,1,10000);
maxErr = max(sqrt(x)-(.5 + (x-.25)*(2/3)))