recsSolveDeterministicPbSP Solves a perfect foresight problem
 
  X = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS) tries to find the perfect
  foresight solution of the model defined in the object MODEL. The initial
  values of state variable are provided by the 1-by-d vector S0. Time horizon
  (number of periods) is given by the integer T. XSS, ZSS and SSS are,
  respectively, a 1-by-m vector containing the values of response variables at the
  deterministic steady state, a 1-by-p vector containing the values of expectations
  variables at steady state, and a 1-by-d vector containing the values of the state
  variables at steady state. recsSolveDeterministicPbSP returns X, a T-by-m matrix,
  the value of response variables over the time horizon.
 
  X = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS,OPTIONS) solves the problem with the
  parameters defined by the structure OPTIONS. The fields of the structure are
     eqsolver         : 'fsolve', 'lmmcp' (default), 'ncpsolve' or 'path'
     eqsolveroptions  : options structure to be passed to eqsolver (default:
                        empty structure)
 
  [X,S] = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS,...) returns S, a T-by-d
  matrix, containing the value of state variables over the time horizon.
 
  [X,S,Z] = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS,...) returns Z, a
  T-by-p matrix, containing the value of expectations variables over the time
  horizon.
 
  [X,S,Z,F] = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS,...) returns F, a
  T-by-m matrix, containing the values of equilibrium equations over the time
  horizon.
 
  [X,S,Z,F,EXITFLAG] = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS,...)
 
  [X,S,Z,F,EXITFLAG,N] = recsSolveDeterministicPbSP(MODEL,S0,T,XSS,ZSS,SSS,...)