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