‘äŒ`‘¥
<![CDATA[
program daikeisoku(input,output);
const a=0;
b=1;
pi=3.1415;
var I1,I2,h:real;
N,x,i:integer;
function f1(x:real):real;
begin
f1:=4/(1+sqr(x))
end; { f1 }
function f2(N:integer):real;
begin
h:=(b-a)/N;
I1:=(f1(a)+f1(b))/2;
for x:=1 to N-1 do I1:=I1+f1(a+x*h);
f2:=h*I1
end; { f2 }
function g1(x:real):real;
begin
x:=sin(2*pi*x);
g1:=exp(x)
end; { g1 }
function g2(N:integer):real;
begin
h:=(b-a)/N;
I2:=(g1(a)+g1(b))/2;
for x:=1 to N-1 do I2:=I2+g1(a+x*h);
g2:=h*I2
end; { g2 }
begin
N:=2;
while N<=16 do
begin
N:=2*N;
writeln('N=',N:1,' ','I1=',f2(N),' ','I2=',g2(N))
end;
end.
]]>
Last modified: 2000$BG/(B1$B7n(B20$BF|(B 13:51