function A=secdim2
% 本程序可实现两种可选边界条件的对流扩散实时动态模拟
% 采用乘方格式
% 初始温度场为 10，速度可选
% 时间差分采用全隐式格式
clc
dt=0.001;
bc=input('南北边界条件：1-第一类，2-第二类: ');
q=0;
ddt=input('迭代次数： ');
L1=0.02;r=1000;k=0.017;sp=0;L2=0.01;
u=input('横向速度： ');
if bc==2
v=0;
else
v=input('纵向速度 ： ');
end
n=100;m=100;
dx=L1/n;dy=L2/m;F1=r*u*dy;F2=r*v*dx;
ap0=r*dx*dy/dt;
D1=k*dy/dx;pe1=F1/D1;D2=k*dx/dy;pe2=F2/D2;
tp0=10*ones(n*m,1);%初场温度
