【文件属性】:
文件名称:一维非定常热传导方程
文件大小:245B
文件格式:CPP
更新时间:2015-12-10 17:03:15
一维非定常
//显式计算一维非定常热传导方程
#include
#include
using namespace std;
int main()
{
double T[11]; //存放初始值和边界值
double result[11]; //存放结果
double L=1; //几何模型
double deltaX=0.1; //空间步长
double deltaT=0.001; //时间步长
int n=L/deltaX+1; //空间节点数
for(int i=1;i<10;i++)//初始化
{
if(abs((i+1)*deltaX) - 0.5*L < 1e-9)
{
T[i]=2*(i)*deltaX;
}
else
{
T[i]=2*(1-(i)*deltaX);
}
T[0]=0;
T[n]=0;
}
//显示初始数据
cout<< "time:" << 0*deltaT <