Проблемы в Pascal (для всех компиляторов!) Опции

[ONi++]

Программу одну составил по олимпиадным задачам. Суть задачи состоит в том, что нужно узнать по данным координатам, находится ли точка в треугольнике. Я сравниваю площадь треугольника и площадь четырехугольника (получившегося из соединения линиями треугольника с точкой) с помощью функции round(). Но, когда она [функция] округляет 4.5 два раза, то сначала получается 5, а потом 4. Что делать? Не хотелось бы менять весь алгоритм из-за какой-то ошибки в компиляторах (кстати, я проверял и на PascalABC, и на TurboPascal, и на FreePascal).

[ONi++]

var
input,output:text;
x,y:array[1..4]of integer;
i,n,xx,yy:integer;
a,b,c,p,s1,s2:real;
begin
  assign(input,'tria-pt.in');
  assign(output,'tria-pt.out');
  reset(input);
  rewrite(output);
  n:=3;
  i:=1;
  while i<=n do
	 begin
		readln(input,x[i],y[i]);
		i:=i+1;
	 end;
  readln(input,xx,yy);
  a:=sqrt(sqr(x[1]-x[2])+sqr(y[1]-y[2]));
  b:=sqrt(sqr(x[2]-x[3])+sqr(y[2]-y[3]));
  c:=sqrt(sqr(x[3]-x[1])+sqr(y[3]-y[1]));
  p:=(a+b+c)/2;
  s1:=sqrt(p*(p-a)*(p-b)*(p-c));
  x[n+1]:=x[1];
  y[n+1]:=y[1];
  i:=1;
  while i<=n do
	 begin
		a:=sqrt(sqr(xx-x[i])+sqr(yy-y[i]));
		b:=sqrt(sqr(x[i]-x[i+1])+sqr(y[i]-y[i+1]));
		c:=sqrt(sqr(x[i+1]-xx)+sqr(y[i+1]-yy));
		p:=(a+b+c)/2;
		s2:=s2+sqrt(p*(p-a)*(p-b)*(p-c));
		i:=i+1;
	 end;
  if round(s1)=round(s2) then
	 writeln(output,'In')
  else
	 writeln(output,'Out');
  close(input);
  close(output);
end.