[AS3]辛普森数值积分

/**
* 在[a, b]间的辛普森数值积分
* @param f 积分函数
* @param a 左界限
* @param b 右界限
* @param n 分段数
* @return [a, b]间的定积分
*/

public static function simpsonInte(f:Function, a:Number, b:Number, n:int = 1024):Number {
 n*=2;//2n => n
 var h:Number = (b-a) / (n);
 var x1:Number = 0;
 var x2:Number = 0;
 n--;//2n-1 => n
 for(var i:int = 1; i i++*/) {
  x1 += f(a + i*h);
  i++;//i+1 => i
  x2 += f(a + i*h);
 }
 return (h/3)*(f(a) + f(b) + 4*x1 + 2*x2);
}

本文链接:https://assnippets.blogspot.com/2008/01/as3_30.html转载请注明出处。

相关文章：