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)

conc = e*Bconc + concBC;

% Interior residuals

resid = wip( Dab*(DDx*conc + DDy*conc) - vx*(Dx*conc) - vy*(Dy*conc), Bconc );

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