Вычисление площади контура и координаты центра его веса

[lilywka]

Заданный замкнутый контур нарисовать на миллиметровке, вычислить его площадь и координаты центра тяжести с помощью простых геометрических формул. Создать программу для вычисления с помощью модифицированной формулы Грина, сравнить результаты (они должны совпадать).
10. Заданная функция F (x) = x - sin x на интервале (-2π ÷ +2 π), определить площадь и координаты центра тяжести контура между функцией и осью Х с помощью модифицированной формулы Грина. Сравнить с ручным расчетом (275, с. 228).

Пыталась сделать но не получилось. Что у меня не правильно?

program Project1;

{$APPTYPE CONSOLE}

uses
  SysUtils,
  Windows,
  math;

function f(x:real):real;

begin
  f:=x-sin(x);
end;
const
  Pi=3.14;
  step: real=0.1;
  Xmin:real=-2*Pi;
  Xmax:real=2*Pi;
var
  steps,i:integer;
  Suma,x,x1,y,y1,XC,YC:real;
begin
  SetConsoleCP(1251);
  SetConsoleOutputCP(1251);
  steps:=round((Xmax-Xmin)/step);
  Suma:=0;
  for i:=0 to steps -1 do 
    begin
      x:= Xmin+step*i;
      x1:=Xmin+step*(i+1);
      writeln('X   ','X1   ','Y   ','Y1   ');
      if(i<>60) then
        begin
          Y:=f(x);   
          Y1:=f(x1);
          Suma:=Suma+0.5*(x*y1-y*x1);
          XC:=XC+(x+x1)/3*0.5*(x*y1-y*x1);
          YC:=YC+(y+y1)/3*0.5*(x*y1-y*x1);
          writeln(x:4:2, ' ', x1:4:2, ' ', y:4:2, ' ',y1:4:2);
        end;
    end;
end;
writeln('площадь=', Abs(Suma):6:2);
writeln('центр мас, X: ', (XC/Suma):6:2);
write('            Y:',(YC/Suma):6:2);
readln;
  { TODO -oUser -cConsole Main : Insert code here }
end.