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-2006unpack(A,'var')unpack(A,'param')% Interior residualRT = DDz*T + u;% Boundary condition residualsRT = setbval(RT,T,'z','min');RT = setbval(RT,T,'z','max');A = set(A,'resid',{RT});