function out = kysu2(y)

out = sqrt((1+0.5.*y)./(1-0.5.*y));