- 541 名前：デフォルトの名無しさん mailto:sage [2009/07/22(水) 14:30:31 ]
- >>540 (1)
% Prolog
f(X,Y) :- Y is 1 / (1 + x ^ 2).
シンプソンの公式(A,B,M,X) :-
D is (B - A) / (2 * M),
f(A,S),
N is 2 * M - 1,
シンプソン公式の二(1,N,A,D,S,S2),
f(B,Z),
X is (S2 + Z) * D / 3.
シンプソンの公式の二(M,N,_,_,X,X) :- M > N,!.
シンプソンの公式の二(M,N,A,D,Y,X) :-
U is A + M * D,
Mod is M mod 2,
シンプソンの公式の三(Mod,U,Y,Z),
M2 is M + 1,
シンプソンの公式の二(M2,N,A,D,Z,X).
シンプソンの公式の三(0,X,S,S2) :- f(X,Y),S2 is S + 2 * Y.
シンプソンの公式の三(1,X,S,S2) :- f(X,Y),S2 is S + 4 * Y.
