{$N+}PROGRAM MC;

CONST
	pi=3.141592654;
	kroku=5000;
        range=65000;

FUNCTION f(x,y:double):double;
BEGIN
	f:=abs(sin(sqr(pi)*(sqr(x-0.5)+sqr(y-0.5))));
END;

FUNCTION rnd:double;
BEGIN
     rnd:=RANDOM(range)/range;
END;

VAR
	i,uspech:longint;
	x,y,z,integral:double;

BEGIN
	RANDOMIZE;
	uspech:=0;

	{Integral z funkce f v intervalu <0,1> x <0,1> prvni metodou}

	FOR i:=1 TO kroku DO BEGIN
		x:=rnd;
		y:=rnd;
		z:=rnd;
		IF (z<f(x,y)) THEN uspech:=uspech+1;
	END;
	integral:=uspech/kroku;

	WRITELN('Monte Carlo integral prvni metodou ',integral:5:4);

	{Integral z funkce f v intervalu <0,1> x <0,1> druhou metodou}

	z:=0;

	FOR i:=1 TO kroku DO BEGIN
		x:=rnd;
		y:=rnd;
		z:=z+f(x,y);
	END;
	integral:=z/kroku;

	WRITELN('Monte Carlo integral druhou metodou ',integral:5:4);

	{Integral z funkce f uvnitr kruhu o polomeru 0.5 se stredem v (0.5,0.5) prvni metodou}

	uspech:=0;
	FOR i:=1 TO kroku DO BEGIN
		x:=rnd;
		y:=rnd;
		IF (sqr(x-0.5)+sqr(y-0.5))<sqr(0.5) THEN BEGIN
			z:=rnd;
			IF (z<f(x,y)) THEN uspech:=uspech+1;
		END;
	END;
	integral:=uspech/kroku;

	WRITELN('Monte Carlo integral v kruhu prvni metodou ',integral:5:4);

	{Integral z funkce f uvnitr kruhu o polomeru 0.5 se stredem v (0.5,0.5) druhou metodou}

	z:=0;
	uspech:=0;
	FOR i:=1 TO kroku DO BEGIN
		x:=rnd;
		y:=rnd;
		IF (sqr(x-0.5)+sqr(y-0.5))<sqr(0.5) THEN z:=z+f(x,y);
	END;
	integral:=z/kroku;

	WRITELN('Monte Carlo integral v kruhu druhou metodou ',integral:5:4);

	READLN;

END.