<<" it_ : Int, time period for exogenous variables for which to evaluate the model"<<endl

<<" residual : Vector{Float64}(model_.eq_nbr), residuals of the dynamic model equations in order of declaration of the equations."<<endl

<<" g1 : Matrix{Float64}(model_.eq_nbr, num_dynamic_vars), Jacobian matrix of the dynamic model equations"<<endl

<<" The rows and columns respectively correspond to equations in order of declaration and variables in order"<<endl

<<" stored in model_.lead_lag_incidence"<<endl

<<" g2 : spzeros(model_.eq_nbr, (num_dynamic_vars)^2) Hessian matrix of the dynamic model equations"<<endl

<<" The rows and columns respectively correspond to equations in order of declaration and variables in order"<<endl

<<" stored in model_.lead_lag_incidence"<<endl

<<" g3 : spzeros(model_.eq_nbr, (num_dynamic_vars)^3) Third order derivative matrix of the dynamic model equations;"<<endl

<<" The rows and columns respectively correspond to equations in order of declaration and variables in order"<<endl

<<" stored in model_.lead_lag_incidence"<<endl<<endl

<<"## Remarks ##"<<endl

<<" [1] `num_dynamic_vars` is the number of non zero entries in the lead lag incidence matrix, `model_.lead_lag_incidence.`"<<endl

<<" [2] The size of `T`, ie the value of `num_temp_terms`, depends on the version of the dynamic model called. The number of temporary variables"<<endl

<<" used for the different returned objects (residuals, jacobian, hessian or third order derivatives) is given by the elements in `tmp_nbr`"<<endl

<<" exported vector. The first element is the number of temporaries used for the computation of the residuals, the second element is the"<<endl

<<" number of temporaries used for the evaluation of the jacobian matrix, etc. If one calls the version of the dynamic model computing the"<<endl

<<" residuals, the jacobian and hessian matrices, then `T` must have at least `sum(tmp_nbr[1:3])` elements."<<endl