中点公式.pptVIP

  1. 1、本文档共27页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
中点公式

* * Quadrature 二阶: 中点公式: 梯形公式: 四阶公式: Simpson’s 1/3rd Rule Boole’s rule,The 6-th Newton-Cotes rule (the first step of Romberg integration) The extrapolated Simpson’s rule. #include stdio.h #include stdlib.h #include math.h #define TRUE 1 struct t_bc{float a, b, h;} bc; void main() { int isimp, k, n; float s; void simps(), trapz(); printf( \nComputer Soft/C7-1 Trapezoidal/Simpsons Rule \n\n ); while( TRUE ){ printf( Type 0 for trapezoidal, or 1 for Simpsons\n ); scanf( %d, isimp ); printf( Number of intervals ? ); scanf( %d, n ); printf( Lower limit of integration? ); scanf( %f, bc.a ); printf( Upper limit of integration? ); scanf( %f, bc.b ); bc.h = (bc.b - bc.a)/n; if( isimp == 0 ){ trapz( s, n ); /*-- Trapezoidal rule */ } else{ simps( s, n ); /*-- Simpsons rule */ } printf( -----------------------------------------\n ); printf( Result = %g \n, s ); printf( -----------------------------------------\n ); printf( \nType 1 to continue, or 0 to stop.\n); scanf( %d, k ); if( k != 1 ) exit(0); } } void simps(ss, n) /* Simpsons rule*/ float *ss; int n; { int i, ls; float sum, s, w, x; double func(); s = 0; sum = 0; if( n/2*2 == n ) ls = 0; else { ls = 3; for( i = 0; i = 3; i++ ) { /* Simpsons 3/8 rule if n is odd */ x = bc.a + bc.h*i; w = 3; if( i == 0 || i == 3 ) w = 1; sum = sum + w*func( x ); } sum = sum*bc.h*3/8L; if( n == 3 ) return; } for( i = 0; i = (n - ls); i++ ){ /* Simpsons 1/3 rule */ x = bc.a + bc.h*(i + ls); w = 2; if( (int)( i/2 )*2 + 1 == i ) w = 4; if( i == 0 || i == n - ls ) w = 1; s = s + w*func( x ); } *ss = sum + s*bc.h/3; return; } void trapz(ss, n) /* Trapezoidal rule */ float *ss; int

文档评论(0)

zijingling + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档