function y = phinv(x,a,b)
y = (x-a)/(b-a);