このディレクトリの索引
http://hibari.2ch.net/test/read.cgi/tech/1291471791/483
#  [1] 授業単元:数値計算 
#  [2] ∫_0^π?sinxdxにおいて2m=4とし、 10分割した際の数値積分をシンプソンの公式により求めなさい 
#  

f(X,Y) :- Y is sin(X).

'10分割した際の数値積分をシンプソンの公式により求める'(_定積分の近似値) :-
        シンプソンの公式(10,0,pi,_定積分の近似値).

シンプソンの公式(_分割数,_a,_b,_定積分の近似値) :-
        _h is ( _b - _a) / (2 * _分割数),
        _m is _n - 1,
        findsum(S,(
                    for(0,_i,_m),
                    f(_a+(2*_i)*_h,Y1),
                    f(_a+(2*_i+1)*_h,Y2),
                    f(_a+(2*_i+2)*_h,Y3),
                    S is Y1 + 4 * Y2 + Y3),
                Sum),
        _定積分の近似値 is (_h / 3) * Sum.