function y=f(x)

y = sqrt( (x-10).^2 + (x.^2+10000-1000).^2 );
