New version of recsAuxiliary:
 - Faster because it uses direct multiplication with the basis matrix and a Newton-Krylov solver.
 - Require independent functions for the definition of auxiliary variables, so no need to extend the usual model file.
Adjust recsinterinit.m so that it outputs s as a field of the structure
Compatibility of recsSolveREE.m with the new behavior of recsinterpinit.m
Bug correction in runeqsolver.
Fix compatibility with MATLAB R2013a in catstruct.m