for i=1:100;
   zet=-i/10;
   zt(i)=-zet;
   gt(i)=-2/3/(1-1/(8*zet));
	ge(i)=(20/3/(1-20*zet)^(4/3)-1-1/(7-zet)^2)*zet/psi_eps(zet);
   bet(i)=1+.5*gt(i)-2/3*ge(i);
   jam(i)=(1+1.5*gt(i))/bet(i);
   jax(i)=-.5*(gt(i)+2/3*ge(i))/bet(i);  
   jar(i)=(1-ge(i))/bet(i);
end;

plot(zt(1:100),jax(1:100),zt(1:100),jam(1:100), zt(1:100),jar(1:100));
   
