function g = rhs(A)% residual.m Nonlinear equation model object residual update.%% INPUT/OUTPUT PARAMETERS%        A : an object derived from the neqmodel 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;