function phi=phiuk(zet)
if zet<0
  	phi=1/(1-16*zet)^.25;
else
	phi=1+4.7*zet;
	
end;