Const
 a=0.0;
 b=3.0;
 h=0.5;

Function F(x:Real):Real;
begin
 F:=Cos(exp(x/3)+x);
end;

Procedure Rectangle;
var
 i,N:Integer;
 Sum:Real;
begin
 N:=Round((b-a)/h);
 Sum:=0;
 for i:=0 to N-1 do
  Sum:=Sum+h*F(a+h*i);
 WriteLn('Rectangle:     ', Sum:0:5);
end;

Procedure Trapezium;
var
 i,N:Integer;
 Sum:Real;
begin
 N:=Round((b-a)/h);
 Sum:=(F(a)+F(b))/2*h;
 for i:=1 to N-1 do
  Sum:=Sum+h*F(a+h*i);
 WriteLn('Trapezium:     ', Sum:0:5);
end;

Procedure Simpson;
var
 i,N:Integer;
 Sum,hs:Real;
begin
 hs:=h;
 N:=Round((b-a)/hs) div 2;
 If (N mod 2)=1 then
  begin
   N:=N*2;
   hs:=hs/2;
  end;
 Sum:=(F(a)+F(b))/3*hs;
 for i:=1 to N do
  Sum:=Sum+4.0*hs/3.0*F(a+hs*(2*i-1));
 for i:=2 to N do
  Sum:=Sum+2.0*hs/3.0*F(a+hs*(2*i-2));
 WriteLn('Simpson:       ', Sum:0:5);
end;

Begin
 Rectangle;
 Trapezium;
 Simpson;
 ReadLn;
End.