clear all; 
close all;

%mu_params = 0.08:0.001:0.12;



%for i=1:length(mu_params)


n_initial = 0.5 ;  			 		%	//labor growth rate
beta_initial = 0.9915;					%//discount factor
alpha = 0.35;					%//input share of capital in Cobb-Douglas function
theta = 0.13;			%//tax rate
A_initial = 1.168; 				%//initial TFP
mu_initial = 0.090;				%//initial proportion of GDP spent on healthcare
%mu_initial = mu_params(i);
sigma_initial = 1;			%//initial value for sigma: когда все шоки кроме амортизации, то лучше ставить сигму=1
z = 0.98;						%// Autoregressive parameter of TFP
q = 0.95;						%// Autoregressive parameter of proportion of GDP spent on healthcare
sig_mu = 0.0001;					%//sigma shock size of proportion of GDP spent on healthcare
sig_A = 0.0001;					%//sigma shock size of TFP
sig_n = 0.0001;					%//sigma for pop size
sig_sigma = 0.0001;				%//sigma for sigma depreciation
sig_beta = 0.0001;				%//sigma for disc factor
sig_rho = 0.0001;
sig_theta = 0.001;				%//sigma shock size of tax rate
tau = 0.95;						%// Autoregressive parameter of tax rate

save param_nc n_initial beta_initial alpha theta A_initial mu_initial sigma_initial z q sig_A sig_mu sig_n sig_sigma sig_beta sig_rho sig_theta tau 

%try
	dynare OLG_dip2_amend noclearall nolog		
%	u_lifetime_eps_mu_mat(i,:) = oo_.mean(20);
%catch
%	u_lifetime_eps_mu_mat(i,:) = NaN;
%end

%end

%plot(mu_params, u_lifetime_eps_mu_mat); xlabel('\mu'); ylabel('lifetime u (mean)');