function g = rhs(A)% rhs.m Nonlinear equation model object residual update.%% INPUT/OUTPUT PARAMETERS%        A : an object derived from the naemodel classunpack(A)g{1} = e_wq*F*(Tq^4-Tw^4) + e_ws*(2-F)*(1-Tw^4) ...        + k_w/alf*(Tq-Tw) + q_w*u;g{2} = e_qw*F*(Tw^4-Tq^4) + e_qs*(1-F)*(1-Tq^4) ...        + nu_q*(1-Tq) + q_q*u;