一维稳态导热第二类边界条件的TDMA解法

代码说明：

一根长6cm的金属棒（下面的数据大致为金属铝的数据）水平放置，从左侧端面施加q=10w/(cm)2的热流，右侧端面的对流换热系数h=1000w/(m*m*k)，Tf=20摄氏度，导热系数可以选取lamada=200w/(m*k)。求稳态时各点温度。 采用TDMA算法求解线性方程组。