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