5.二维显式导热程序（存在问题，没找到）

//存在一些问题
#include <iostream>

	double length,height;  //导热体的长度,高度
	double TL,TR,TU,TD;  //左侧温度,右侧温度
	double den,cv,k;  //物质密度,比热,导热系数
	double s;  //热源
	double dt;  //时间间隔

    int n,m;   //x方向和y方向控制体的个数
	int tmaxstep; 

	int tr;  //0,代表显示计算；1代表隐式计算
//-------------以上为和物理问题有关的变量--------------------

	double dx,dy;
	double T0[103][103],T[103][103];  
	double ae0[103][103],aw0[103][103],an0[103][103],as0[103][103],ap0[103][103],ap1[103][103],b[103][103];
    double ae1[103][103],aw1[103][103];

//double a[103],b0[103],c0[103],d[103],x[103],p[103],q[103];  //前三为系数矩阵系数，d为等式右侧矩阵系数,x为未知数
//三对角子函数
/*
void TDMA()  
{
	int i;
	i=1;
		p[i]=-c0[i]/a[i];
		q[i]=d[i]/a[i];
	for(i=2;i<n;i++)
	{
		p[i]=-c0[i]/(a[i]+b0[i]*p[i-1]);
		q[i]=(d[i]-b0[i]*q[i-1])/(a[i]+b0[i]*p[i-1]);
	}
	i=n;
		x[i]=(d[i]-b0[i]*q[i-1])/(a[i]+b0[i]*p[i-1]);
	for(i=n-1;i>=1;i--)
	{
		x[i]=p[i]*x[i+1]+