
function psi=psi_ft(zet)
	
	psi=5.5./(1-7.0*zet).^.667;                                               
   uns=find(zet>0);
if ~isempty(uns)

	psi(uns)=5.5*(1+2.5*zet(uns).^.667);

end;