高手帮忙看一下--结果输出问题

li343

程序为何解不出t[n],v[n]和m[n]的值?
程序如下:

1. #include<stdio.h>
   
2. #include<iostream.h>
   
3. #include<math.h>
   
4. void main()
   
5. {
   
6. int i,j;
   
7. double p[100],t[100],m[100],v[100]; // p为密度，t为温度，m为质量流量
   
8. double b[13]={20,25,30,35,40,45,50,55,60,65,70,75,80},G[13]={697.78,705.04,708.51,708.15,703.98,696.02,684.33,669.03,650.23,628.11,602.86,574.72,543.96}; //输入角度和G
   
9.     double cp=1005,ta=293,pa=1.205,ycol=0.6,a=1378,Hchi=123,h=10,L=1501,Achi=2290.2;
   
10. double f=0.008428,ein=0.056;  //损失系数
    
11.   //double f=0.056,ein=1.5;  //损失系数
    
12. double D=54;  //D为集热棚的水力直径
    
13. double Acol=0.5*a*a;
    
14. for(j=0;j<13;j++)
    
15. {
    
16.   //t[0]=300,p[0]=1.0,m[0]=10000,v[0]=1.0; //赋初值
    
17.   t[0]=0,p[0]=1.0,m[0]=10000,v[0]=1.0; //赋初值
    
18.   int n=0;
    
19.   for(i=0;1;i++)
    
20.   {
    
21.   //```````````````````````方程的迭代格式``````````````````````````````````````````````//
    
22.             t[i+1]=(ycol*Acol/cp)*G[j]/m+ta;
    
23.             p[i+1]=2.2204-0.00353*t;
    
24.   
    
25.             m[i+1]=p*v*Achi;
    
26.       v[i+1]=sqrt((2*(pa-p)*9.8*(Hchi+0.5*a*sin(b[j]*3.1415926/180))-0.5*f*L*p*pow(v,2)/D-0.5*ein*pa*pow((m/a/h),2))/p);
    
27.   //```````````````````````````````````````````````````````````````````````````````````//
    
28.    if((fabs(t[i+1]-t)<10e-6) && (fabs(p[i+1]-p)<10e-6)&&(fabs(m[i+1]-m)<10e-6)&&(fabs(v[i+1]-m)<10e-6))// 判断跌代的值
    
29.     break;  
    
30.    n++;
    
31.   }
    
32. printf("%f\t",p[n]);
    
33. printf("%f\t",t[n]);
    
34.     printf("%f\t",v[n]);
    
35. printf("%f\t",m[n]);
    
36. 
    
37. printf("\n");
    
38. }
    
39. 
    
40. }

复制代码

[ 本帖最后由 风花雪月 于 2008-4-9 08:24 编辑 ]

风花雪月

是否有什么提示？

cooller

你的程序在vc中错误很多。我想帮助你，但是太多错误了。你先调试一下吧。