DEVOIR PRATIQUE BAC 2013 G1

program aires;
uses wincrt;
type
enreg= record
rect,trapez:real;
nombre:integer;
end;
fenreg=file of enreg;
var
f:fenreg;
eps:real;
resultat:integer;
function fct (x:real):real;
begin
fct:=sqr(x);
end;
procedure saisirepsilon (var eps:real);
begin
repeat
writeln('entrer epsilon');
readln(eps);
until (eps>=0.001) and (eps<=0.1);
end;
function rectangle (a,b:real; n:integer):real;
var
k,s:real;
i:integer;
begin
s:=0;
k:=(b-a)/n;
for i:= 1 to n do
begin
s:=s+ fct(a+k/2)*k;
a:=a+k;
end;
rectangle:=s;
end;
function trapeze (a,b:real; n:integer):real;
var
i:integer;
s,k:real;
begin
s:=0;
k:=(b-a)/n;
for i:= 1 to n do
begin
s:=s+ (fct(a)+fct(a+k))*k/2;
a:=a+k;
end;
trapeze:=s;
end;
procedure traitement (var f:fenreg; eps:real; var resultat:integer);
var
n:integer;
e:enreg;
k,j,rec,trap:real;
begin
rewrite(f);
n:=1;
k:=0;
j:=0;
rec:=1;
trap:=1;
while(abs(rec-k)>eps)and(abs(trap-j)>eps)do
k:=rec;
j:=trap;
n:=n+1;
rec:=rectangle(0,3,n);
trap:=trapeze(0,3,n);
e.nombre:=n;
e.rect:=rec;
e.trapez:=trap;
write(f,e);
writeln('n=',n,' | ', rec:9:5,' | ', trap:9:5);
if(abs(rec-k)resultat:=1
else if (abs(rec-k)>abs(trap-j)) then
resultat:=2
else
resultat:=3;
writeln('la methode est ',resultat);
close(f);
end;
procedure afficher (var f:fenreg;resultat:integer);
var
e:enreg;
begin
writeln('| N |','| Rectangle |','| trapeze |');
reset(f);
while not(eof(f))do
begin
read(f,e);
writeln('|',e.nombre,'|','| ',e.rect:9:5,' |','| ',e.trapez:9:5,' |');
end;
case resultat of
1:writeln('la methode qui converge la premiere vers l''aire exacte est la methode de rectangle');
2:writeln('la methode qui converge la premiere vers l''aire exacte est la methode de trapeze');
3:writeln('les deux methodes convergent en meme temp vers 9 qui est l''aire exacte');
end;
end;
begin
assign(f,'d:\calcul.dat');
saisirepsilon(eps);
traitement(f,eps,resultat);
afficher(f,resultat);
end.

La méthode des trapèzes :
s:=s+(fct(a)+fct(a+h))*k/2;

Il y a une ligne incomplète :

if(abs(rec-k)