function  param_save(beta, phiu,fracb,fracf)

beta=0.995;
sp=0.03; % steady state net spread
Ls=1-0.058; 
bh=0.7; % habit formation
eta=1; % inverse frisch elasticity 
Bw=0.1; % unemployment subsidy over steady state wage
kappa=1-40^(-1); % decay parameter for both private and goverment
sigma=0.95; % survival probability
Delta=1/3; % steady state difference in recoverability government vs private
deltaN=0.005; % seperation rate
mubar=1; % maching and searching
epsilon=11; %elasticity good
phip=0.75; % price stickiness
gammap=0; % price indexation
alpha = 0.33; %production function
delta=0.025; % depreciation rate
xi=1; % vancancy cost
kappaI=2; % capital adjustment cost
bgy=0.2; % the share of government bond over gdp
gy=0.2; % government spending share
Pis=1.02; % steady state inflation
Qgammas=0.9; % vacancy filled probability
As=1;



% Monetary policy rate
phipi=1.5; 
rhor=0.8;
sr=0.0025;

% shock parameters
rhoA=0.95;
rhoG=0.095;
rhot=0.098;
rhof=0.8;
rhob=0.8;
sA = 0.0065;
sG = 0.01;
st = 0.04;
sf = sr; 
sb = sr;




pk = 1;

if nargin < 4
    phiu=0.25;
	if nargin < 3
        fracb = 0.06*4;
        fracf = 0;
    end
end


save param beta sp Ls bh eta Bw kappa sigma Delta deltaN mubar epsilon...
           phip gammap alpha delta xi kappaI bgy gy Pis Qgammas phipi rhor...
           sr rhoA rhoG rhot rhof rhob sA sG st sf sb pk phiu fracb fracf As
      