format long;

n=2048 
h=1:h:2;
h=1/n;
x=1:h:2;
y=abs( (x-sqrt(2)) ).*(x-sqrt(2)).^2;
hstar = h/100;
xstar= 1:hstar:2;
ystar = spline(x,y,xstar);
max( abs(ystar-abs((xstar-sqrt(2))).*(xstar-sqrt(2)).^2 ) )