function y = fp(x)

y = 1/2*x^(-1/2);