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

end;