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')

r = get(S,'qp');
r = r{1};

theta = [-pi:0.01:pi];
for i = 1:N
   qlamp(i,1) = (h_lamp*Q_lamp/8/pi^2) ...
   *sum( (h_lamp^2+r(i)^2+r_lamp^2-2*r_lamp*r(i)*cos(theta)).^(-3/2) )*0.01;
   end

qlamp = 0.70*qlamp; % 0.3 is absorbed by quartz window
qlamp = u*qlamp;    % Operate at a fraction of u

qlamp = Q - scalarfield(S,qlamp);

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