function p=pn(x)

p=1/x;
% p=1/(x.^2);