Hi there,

I have a model with an occasionally binding constraint and I am using the extended_path algorithm (1600 periods, just one shock with sd 0.002) Unfortunately I get the error message reported below. Any clues on how I can proceed? Any hint and/or suggestion will be very appreciated.

Barbara

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In extended_path at 213

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 63

In extended_path at 344

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 63

In extended_path at 344

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 63

In extended_path at 344

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 63

In extended_path at 344

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 63

In extended_path at 344

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 63

In extended_path at 344

In binding at 345

In dynare at 120

Warning: Matrix is singular to working precision.

In solve_perfect_foresight_model at 77

In homotopic_steps at 136

In extended_path at 344

In binding at 345

In dynare at 120

??? Error using ==> homotopic_steps at 151

extended_path::homotopy: Oups! I did my best, but I

am not able to simulate this model…

Error in ==> extended_path at 344

[INFO,tmp] =

homotopic_steps(endo_simul,exo_simul_1,.5,.01,pfm1);

Error in ==> binding at 345

oo_.endo_simul = oo_.steady_state,

extended_path(], 1600) ];

Error in ==> dynare at 120

evalin(‘base’,fname) ;