Методом трапеций высислить значение интеграла

[zloygeniyrus]

На заданном отрезке интегрирования[a,b].Cчитать заданным численный метод решения.
Включить в программу вычисление точного значения интеграла.На печать вывести количество отрезков разбиения интервала,приближенные значения интеграла,соответствующие количеству отрезков разбиения,точное значение интеграла и относительную погрешность вычисления.
Подынтегральная ф-ция :sin квадрат X;
Интервал интегрирования [0,pi/2];
Требуемая точность: 0.5*10 в -3 степени.

program IntTrap;
const
     e=0.5e-3;
     a=0;
     b=pi/2;
function F(x:real):real;
begin
  F:=sin(x)*sin(x);
end;

function pervF(x:real):real;
begin
  pervF:=1;{ya neznayu}
end;

var
  x,dx,s,trg,de:real;
  i,n:integer;
begin
  n:=1;
  repeat
    s:=0;
    n:=n+1;
    x:=a;
    dx:=(b-a)/n;
    de:=0;
    for i:=0 to n-1 do
      begin
        trg:=(F(x+dx)-F(x))*dx/2;
        s:=s+F(x)*dx+trg;
        de:=de+abs(trg);
        x:=x+dx;
      end;
  until de<e;

  writeln('kol-vo otrezkov      :',n);
  writeln('dostignutaya to4nost :',de:0:8);
  writeln('pribliz integral     :',s:0:8);
  writeln('tochniy integral     :',pervF(b)-pervF(a):0:8);
end.

Помогите плз. дописать первообразную и вообще дописать текст проги)
Заранее спасибо