传热学非稳态热传导问题网格划分数值计算传热学第四版chapter4 4-6）C++源文件

传热学第四版 高等教育出版社 杨世铭 陶文铨 编著 第四章非稳态热传导例题4-6的数值解法的C++源文件，老师要做的大作业，结果可以做出例题里面的等温线图，给以后要做大作业的同学提供一个参考

#include
#include
#include
using namespace std;
int main（）
{
	ofstream outfile（“4-6.xls“ios::out）;//输出数据文件
	double d=0.1t0=80Ttf=300h=1163na=50a=1.39e-5dxdt=0.18n;
	
	cout<<“Please enter the time（min）:“<	cin>>n;
	int m=11ij;//节点数
	double t[20][2]FoBi;
	dx=d/m;Fo=a*dt/（dx*dx）;Bi=h*dx/na;T=n*60/dt;
	for（i=0;i<=m;i++）
	{
		t[i][0]=t0;
	}
	for（j=0;j<=T;j++）
	{
		for（i=1;i<=m-1;i++）
		{
		t[i][1]=Fo*（t[i+1][0]+t[i-1][0]）+（1-2*Fo）*t[i][0];
		}
    t[0][1]=t[0][0]+2*Fo*（t