
function psi=psi_ft(zet)
if zet<0
	
	psi=5.5/(1-7.0*zet)^.667;                                               
else
	psi=5.5*(1+2.5*zet^.667);

end;