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 residualRp = delz*k*DDr*T - b*emiss*(T^4-300^4) + Q;% Boundary condition residuals% symmetry at r=0Rp = setbval(Rp,Dr*T,'r','min');% radiation from wafer edgeRp = setbval(Rp,k*Dr*T-b*emiss*(300^4-T^4),'r','max');A = set(A,'resid',Rp);