% step2s.m 
% compute the eigenvalues of the Jacobian 
% simplified version

beta=input('give beta: ');
arr=arrhen(Xbar(2)); 
arrprime=arr*gamma*Xbar(2)^(-2);
A=[-1-arr -Xbar(1)*arrprime; alpha*arr -1-beta+alpha*Xbar(1)*arrprime] 
[vv,lambda]=eig(A)

% compute the corresponding value of Ubar 
Ubar=Xbar(2)-(1-Xbar(2)+alpha*Xbar(1)*arrhen(Xbar(2)))/beta
