residual.m

MDPSAS工具箱是马里兰大学开发的

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-2003

unpack(A,'var')
unpack(A,'param')

% Interior residual

Rp = delz*k*DDr*T - b*emiss*(T^4-300^4) + Q;

% Boundary condition residuals

% symmetry at r=0

Rp = setbval(Rp,Dr*T,'r','min');

% radiation from wafer edge

Rp = setbval(Rp,k*Dr*T-b*emiss*(300^4-T^4),'r','max');

A = set(A,'resid',Rp);