function y = myrat(K,phi)
%yF = fnorm(K);

yphi = phi*K;

yPsi = psinorm(K);

y = yphi/(yPsi+realmin);

return