{$N+}PROGRAM GINT;

CONST
	pi=3.141592654;
	m=15;

TYPE
	pole=ARRAY[1..m] OF double;

FUNCTION f(x:double):double;
BEGIN
	f:=sin(x)*exp(cos(x));
        {f:=x*x*x*x*ln(10+abs(x))*cos(x)*exp(-abs(x));}
END;


VAR
	a,b,integral:double;
	body,vahy:pole;
	i:integer;
        filein:text;
BEGIN
	a:=0;
	b:=pi;
        integral:=0;

        ASSIGN(filein,'soubor.dat');
        RESET(filein);
        
        FOR i:=1 TO m DO BEGIN
                READ(filein,body[i]);
                READ(filein,vahy[i]);
        END;
        
	FOR i:=1 TO m DO BEGIN
                integral:=integral+vahy[i]*f(((b-a)*body[i]+b+a)/2);
	END;
        integral:=integral*(b-a)/2;

        WRITELN('Integral z funkce sin(x)*exp(cos(x)) na intervalu <',a:2:1,',',b:2:1,'> je ',integral:8:7);
        CLOSE(filein);
        READLN;
END.