function y = f521(x)
y = sqrt(x-2);