Define the function "calc_nodal_eos_sys" to solve the EOS on each node