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-2003unpack(A,'var')unpack(A,'param')% Interior residualRT = kwaf*DDr*T + kwaf*DDz*T;% Boundary condition residuals% symmetry at centerRT  = setbval(RT,Dr*T,'r','min');% radiation from edgeRT  = setbval(RT, kwaf*Dr*T ...    - Esat*Ewaf*alph*sig*(Tsus^4-T^4),'r','max');% radiation from bottomRT  = setbval(RT, -kwaf*Dz*T ...    - Esat*Ewaf*alph*sig*(Tsus^4-T^4),'z','min');% radiation/conduction to reactor topRT = setbval(RT, -kwaf*Dz*T ...    - kH2*(T-Troof)/Lrxr ...    - Eroof*Ewaf*alph*sig*(T^4-Troof^4),'z','max');% Wafer bow modeldelT = getbval(T,'z','min') - getbval(T,'z','max');Ry = DDr*y - betaSapp*delT/Lw;Ry = setbval(Ry,Dr*y,'r','min');Ry = setbval(Ry,y,'r','max');A = set(A,'resid',{RT, Ry});