function phix=phiuv(zet)
%modified 3/19/98 to allow vector input/output by Peter Guest 

%stable cases
	x=1+4.7*zet;
   phix=x;
	

%unstable cases
        uns=find(zet<0 &zet>-.2);    % &zet>-.2
        if ~isempty(uns);
            y=1./(1-15.*zet(uns)).^.25;
            phix(uns)=y;
           
        end; 
		  unx=find(zet<=-.2);
			 if ~isempty(unx);
           yy=1./(1.26-8.38.*zet(unx)).^.33;
            phix(unx)=yy;
           
        end; 
     
