Add missing parenthesis
Header fix
Fix bug in linsolve when second argument is a zero-column matrix

Fixes crash of 1st-order purely forward models under Octave
Deactivate 2nd and 3rd order approximation of purely forward models

There are some strange crashes in both normal mode (gensylv) and k_order_solver
mode