diff --git a/matlab/dsge_likelihood.m b/matlab/dsge_likelihood.m
index a646f82af4a54546efdc534282f7314f7e633e1c..3385bb253d03d6e781fed17587fc82b6c7382353 100644
--- a/matlab/dsge_likelihood.m
+++ b/matlab/dsge_likelihood.m
@@ -145,7 +145,7 @@ function [fval,DLIK,Hess,exit_flag,ys,trend_coeff,info,Model,DynareOptions,Bayes
 % In summary, an initial call to the present function, without triggering
 % any condition, guarantees that 'penalty' is properly initialized when needed.
 
-persistent penalty prior_penalty
+persistent penalty
 
 % Initialization of the returned variables and others...
 fval        = [];
@@ -189,10 +189,6 @@ end
 if ~isequal(DynareOptions.mode_compute,1) && any(xparam1<BayesInfo.lb)
     k = find(xparam1<BayesInfo.lb);
     fval = penalty+sum((BayesInfo.lb(k)-xparam1(k)).^2);
-    xtemp=xparam1;
-    xtemp(k)=BayesInfo.lb(k);
-    prior_correction = max(0,-prior_penalty-priordens(xtemp,BayesInfo.pshape,BayesInfo.p6,BayesInfo.p7,BayesInfo.p3,BayesInfo.p4));
-    fval = fval+prior_correction;
     exit_flag = 0;
     info = 41;
     if analytic_derivation,
@@ -205,10 +201,6 @@ end
 if ~isequal(DynareOptions.mode_compute,1) && any(xparam1>BayesInfo.ub)
     k = find(xparam1>BayesInfo.ub);
     fval = penalty+sum((xparam1(k)-BayesInfo.ub(k)).^2);
-    xtemp=xparam1;
-    xtemp(k)=BayesInfo.ub(k);
-    prior_correction = max(0,-prior_penalty -priordens(xtemp,BayesInfo.pshape,BayesInfo.p6,BayesInfo.p7,BayesInfo.p3,BayesInfo.p4));
-    fval = fval+prior_correction;
     exit_flag = 0;
     info = 42;
     if analytic_derivation,
@@ -777,7 +769,6 @@ DynareOptions.kalman_algo = kalman_algo;
 
 % Update the penalty.
 penalty = fval;
-prior_penalty = -lnprior;
 
 if analytic_derivation==0 && nargout==2,
     lik=lik(start:end,:);