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