function  param_save(psif,psib,fracb,fracf,rhof_QT,rhob_QT)

alpha = 0.33; % production function
beta = 0.995; % discount factor
delta0 = 0.025; % steady state depreciation rate
Deltas = 1/3; % steady state difference in recoverability government vs private
sp = 0.03; % steady state net spread
Pis = 1; % steady state gross inflation
kappaf = 1-40^(-1); % decay parameter private
kappab = 1-40^(-1); % decay parameter government
epsip = 11; % elasticity goods
epsiw = 11; % elasticity labor
phiw = 0.75; % wage stickiness
gammap = 0; % price indexation
gammaw = 0; % wage indexation
gy = 0.2; % government spending share 
bh = 0.7; % habit formation 
eta = 1; % inverse Frisch elasticity
ff = 0.1; % fraction of private debt held by central bank
levs = 4; % steady state leverage 
sigma = 0.95; % survival probability
psi = .81;
bgy = 0.41*4; % steady state debt-gdp ratio
rotemberg=118;

% Monetary policy rule
rhor = 0.8;
phipi = 1.5;
phiy = 0.25;

% investment and utilization adjustment costs
kappai = 2;
delta2 = 0.01;

% shock parameters
rhoA = 0.95;
rhom = 0.9;
rhoc = 0.9;
rhol = 0.9;
rhoG = 0.95;
rhoB = 0.98;
rhof = rhor;
rhob = rhor;
rhot = 0.98;
rhoD = 0.98;
sA = 0.0065;
sm = 0.00;
sc = 0.01;
sl = 0.00;
sG = 0.01;
st = 0.04;
sD = 0.00; 
sB = 0.00;

Ld = 1;
u = 1;
pk = 1;

if nargin < 5
    rhof_QT = rhor;
    rhob_QT = rhor;
	if nargin < 3
        fracb = 0.06*4;
        fracf = 0;
    end
end

sr =  0.0025; 
sf = sr; 
sb = sr;

save param alpha beta delta0 Deltas sp Pis kappaf kappab epsip epsiw ...
             phiw gammap gammaw gy bh eta ff levs sigma ...
            psi bgy rotemberg rhor phipi phiy kappai delta2 rhoA rhom rhoc rhol ...
            rhoG rhoB rhof rhob rhot rhoD sr sA sm sc sl sG sf st sD sb ...
            sB psif psib Ld u pk rhof_QT rhob_QT fracb  fracf 