青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

隨筆-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)論
  
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            欧美在线关看| 性欧美1819sex性高清| 久热re这里精品视频在线6| 午夜日韩激情| 伊甸园精品99久久久久久| 免费观看日韩av| 你懂的视频一区二区| 日韩午夜电影av| 亚洲午夜国产一区99re久久| 国产美女精品免费电影| 毛片av中文字幕一区二区| 免费久久99精品国产| 国产精品99久久久久久www| 亚洲已满18点击进入久久| 伊人久久婷婷色综合98网| 亚洲国产三级网| 欧美日韩亚洲综合| 久久精品二区三区| 欧美精品色网| 久久久久久国产精品一区| 麻豆九一精品爱看视频在线观看免费| 亚洲精品视频免费| 午夜精品久久久| 亚洲卡通欧美制服中文| 亚洲摸下面视频| 91久久久久久国产精品| 亚洲男人第一av网站| 亚洲国产小视频在线观看| 夜久久久久久| 91久久久亚洲精品| 亚洲欧美中文另类| 一区二区三区高清不卡| 久久精品视频网| 亚洲专区一区| 欧美精品成人| 免费一区二区三区| 国产欧美日韩视频在线观看| 亚洲激情av| 在线观看一区二区视频| 亚洲你懂的在线视频| 99视频一区二区| 久久久久久综合网天天| 午夜日韩av| 欧美天天在线| 亚洲人成小说网站色在线| 影音先锋欧美精品| 亚洲欧美日韩区| 亚洲免费综合| 欧美日韩精品二区| 亚洲国产综合91精品麻豆| 国外视频精品毛片| 欧美一区二区三区的| 亚洲欧美成人| 欧美日韩国产在线看| 欧美福利网址| 在线观看亚洲精品| 久久婷婷麻豆| 美女脱光内衣内裤视频久久影院| 国产精品综合久久久| 亚洲午夜精品福利| 亚洲综合首页| 国产精品久久久久9999吃药| 日韩一区二区免费看| 一本色道综合亚洲| 欧美日韩中文另类| 亚洲欧美bt| 久久先锋影音av| 国产日韩精品在线观看| 午夜精品视频| 久久免费国产精品| 在线观看日韩av| 久久阴道视频| 最新国产拍偷乱拍精品| 日韩视频免费观看| 欧美视频在线观看 亚洲欧| 一区二区激情| 久久国产福利国产秒拍| 国内精品**久久毛片app| 久久久91精品国产| 欧美激情精品久久久久久大尺度 | 篠田优中文在线播放第一区| 亚洲欧美影院| 国产区欧美区日韩区| 久久久精品动漫| 亚洲国产精品一区二区三区| 一区二区三区高清| 国产日韩三区| 久久综合网络一区二区| 亚洲激情欧美| 午夜精品一区二区三区在线视| 国产亚洲精品自拍| 欧美bbbxxxxx| 亚洲性视频h| 老司机精品福利视频| 99re6热只有精品免费观看| 欧美视频观看一区| 久久精品免费电影| 日韩亚洲欧美成人一区| 久久精品亚洲一区| 亚洲美女一区| 国产日本欧美一区二区| 欧美ab在线视频| 亚洲综合999| 亚洲国产高清视频| 久久成年人视频| 亚洲日韩视频| 国产一区二区三区高清播放| 欧美激情中文字幕一区二区| 先锋影音国产一区| 日韩午夜av电影| 免费一级欧美片在线观看| 亚洲欧美www| 日韩视频欧美视频| 激情久久五月天| 国产精品久久久对白| 免费亚洲电影在线观看| 午夜精品一区二区三区在线| 亚洲区一区二区三区| 久久一区二区三区超碰国产精品| 亚洲午夜视频| 日韩视频一区二区三区在线播放免费观看| 国产毛片久久| 国产精品igao视频网网址不卡日韩 | 国产视频精品xxxx| 国产精品第一页第二页第三页| 看片网站欧美日韩| 久久成人精品| 亚洲综合另类| 亚洲男女毛片无遮挡| 一本色道久久综合狠狠躁篇怎么玩| 欧美成人免费小视频| 麻豆国产精品777777在线| 欧美专区在线播放| 欧美亚洲三区| 亚洲综合色自拍一区| 亚洲视频在线观看网站| 国产精品久久二区| 欧美日韩精品一区二区天天拍小说 | 欧美精品v日韩精品v国产精品| 久久久久久久91| 欧美在线一二三四区| 香蕉久久夜色精品| 亚洲欧美一区二区视频| 亚洲一二三区视频在线观看| 一区二区三区四区五区精品视频| 亚洲日本中文字幕| 亚洲乱码日产精品bd| 亚洲精品色婷婷福利天堂| 91久久精品一区二区别| 亚洲欧洲精品一区| 日韩视频免费观看高清在线视频| 亚洲精品一二| 亚洲视频二区| 午夜精品美女久久久久av福利| 亚洲综合视频一区| 久久爱www久久做| 久久天天躁狠狠躁夜夜爽蜜月 | 亚洲欧美国产制服动漫| 亚洲免费视频一区二区| 欧美亚洲综合在线| 久久久久久久一区二区| 免费观看久久久4p| 欧美日韩三级在线| 国产精品免费一区豆花| 国产午夜精品理论片a级探花| 好看的av在线不卡观看| 亚洲国产日韩一区二区| 99在线观看免费视频精品观看| 亚洲午夜高清视频| 久久天堂成人| 亚洲三级电影在线观看| 亚洲视频在线一区| 久久黄金**| 欧美国产日本在线| 国产精品三级视频| 亚洲第一网站免费视频| 一区二区日本视频| 久久夜色撩人精品| 亚洲免费观看| 久久经典综合| 欧美视频导航| 在线欧美视频| 亚洲欧美日韩在线一区| 免费在线成人av| 一区二区三区日韩欧美| 久久婷婷国产综合精品青草| 欧美日韩一区在线观看| 精品9999| 香蕉成人伊视频在线观看| 欧美电影免费观看| 亚洲欧美日韩成人| 欧美日本国产| 亚洲高清免费视频| 性色av一区二区怡红| 亚洲精选91| 欧美aa国产视频| 伊人狠狠色j香婷婷综合| 西瓜成人精品人成网站| 亚洲毛片网站| 欧美成人综合|