Added solve_linear_system() funtion to nodal*.py and changed
function names / conventions in _NLFunction and fsolve.