Revision history [back]

Hello every one, I am a user of Sagemath online, I wrote theses below lines to find the steady state of my problem, But at the end I have no result except this []!! I would like to know what this means? It means there is no solution for this system equations? I put my lines here. Thank so much if anyone could help me. var ('Omegabar Qbar ') Ybar=1 pbeta=.99 epsilon=-1 sigma=.03 theta=.3 k=13 delta=.025 rho=1.1 Ybar=1; Lambdabar=pbeta Rlbar=1+1/Qbar Rbar=1/pbeta Rebar=1/pbeta munbar=LambdabarOmegabarRbar musbar=LambdabarOmegabar(Rlbar-Rbar) muebar=LambdabarOmegabar(1-pbeta) epersbar=musbar/muebar xbar=-epersbar+(epersbar^(2)+(2/k)(1-epsilonepersbar))^(1/2) Thetabar=theta(1+epsilonxbar+(k/2)xbar^(2)) phibar=Omegabar/(Thetabar-musbar-xbarmuebar) Ebar=xbarQbar Nbar=Qbar/phibar Dbar=Qbar-Nbar-Ebar Nnbar=delta(Qbar+Ybar)

eq1=Omegabar==(sigma+(1-sigma)(munbar+phibarmusbar+phibarxbarmuebar)) eq2=Nbar==(1-sigma)(Qbar+Ybar-RebarEbar-Rbar*Dbar)+Nnbar

solve([eq1,eq2], Qbar)

(Omegabar, Qbar) []

Hello every one, I am a user of Sagemath online, I wrote theses below lines to find the steady state of my problem, But at the end I have no result except this []!! I would like to know what this means? It means there is no solution for this system equations? I put my lines here. Thank so much if anyone could help me.

var ('Omegabar Qbar ')
Ybar=1
pbeta=.99
epsilon=-1 
sigma=.03 
theta=.3 
k=13
delta=.025
rho=1.1
Ybar=1;
Lambdabar=pbeta
Rlbar=1+1/Qbar
Rbar=1/pbeta
Rebar=1/pbeta
munbar=LambdabarOmegabarRbar
musbar=LambdabarOmegabar(Rlbar-Rbar)
muebar=LambdabarOmegabar(1-pbeta)
munbar=Lambdabar*Omegabar*Rbar
musbar=Lambdabar*Omegabar*(Rlbar-Rbar)
muebar=Lambdabar*Omegabar*(1-pbeta)
epersbar=musbar/muebar
xbar=-epersbar+(epersbar^(2)+(2/k)(1-epsilonepersbar))^(1/2)
Thetabar=theta(1+epsilonxbar+(k/2)xbar^(2))
phibar=Omegabar/(Thetabar-musbar-xbarmuebar)
Ebar=xbarQbar
xbar=-epersbar+(epersbar^(2)+(2/k)*(1-epsilon*epersbar))^(1/2)
Thetabar=theta*(1+epsilon*xbar+(k/2)*xbar^(2))
phibar=Omegabar/(Thetabar-musbar-xbar*muebar)
Ebar=xbar*Qbar
Nbar=Qbar/phibar
Dbar=Qbar-Nbar-Ebar
Nnbar=delta(Qbar+Ybar)
Nnbar=delta*(Qbar+Ybar)
  
eq1=Omegabar==(sigma+(1-sigma)(munbar+phibarmusbar+phibarxbarmuebar))
eq2=Nbar==(1-sigma)(Qbar+Ybar-RebarEbar-Rbar*Dbar)+Nnbar
 
eq1=Omegabar==(sigma+(1-sigma)*(munbar+phibar*musbar+phibar*xbar*muebar))
eq2=Nbar==(1-sigma)*(Qbar+Ybar-Rebar*Ebar-Rbar*Dbar)+Nnbar
solve([eq1,eq2], Qbar)
 
Qbar)
-----------
(Omegabar, Qbar)
[][]