diff --git a/src/auxiliary_initialization.m b/src/auxiliary_initialization.m
index 98f79aecb36478b5b7a7c72eedb091a037495280..002a8038af4ba412d0bcc956964570cee6aaf4ac 100644
--- a/src/auxiliary_initialization.m
+++ b/src/auxiliary_initialization.m
@@ -2,7 +2,7 @@ function initial_distribution = auxiliary_initialization(ReducedForm,Y,start,Par
 
 % Evaluates the likelihood of a nonlinear model with a particle filter allowing eventually resampling.
 
-% Copyright (C) 2011-2014 Dynare Team
+% Copyright (C) 2011-2017 Dynare Team
 %
 % This file is part of Dynare (particles module).
 %
diff --git a/src/auxiliary_particle_filter.m b/src/auxiliary_particle_filter.m
index 42d57ca450fda4de01abdb41d85f4c485bb6840b..61a43bc13353c6f2f3acf33ce6be8e13fb7cffeb 100644
--- a/src/auxiliary_particle_filter.m
+++ b/src/auxiliary_particle_filter.m
@@ -3,7 +3,7 @@ function [LIK,lik] = auxiliary_particle_filter(ReducedForm,Y,start,ParticleOptio
 % Evaluates the likelihood of a nonlinear model with the auxiliary particle filter
 % allowing eventually resampling.
 %
-% Copyright (C) 2011-2015 Dynare Team
+% Copyright (C) 2011-2017 Dynare Team
 %
 % This file is part of Dynare (particles module).
 %
diff --git a/src/conditional_filter_proposal.m b/src/conditional_filter_proposal.m
index 0c4121884b74a3971a19a4464b3a2a4c1dc6cd42..5137b83f7036fac713b43087d16b3928ec64c50b 100644
--- a/src/conditional_filter_proposal.m
+++ b/src/conditional_filter_proposal.m
@@ -18,7 +18,7 @@ function [ProposalStateVector,Weights] = conditional_filter_proposal(ReducedForm
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2012-2013 Dynare Team
+% Copyright (C) 2012-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/conditional_particle_filter.m b/src/conditional_particle_filter.m
index 01437773e74f9b928c2545a1597276d26028df47..6216cf15c928c07163876e038ff421db376cc481 100644
--- a/src/conditional_particle_filter.m
+++ b/src/conditional_particle_filter.m
@@ -37,7 +37,7 @@ function [LIK,lik] = conditional_particle_filter(ReducedForm,Y,start,ParticleOpt
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2010 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/fit_gaussian_mixture.m b/src/fit_gaussian_mixture.m
index cc34e66b6e10d21a764215488506e0a7648e2180..672c90e68104cfcee4431d7ab0a16f2ff0e14e69 100644
--- a/src/fit_gaussian_mixture.m
+++ b/src/fit_gaussian_mixture.m
@@ -1,6 +1,6 @@
 function [StateMu,StateSqrtP,StateWeights] = fit_gaussian_mixture(X,X_weights,StateMu,StateSqrtP,StateWeights,crit,niters,check)
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/gaussian_densities.m b/src/gaussian_densities.m
index 2246e45e15d504a7616457e105ccc0cbdb495819..889a5bbe239d6d22ff60d441ce8f0392a135f245 100644
--- a/src/gaussian_densities.m
+++ b/src/gaussian_densities.m
@@ -19,7 +19,7 @@ function IncrementalWeights = gaussian_densities(obs,mut_t,sqr_Pss_t_t,st_t_1,sq
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2010 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/gaussian_filter.m b/src/gaussian_filter.m
index d7e82df1bd218fbd1bd32c209d7016633fc27b6c..e885cce725e6463445ccb91476d21a1e31058b57 100644
--- a/src/gaussian_filter.m
+++ b/src/gaussian_filter.m
@@ -30,7 +30,7 @@ function [LIK,lik] = gaussian_filter(ReducedForm, Y, start, ParticleOptions, Thr
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2015 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/gaussian_filter_bank.m b/src/gaussian_filter_bank.m
index 3e9deb31cdae9665d09d73afe63d2dfc898445db..fb229963171752159b3d3aa5d5eaff50660fa9a4 100644
--- a/src/gaussian_filter_bank.m
+++ b/src/gaussian_filter_bank.m
@@ -19,7 +19,7 @@ function [PredictedStateMean,PredictedStateVarianceSquareRoot,StateVectorMean,St
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2012 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/gaussian_mixture_densities.m b/src/gaussian_mixture_densities.m
index 7efe1ce015b8c60a9b1874b571a4e32ccbc0e42f..c826059345220f89b85bdd0908c17d9daba76c4d 100644
--- a/src/gaussian_mixture_densities.m
+++ b/src/gaussian_mixture_densities.m
@@ -21,7 +21,7 @@ function  IncrementalWeights = gaussian_mixture_densities(obs,StateMuPrior,State
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2012 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/gaussian_mixture_filter.m b/src/gaussian_mixture_filter.m
index eaab0070f67feaf05d667e4869175e35ca25353a..f744210d8111b654efa929f180ce51dd47b4c0ee 100644
--- a/src/gaussian_mixture_filter.m
+++ b/src/gaussian_mixture_filter.m
@@ -35,7 +35,7 @@ function [LIK,lik] = gaussian_mixture_filter(ReducedForm,Y,start,ParticleOptions
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2013 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/gaussian_mixture_filter_bank.m b/src/gaussian_mixture_filter_bank.m
index 0171927126c0f629c10bcd7a6988597ed313dc8e..9c47b102145ee3b55e1aea5d3cd37d6532d88fdf 100644
--- a/src/gaussian_mixture_filter_bank.m
+++ b/src/gaussian_mixture_filter_bank.m
@@ -23,7 +23,7 @@ function [StateMuPrior,StateSqrtPPrior,StateWeightsPrior,StateMuPost,StateSqrtPP
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2012 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/importance_sampling.m b/src/importance_sampling.m
index b4a6daf4a55266a2322ccb3d30a6ecd80c5a198f..ae20b02a79ce01f6c1d6266f8e8fb1b362b94d8b 100644
--- a/src/importance_sampling.m
+++ b/src/importance_sampling.m
@@ -1,6 +1,6 @@
 function State_Particles = importance_sampling(StateMuPost,StateSqrtPPost,StateWeightsPost,numP)
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/measurement_equations.m b/src/measurement_equations.m
index 7cf4d9320b76f4c13db533b5fbf1818ae0822e7a..1c386b6fd7b1b541ffeec29fc2c8406b6372737a 100644
--- a/src/measurement_equations.m
+++ b/src/measurement_equations.m
@@ -1,6 +1,6 @@
 function measure = measurement_equations(StateVectors,ReducedForm,ThreadsOptions)
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/multivariate_smooth_resampling.m b/src/multivariate_smooth_resampling.m
index e62111cc1642cd86ccdf56de864d74822fead1e8..01c9208a7eb8d48c1013d4e9b2bc3ec4471cf908 100644
--- a/src/multivariate_smooth_resampling.m
+++ b/src/multivariate_smooth_resampling.m
@@ -38,7 +38,7 @@ function new_particles = multivariate_smooth_resampling(particles,weights)
 %! @end deftypefn
 %@eod:
 
-% Copyright (C) 2012-2013 Dynare Team
+% Copyright (C) 2012-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/mykmeans.m b/src/mykmeans.m
index 4506aca25f85c9d1ec032dd4e8f25a4d5ba1c9b7..c39afedb9d9adade82fe846299babcf5db6f2b16 100644
--- a/src/mykmeans.m
+++ b/src/mykmeans.m
@@ -1,6 +1,6 @@
 function [c,SqrtVariance,Weights] = mykmeans(x,g,init,cod)
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/neff.m b/src/neff.m
index 9e02b679e7519765f8a108d63548514e15794928..e0fb4f70ea431d6b604898f142200e4950acaf75 100644
--- a/src/neff.m
+++ b/src/neff.m
@@ -1,7 +1,7 @@
 function n = neff(w)
 % Evaluates the criterion for resampling
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2014 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/nonlinear_kalman_filter.m b/src/nonlinear_kalman_filter.m
index 608fbacec644b37fa2fb8c40adbadb9167e0ace2..6710e69d477aca23dfdd945a03ca3db53e928139 100644
--- a/src/nonlinear_kalman_filter.m
+++ b/src/nonlinear_kalman_filter.m
@@ -30,7 +30,7 @@ function [LIK,lik] = nonlinear_kalman_filter(ReducedForm, Y, start, ParticleOpti
 %
 % NOTES
 %   The vector "lik" is used to evaluate the jacobian of the likelihood.
-% Copyright (C) 2009-2015 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/probability.m b/src/probability.m
index 94f09b9cb5b6aa86e170aaaeef9cd4331ff0fcab..68efb3fa1e874611d9f3e84951409df350c96049 100644
--- a/src/probability.m
+++ b/src/probability.m
@@ -1,6 +1,6 @@
 function [prior,likelihood,C,posterior] = probability(mu,sqrtP,prior,X)
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/probability2.m b/src/probability2.m
index dfc71ff24b7c16e3a7304b7a339865d2d9cf58a2..268cd03941ec9fc8934988f114caad223bb313d1 100644
--- a/src/probability2.m
+++ b/src/probability2.m
@@ -14,7 +14,7 @@ function [density] = probability2(mu,S,X)
 %
 % NOTES
 %
-% Copyright (C) 2009-2012 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/probability3.m b/src/probability3.m
index 5e2a9b4a8af1edf41c1f4db9f056a479b115d232..710a84267effdcdd1383fe60e26a139113eedd59 100644
--- a/src/probability3.m
+++ b/src/probability3.m
@@ -1,6 +1,6 @@
 function [prior,likelihood,C,posterior] = probability3(mu,sqrtP,prior,X,X_weights)
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/resample.m b/src/resample.m
index 3fee0c355f34cf73ef5e5460ecbb6e28e504cad5..e29d62e4f2d3ad03e51933bd36419602f6992022 100644
--- a/src/resample.m
+++ b/src/resample.m
@@ -36,7 +36,7 @@ function resampled_output = resample(particles,weights,ParticleOptions)
 %! @end deftypefn
 %@eod:
 
-% Copyright (C) 2011-2013 Dynare Team
+% Copyright (C) 2011-2014 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/residual_resampling.m b/src/residual_resampling.m
index 9a7d8cce84fea3fb1c678af737a730223d1184c6..d4d250309ce354b4e6781a829ab057c62f93e03b 100644
--- a/src/residual_resampling.m
+++ b/src/residual_resampling.m
@@ -30,7 +30,7 @@ function return_resample = residual_resampling(particles,weights,noise)
 %! @end deftypefn
 %@eod:
 
-% Copyright (C) 2011-2013 Dynare Team
+% Copyright (C) 2011-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/sequential_importance_particle_filter.m b/src/sequential_importance_particle_filter.m
index 51e0bd7a1f192c6e6be07e77833bcb36e2cdd72f..72a789ccd62e0b7d2e4d50076ce66cabc125a9d8 100644
--- a/src/sequential_importance_particle_filter.m
+++ b/src/sequential_importance_particle_filter.m
@@ -2,7 +2,7 @@ function [LIK,lik] = sequential_importance_particle_filter(ReducedForm,Y,start,P
 
 % Evaluates the likelihood of a nonlinear model with a particle filter (optionally with resampling).
 
-% Copyright (C) 2011-2014 Dynare Team
+% Copyright (C) 2011-2015 Dynare Team
 %
 % This file is part of Dynare (particles module).
 %
diff --git a/src/solve_model_for_online_filter.m b/src/solve_model_for_online_filter.m
index 1ea2e67a474791900ade4a34db94771425d1feca..74d2fd925544d28c299db7ef792d3e31edacd0a0 100644
--- a/src/solve_model_for_online_filter.m
+++ b/src/solve_model_for_online_filter.m
@@ -101,7 +101,7 @@ function [ys,trend_coeff,exit_flag,info,Model,DynareOptions,BayesInfo,DynareResu
 %! @end deftypefn
 %@eod:
 
-% Copyright (C) 2013 Dynare Team
+% Copyright (C) 2013-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/spherical_radial_sigma_points.m b/src/spherical_radial_sigma_points.m
index db71f61196956cba96aafca0e6b020223a5c9f9d..aeaf7d5c5869c881e80985e29376c8dcbeb61eeb 100644
--- a/src/spherical_radial_sigma_points.m
+++ b/src/spherical_radial_sigma_points.m
@@ -15,7 +15,7 @@ function [nodes,weights] = spherical_radial_sigma_points(n)
 %
 % NOTES
 %
-% Copyright (C) 2009-2012 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/traditional_resampling.m b/src/traditional_resampling.m
index 87f5503681577f95a0cc8ae32c57fd4464fe0e3c..ce6ed90c45b6cb7db5b272fd244868fd8b40019f 100644
--- a/src/traditional_resampling.m
+++ b/src/traditional_resampling.m
@@ -33,7 +33,7 @@ function return_resample = traditional_resampling(particles,weights,noise)
 %! @end deftypefn
 %@eod:
 
-% Copyright (C) 2011-2013 Dynare Team
+% Copyright (C) 2011-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/univariate_smooth_resampling.m b/src/univariate_smooth_resampling.m
index e89b2cc6dad646fe9123b9cad5f27525ef97554a..1a8d10439d09fa48af85ed373b686d1b2e750aa6 100644
--- a/src/univariate_smooth_resampling.m
+++ b/src/univariate_smooth_resampling.m
@@ -34,7 +34,7 @@ function new_particles = univariate_smooth_resampling(weights,particles,number_o
 %! @end deftypefn
 %@eod:
 
-% Copyright (C) 2012 Dynare Team
+% Copyright (C) 2012-2017 Dynare Team
 %
 % This file is part of Dynare.
 %
diff --git a/src/unscented_sigma_points.m b/src/unscented_sigma_points.m
index d1cba6eb59e10f3c5984a3945d1d61a9f02693bc..54c5b0b0de690bfd7ae2eea4ab306870e0d55493 100644
--- a/src/unscented_sigma_points.m
+++ b/src/unscented_sigma_points.m
@@ -14,7 +14,7 @@ function [nodes,W_m,W_c] = unscented_sigma_points(n,ParticleOptions)
 %
 % NOTES
 %
-% Copyright (C) 2009-2012 Dynare Team
+% Copyright (C) 2009-2017 Dynare Team
 %
 % This file is part of Dynare.
 %