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 = DDx*T + DDy*T;

% Boundary condition residuals

% radiation from bottom

Rp  = setbval(Rp,Dx*T-(T^4-0),'x','min');

% radiation from top

Rp  = setbval(Rp,Dx*T-(0-T^4),'x','max');

% radiation from left side

Rp  = setbval(Rp,Dy*T-(T^4-0),'y','min');

% constant temp BC

Rp = setbval(Rp,T-Tbc,'y','max');

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