
function psi=psi_eps(zet)
if zet<0
	
	psi=(1/(1-20*zet).^.333-zet-1/(7-zet));                                               
else
	psi=(1+2.5*zet.^.667).^1.5;

end;