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

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

if colloc

   c = a;
   Rp = DDr*c - tm^2*c^2;
   Rp = setbval(Rp,Dr*c,'r','min'); % BC residuals
   Rp = setbval(Rp,c-1, 'r','max');

else

   c  = One + a*P;
   R  = DDr*c - tm^2*c^2;
   Rp = wip(R,P);

   end

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