diff --git a/matlab/compute_trend_coefficients.m b/matlab/compute_trend_coefficients.m
index 48f08bbaf377096a6aa7ce0747ca82fc25e8aa00..2dbed2ae2e1e183f2f39878966ef16b25d993485 100644
--- a/matlab/compute_trend_coefficients.m
+++ b/matlab/compute_trend_coefficients.m
@@ -1,5 +1,5 @@
 function [trend_addition, trend_coeff]=compute_trend_coefficients(M_,DynareOptions,nvarobs,ntobs)
-% [trend_addition, trend_coeff]=compute_trend_coefficients(DynareOptions,nvarobs,ntobs)
+% [trend_addition, trend_coeff]=compute_trend_coefficients(M_,DynareOptions,nvarobs,ntobs)
 % Computes the trend coefficiencts and the trend, accounting for
 % prefiltering
 %
diff --git a/tests/kalman_filter_smoother/fs2000_smoother_only_ns.mod b/tests/kalman_filter_smoother/fs2000_smoother_only_ns.mod
index 59e1980f5f4f63ef43fee01568cb46ad62f6bb89..5c49531bd80f30ea1262b788242135edf83f1711 100644
--- a/tests/kalman_filter_smoother/fs2000_smoother_only_ns.mod
+++ b/tests/kalman_filter_smoother/fs2000_smoother_only_ns.mod
@@ -107,9 +107,24 @@ P_obs (log(mst)-gam);
 Y_obs (gam);
 end;
 
+fsdat_simul;
+P_obs_data=log(P_obs(1:192));
+Y_obs_data=log(Y_obs(1:192));
+
 estimation(order=1, datafile=fsdat_simul, mode_compute=0,nobs=192, loglinear,diffuse_filter, smoother,smoothed_state_uncertainty) m P c e W R k d n l gy_obs gp_obs y dA;
+if max(abs(oo_.SmoothedVariables.Y_obs-Y_obs_data))>1e-10 || max(abs(oo_.SmoothedVariables.P_obs-P_obs_data))>1e-10
+    error('Smoothed observables do not match')
+end
 estimation(order=1, datafile=fsdat_simul, mode_compute=0,nobs=192, loglinear,diffuse_filter, smoother,kalman_algo=3,smoothed_state_uncertainty) m P c e W R k d n l gy_obs gp_obs y dA;
+if max(abs(oo_.SmoothedVariables.Y_obs-Y_obs_data))>1e-10 || max(abs(oo_.SmoothedVariables.P_obs-P_obs_data))>1e-10
+    error('Smoothed observables do not match')
+end
 estimation(order=1, datafile=fsdat_simul, mode_compute=0,nobs=192, loglinear,diffuse_filter, smoother,kalman_algo=4,smoothed_state_uncertainty) m P c e W R k d n l gy_obs gp_obs y dA;
+if max(abs(oo_.SmoothedVariables.Y_obs-Y_obs_data))>1e-10 || max(abs(oo_.SmoothedVariables.P_obs-P_obs_data))>1e-10
+    error('Smoothed observables do not match')
+end
+
+
 
 /*
  * The following lines were used to generate the data file. If you want to