• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            隨筆-48  評(píng)論-259  文章-1  trackbacks-0

            #include"iostream.h"
            #include"math.h"

            double c[10][10];

            double f(double x)
            {
              double sum=0;
              if(x==0) return 1;
              sum=sin(x)/x;
              return sum;
            }
            void initcotes(double c[][10])
            {
              c[1][0]=c[1][1]=0.5;
              c[2][0]=c[2][2]=1.0/6.0;c[2][1]=2.0/3.0;
              c[3][0]=c[3][3]=1.0/8.0;c[3][1]=c[3][2]=3.0/8.0;
              c[4][0]=c[4][4]=7.0/90.0;c[4][1]=c[4][3]=16.0/45.0;c[4][2]=2.0/15.0;
              c[5][0]=c[4][5]=19.0/288.0;c[5][1]=c[5][4]=25.0/96.0;c[5][2]=c[5][3]=25.0/144.0;
            }
            void Trapezoid(double a,double b)
            {
              cout<<"梯形公式的結(jié)果:"<<(b-a)*(f(a)+f(b))/2<<endl;
            }
            void MidRect(double a,double b)
            {
              cout<<"中矩形公式的結(jié)果:"<<(b-a)*f((b+a)/2)<<endl;
            }
            void NewtonCotes(double a,double b)
            {
              int n,k;double h;
              cout<<"請(qǐng)輸入n的值:";
              cin>>n;
              h=(b-a)/double(n);
              double sum=0;
              for(k=0;k<=n;k++)
               sum+=c[n][k]*f(a+k*h);
              cout<<"牛頓-柯特斯公式的結(jié)果:"<<(b-a)*sum<<endl;
            }
            int STrapezoid(double a,double b)

              int n,k,q;double h;
              cout<<"1--復(fù)化梯形公式"<<endl;
              cout<<"2--復(fù)化辛普森求積公式"<<endl;
              cout<<"輸入你想進(jìn)行的操作:";
              cin>>q; 
              cout<<"請(qǐng)輸入n的值:";
              cin>>n;
              h=(b-a)/double(n);
              double sum=0;
              sum+=(f(a)+f(b));
              for(k=1;k<=n-1;k++)  sum+=2*f(a+k*h);
             if(q==1)
             {
              cout<<"復(fù)化梯形公式的結(jié)果:"<<(h/2)*sum<<endl;
              return 1;
             }


              for(k=0;k<n;k++)
               sum+=4*f(a+(k+0.5)*h);
              cout<<"復(fù)化辛普森求積公式的結(jié)果:"<<(h/6)*sum<<endl;
              return 1;

            }

            void main()
            {
              double a,b;
              int p;
              cout<<"請(qǐng)輸入積分的下、上限:";
              cin>>a>>b;
              initcotes(c);
              while(1)
              {
                cout<<"0--退出"<<endl;
             cout<<"1--梯形公式"<<endl;
             cout<<"2--中矩形公式"<<endl;
             cout<<"3--牛頓柯特斯公式:"<<endl;
             cout<<"4--復(fù)化公式"<<endl;
             cout<<"輸入你想進(jìn)行的操作:";
             cin>>p;
             
             switch(p)
             {
             case 1:Trapezoid(a,b);
             case 2:MidRect(a,b);
             case 3:NewtonCotes(a,b);
                case 4:STrapezoid(a,b);
             }
             if(p==0)  break;
              }
            }


            評(píng)論:
            # re: 數(shù)值分析--數(shù)值積分公式 2007-12-05 10:50 |
            頂!!!!!!!!!!! 太強(qiáng)了   回復(fù)  更多評(píng)論
              
            # re: 數(shù)值分析--數(shù)值積分公式 2008-12-02 11:04 | kissi
            好厲害  回復(fù)  更多評(píng)論
              
            久久精品女人天堂AV麻| 精品久久久一二三区| 婷婷久久香蕉五月综合加勒比| 日韩精品久久久久久久电影| 久久精品国产2020| 一本久久a久久精品综合夜夜| 久久久久亚洲爆乳少妇无| 亚洲精品tv久久久久久久久 | 97久久天天综合色天天综合色hd| 久久精品国产只有精品2020| 无码国内精品久久综合88| 99re久久精品国产首页2020| 亚洲Av无码国产情品久久| 国产成人久久精品麻豆一区 | 狠狠精品久久久无码中文字幕| 欧美久久综合性欧美| 久久精品人人做人人妻人人玩| 亚洲综合久久夜AV | 国产伊人久久| 久久久久久A亚洲欧洲AV冫| 麻豆一区二区99久久久久| 久久久久av无码免费网| 狠狠精品久久久无码中文字幕| 亚洲欧美日韩中文久久| 一级做a爰片久久毛片免费陪 | 日韩精品国产自在久久现线拍| 无码精品久久久久久人妻中字| 欧美精品丝袜久久久中文字幕| 品成人欧美大片久久国产欧美| 久久久久久九九99精品| 亚洲欧美伊人久久综合一区二区 | 一本色综合久久| 久久久久亚洲AV成人网人人软件| 久久电影网2021| 色综合久久综精品| 久久最近最新中文字幕大全| 国产精品久久久久久福利69堂| 久久99国产乱子伦精品免费| 色综合久久综精品| 色综合久久中文色婷婷| 久久精品国产精品亚洲|