标题:C++语言编程,用定步长辛普森公式求积分
只看楼主
zoeyjia
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2013-10-7
结帖率:0
 问题点数:0 回复次数:2 
C++语言编程,用定步长辛普森公式求积分
请问下面的程序哪里错了,急需各位大神帮忙指出修正。程序题目见附件。急!急!急!!!小女子在此不胜感激
C++实验 题目.rar (7.54 KB)

#include <iostream.h>
#include <math.h>
double f(double x)
{
    double t;
    t=cos(x);
    return t;
}
void main()
{
    int n=1,i;
    float m,s,s1,s2=0,s3,h,a=0,b=0.8,c=0.5,p,q,d;
    double eps=1e-8;
    while(fabs(f(n)-f(n-1))>=eps)
    {
        h=(b-a)/(2*n);
        m=(b-a)/(3*n);
        c=0.5;
        d=f(a)-f(b);
        s1=c*d;
        for(i=1;i<=n;i++)
        {
            p=2*f(a+(2*i-1)*h);
            q=f(a+2*i*h);
            s3=p+q;
            s2=s2+s3;
            s1=m*(s1+s2);
        }
        cout<<s<<endl;
    }
}
搜索更多相关主题的帖子: 语言编程 不胜感激 include 辛普森 double 
2013-10-08 10:32
在这里爬起
Rank: 3Rank: 3
等 级:论坛游侠
帖 子:43
专家分:113
注 册:2013-8-9
得分:0 
程序代码:
#include <iostream.h>
#include <math.h>
#define eps 1E-8
double f(double x)
{
    double t;
    t=cos(x);
    return t;
}
void main()
{
    int n,i;
    double m,s1,s3,h,a=0,b=0.8,c1,c=0.5,p,q,d;
        for(n=1;;n++)
        {
          double s2=0,s;
          m=(b-a)/(3*n);
          h=(b-a)/(2*n);
          d=f(a)-f(b);
          c1=c*d;
            for(i=1;i<=n;i++)
            {
            p=2*f(a+(2*i-1)*h);    
            q=f(a+2*i*h);
            s3=p+q;
            s2=s2+s3;
            }
          s=m*(c1+s2);
          s1=((b-a)/(3*n*(n-1)))*(2*f(a+(2*n*h)-h)+f(a+2*n*h));
         if(fabs(s1)<eps)
         {
           cout<<s<<endl;   
              break;
         }
        }   
}
我给你改了下。你自己看看
2013-10-12 13:23
zoeyjia
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2013-10-7
得分:0 
回复 2楼 在这里爬起
多谢啦~~
2013-10-15 10:12



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-421561-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.021653 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved