满意答案
wtweixiao
2013.11.13
采纳率:49% 等级:12
已帮助:11058人
给你个参考,当时我们上数值传热学的时候老师叫我们做的一个平板非稳态导热
不懂再问
#include
main()
{float a[11],t[11],y,f=34.89/(78*7.12),e=2.33/34.89;
int i,j;
long int k;
clrscr();
for(i=0;i<11;i++)
a[i]=1000;
for(k=1;;k++)
{y=0;
for(i=1;i<11;i++)
if(i!=10) t[i]=(f*a[i+1]+f*a[i-1]+(1-2*f)*a[i]);
else t[i]=(1-2*f-2*f*e)*a[i]+2*f*a[i-1]+2*f*e*20;
t[0]=a[0]*(1-2*f)+2*f*a[1];
for(i=0;i<11;i++)
{if(y
a[i]=t[i];}
if(k<=10) printf("The %lds temperature is:%.2f\n",k,a[10]);
if(y<=1e-6) break; }
printf("the total times:%ld.\t t=%.2f",k,a[10]); }
50分享举报