function rbc6_steady=f1(x);
% xの中身はbeta, delta, psiの順番
beta=x(1) ; delta=x(2) ; psi=x(3) ;  

% グローバル変数の定義
global alpha lam rSS IYSS LSS  KSS CSS ISS YSS Y_LSS ;

% 定常状態の条件式、ノートを参照
SS1=rSS-(1/beta-1) ;
x=IYSS/delta ; %x=(K/L)^(1-alpha), ノート参照
SS2=alpha/x-delta-rSS ;
SS3=psi-(1-alpha)*(LSS^(-1-lam))/(1-IYSS) ;

% ついでに諸変数の定常値を計算しておく
K_LSS=x^(1/(1-alpha)) ; %K/L
KSS=K_LSS*LSS ; %K
YSS=(KSS^alpha)*(LSS^(1-alpha)) ; %Y
CSS=(1-IYSS)*YSS ;
ISS=IYSS*YSS ;
Y_LSS=YSS/LSS ;

rbc6_steady  = [SS1 ; SS2 ; SS3];    
