鱼C论坛

 找回密码
 立即注册
查看: 2531|回复: 5

大家帮我找一下这个程序的错误,好纠结啊 ~

[复制链接]
发表于 2011-11-27 17:14:25 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
#include<stdio.h>
#include<string.h>
#include<math.h>
int SGN(int n)
        {
      int m;
      if(n>=0)
                  m=1;
      else
                  m=-1;
      return m;
        }

main()
{   
         
         int i,j,K,k,N=6,M=17,ci=0;
     int L[17],IO[17],JO[17];
     double EP,FFN,D[17],Q[17],S[17],H[17];
     double  F[6],R[6],DQ[6];
     printf("Input EP and FFN:\n\n");
     scanf("%f%f",&EP,&FFN);
     printf("\n Input K L(m) D(m) Q(L/s) IO JO:\n");
     for(i=0;i<M;i++)

     scanf("%d%d%f%f%d%d",&j,&L[i],&D[i],&Q[i],&IO[i],&JO[i]);
     for(i=0;i<M;i++)
       Q[i]=Q[i]/1000;
     for(i=0;i<M;i++)
         {
       S[i]=10.293*FFN*FFN*L[i]/pow(D[i],5.33);
       Q[i]=Q[i]*SGN(IO[i]);
         }
     do
         {
       ci=ci+1;
       for(i=0;i<N;i++)
       H[i]=S[i]*fabs(Q[i])*Q[i];
       for(i=0;i<N;i++)
           {
         F[i]=0;
         R[i]=0;
         for(j=0;j<M;j++)
                 {
           if(fabs(IO[j])-1==i)
                   {
                           F[i]=F[i]+H[j];
               R[i]=R[i]+S[j]*fabs(Q[j]);
                   }
           else
                   {
               if(fabs(JO[j])-1==i)
                           {
                                   F[i]=F[i]+H[j]*SGN(JO[j]);
                       R[i]=R[i]+S[j]*fabs(Q[j]);
                           }
                   else;
                   }
                }
        DQ[i]=-F[i]/R[i]/2;
          }

    for(i=0;i<M;i++)
        {
                if(JO[i]!=0)
        Q[i]=Q[i]+DQ[fabs(IO[i])-1]-DQ[fabs(JO[i])-1];
        else
                        Q[i]=Q[i]+DQ[fabs(IO[i])-1];
                H[i]=S[i]*fabs(Q[i])*Q[i];
        }
    k=0;
    for(i=0;i<N;i++)
        {
                if(fabs(F[i])<=EP);
        else
                        k=1;
        }

    while(k==1);
    printf("\n\nResults:\n");
    printf("%32s%.3f","EP=",EP);  
        printf("%10s%d\n","ci=",ci);
    printf("%33s%d","N=",N);      
        printf("%11s%d\n","M=",M);
    printf("\n");
    printf("%7s%13s%13s%10s%14s%15s\n","K","L(m)","D(m)","S","H(m)","Q(m^3/s)");
    for(i=0;i<M;i++)
    printf("%7d%13d%13.2f%13.3E%11.2f%14.5f\n",i+1,L[i],D[i],S[i],H[i],Q[i]);
}
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
发表于 2011-11-29 00:41:09 | 显示全部楼层
你的程序想要实现什么,具体思路是什么,以及请在适当的地方作注示。

看这么长的代码才叫人纠结呢。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
发表于 2011-11-29 02:16:50 | 显示全部楼层
楼主,你一个注释没有,你当我们是计算机啊,能够解析代码么..看着眼睛疼.
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
发表于 2011-11-30 16:03:49 | 显示全部楼层
花了二十分钟帮你整理格式,一分钟帮你解决两个问题:
1. 少了一个右大括号
2. fabs()返回double型,你用来做数组下标肯定会出错的。

请以后注意代码格式,因为很少有人喜欢帮别人整理内衣裤的。。。。。。


  1. #include<stdio.h>
  2. #include<string.h>
  3. #include<math.h>

  4. int SGN(int n)
  5. {
  6.         int m;
  7.         if(n>=0)
  8.                 m=1;
  9.         else
  10.                 m=-1;
  11.         return m;
  12. }

  13. int main()
  14. {   
  15.         int i,j,K,k,N=6,M=17,ci=0;
  16.         int L[17],IO[17],JO[17];
  17.         double EP,FFN,D[17],Q[17],S[17],H[17];
  18.         double F[6],R[6],DQ[6];
  19.        
  20.         printf("Input EP and FFN:\n\n");
  21.         scanf("%f%f",&EP,&FFN);
  22.         printf("\n Input K L(m) D(m) Q(L/s) IO JO:\n");
  23.        
  24.         for(i=0;i<M;i++)
  25.                 scanf("%d%d%f%f%d%d",&j,&L[i],&D[i],&Q[i],&IO[i],&JO[i]);
  26.        
  27.         for(i=0;i<M;i++)
  28.                 Q[i]=Q[i]/1000;

  29.         for(i=0;i<M;i++)
  30.         {
  31.                 S[i]=10.293*FFN*FFN*L[i]/pow(D[i],5.33);
  32.                 Q[i]=Q[i]*SGN(IO[i]);
  33.         }
  34.        
  35.         do
  36.         {
  37.                 ci=ci+1;
  38.                 for(i=0;i<N;i++)
  39.                         H[i]=S[i]*fabs(Q[i])*Q[i];
  40.                 for(i=0;i<N;i++)
  41.                 {
  42.                         F[i]=0;
  43.                         R[i]=0;
  44.                         for(j=0;j<M;j++)
  45.                      {
  46.                                 if(fabs(IO[j])-1==i)
  47.                                {
  48.                                         F[i]=F[i]+H[j];
  49.                                            R[i]=R[i]+S[j]*fabs(Q[j]);
  50.                              }
  51.                                 else
  52.                                 {
  53.                                            if(fabs(JO[j])-1==i)
  54.                                         {
  55.                                                 F[i]=F[i]+H[j]*SGN(JO[j]);
  56.                                            R[i]=R[i]+S[j]*fabs(Q[j]);
  57.                                }
  58.                                }
  59.                     }
  60.                         DQ[i]=-F[i]/R[i]/2;
  61.                 }
  62.                 for(i=0;i<M;i++)
  63.                 {
  64.                         if(JO[i]!=0)
  65.                                 Q[i]=Q[i]+DQ[(int)fabs(IO[i])-1]-DQ[(int)fabs(JO[i])-1];
  66.                         else
  67.                                 Q[i]=Q[i]+DQ[(int)fabs(IO[i])-1];
  68.                         H[i]=S[i]*fabs(Q[i])*Q[i];
  69.                 }
  70.                 k=0;
  71.                 for(i=0;i<N;i++)
  72.                 {
  73.                         if(fabs(F[i])<=EP)
  74.                                 ;
  75.                         else
  76.                                 k=1;
  77.                 }
  78.         }while(k==1);
  79.                
  80.         printf("\n\nResults:\n");
  81.         printf("%32s%.3f","EP=",EP);  
  82.         printf("%10s%d\n","ci=",ci);
  83.         printf("%33s%d","N=",N);      
  84.         printf("%11s%d\n","M=",M);
  85.         printf("\n");
  86.         printf("%7s%13s%13s%10s%14s%15s\n","K","L(m)","D(m)","S","H(m)","Q(m^3/s)");

  87.         for(i=0;i<M;i++)
  88.                 printf("%7d%13d%13.2f%13.3f%11.2f%14.5f\n",i+1,L[i],D[i],S[i],H[i],Q[i]);
  89.                
  90.         return 0;
  91. }
复制代码



想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
 楼主| 发表于 2011-11-30 16:05:46 | 显示全部楼层
哈哈哈,谢了哈,甲鱼哥~
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
头像被屏蔽
发表于 2011-12-2 19:40:41 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-3-29 22:34

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表