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

end;