diff --git a/DynamicModel.cc b/DynamicModel.cc
index 336f2cd4..694dab45 100644
--- a/DynamicModel.cc
+++ b/DynamicModel.cc
@@ -2390,6 +2390,8 @@ DynamicModel::writeDynamicModel(ostream &DynamicOutput, bool use_dll, bool julia
                     << "                  steady_state::Vector{Float64}, it_::Int, "
                     << "residual::Vector{Float64})" << endl
                     << "#=" << endl << comments.str() << "=#" << endl
+                    << "@assert size(y) == " << dynJacobianColsNbr << endl
+                    << "@assert size(params) == " << symbol_table.param_nbr() << endl
                     << "@assert size(residual) == " << nrows << endl
                     << "fill!(residual, 0.0)" << endl
                     << "#" << endl
diff --git a/StaticModel.cc b/StaticModel.cc
index 1f5ef56e..5407ffcf 100644
--- a/StaticModel.cc
+++ b/StaticModel.cc
@@ -1453,6 +1453,9 @@ StaticModel::writeStaticModel(ostream &StaticOutput, bool use_dll, bool julia) c
                    << "params::Vector{Float64}," << endl
                    << "                 residual::Vector{Float64})" << endl
                    << "#=" << endl << comments.str() << "=#" << endl
+                   << "@assert size(y) == " << symbol_table.endo_nbr() << endl
+                   << "@assert size(x) == " << symbol_table.exo_nbr() << endl
+                   << "@assert size(params) == " << symbol_table.param_nbr() << endl
                    << "@assert size(residual) == " << equations.size() << endl
                    << "fill!(residual, 0.0)" << endl
                    << "#" << endl