function A = residual(A)% residual.m Nonlinear equation model object residual update.%% INPUT/OUTPUT PARAMETERS%        A : an object derived from the neqmodel class%% Copyright (c) by Raymond A. Adomaitis, 1998-2005unpack(A,'var')unpack(A,'param')Rp = Esus*Eroof*sig*(Tsus^4-T^4) + kH2*(Tsus-T)/Lrxr - kRoof*(T-300)/Lins;A = set(A,'resid',Rp);