residual.m

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

function A = residual(A)

unpack(A)

% 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);

% 0: AE, 1: ODE
A = set(A,'dxdtcoeff',{[0 0 0 0]});