I'm analysing a thermodynamic model where I must solve a system of equations for a certain component.
The enthalpy however depends on the temperature and the concentration, which are also unkown variables. is it possible to add the function of the relation between the enthalpy and the concentration in the system of equations?
See below the script that I tried. (eqn 7 and 8 are the ones that don't work currently)
%%variables that are known.
%%determine unknown variables.
syms x_5 y_6 m_5 m_6 h_5 h_6 Q_rect T
[solm5 ,solm6, solh5, solh6 ,Qrect ,T]=vpasolve([eqn1,eqn2,eqn3,eqn4,eqn5,eqn6,eqn7,eqn8],[x_5, y_6,m_5,m_6,h_5,h_6,Q_rect,T,h_5ex,h_6ex]);