二维稳态传热.cpp

传热学计算类程序

#include<iostream.h>
#include<math.h>
double x,y,s;
void main(void)
{
	double l,w,T1,T2,t,N,M,h;        // h,tw,chuanrexishu,houdu,gaodu,n,l,th,p,A,x;
    cin>>l>>w>>T1>>T2>>N>>M>>h;	
	 	for(y=0;y<=w;y+=w/M) 
			for(x=0;x<=l;x+=l/N)
			{
				s=0;
				for(double n=1;n<10000;n++){
	if(n%2==1)s=s+2/double(n)*sin(n*3.1415926*x/l)*sin(h*n*3.1415926*y/l)/sin(h*n*3.1415926*w/l);
				}
				
				t=T1+(T2-T1)*2*s/3.1415926;
	if(x>0&&y>0)cout<<"x="<<x<<"\t"<<"y="<<y<<"\t"<<"t="<<t<<"\n";
	}
	
	}