function [ N ] = N(h) N=((2+h)./(2-h)).^(1./h); end