• <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)論
              
            日韩欧美亚洲综合久久影院Ds| 国产精品美女久久久久久2018| 久久综合欧美成人| 久久99国产精品久久99| 久久精品无码一区二区三区免费| 久久久久亚洲av成人无码电影| 一本久久知道综合久久| 亚洲一区中文字幕久久| 久久精品国产2020| 色综合合久久天天综合绕视看| 亚洲精品国产自在久久| 国产V综合V亚洲欧美久久| 久久久久国产精品麻豆AR影院 | 久久人人爽人人爽人人爽| 亚洲国产精品无码久久SM| 国产精品九九久久免费视频| 久久久亚洲裙底偷窥综合 | 88久久精品无码一区二区毛片 | 久久99国产精品二区不卡| 亚洲七七久久精品中文国产 | 久久综合亚洲色HEZYO社区| 狠狠色丁香婷婷久久综合 | 国内精品久久久久国产盗摄| 香蕉99久久国产综合精品宅男自 | 国产精品熟女福利久久AV| 伊人久久精品无码av一区| 欧美日韩精品久久久久| 亚洲天堂久久精品| 99久久久精品免费观看国产| 久久久久一本毛久久久| 91精品免费久久久久久久久| 精品免费久久久久久久| 婷婷伊人久久大香线蕉AV | 久久天堂电影网| 国产成人久久精品一区二区三区| 久久久久亚洲AV无码观看 | 99久久人人爽亚洲精品美女 | 国内精品伊人久久久久妇| 久久久精品国产| 精品久久久久久中文字幕大豆网| 久久这里都是精品|