Паскаль. Вычисление интегралов.

[jukers]

Составить программу, подсчитывающую с точностью Е=10^-3 значение интеграла : 1)интеграл от Пи/6 do Pi/3(tg^2(x) + ctg^2(x))*d*x) по методу Правых прямоугольников. Для точности использовать метод двойного пересчета. Результат срав-нить со значением интеграла, подсчитанным с помощью первообразной F(x) по формуле Ньютона-Лейбница: F(b)-F(a).
Метод приближенного вычисления определенного интеграла должен быть запро-граммирован с помощью процедуры. Отрезок интегрирования, точность вычисления вводятся с клавиатуры, подынтегральная функция описывается с помощью функции пользователя.
Распечатать: отрезок интегрирования, точность вычисления, приближенное значение определенного интеграла, значение интеграла, подсчитанное по формуле Ньютона-Лейбница.

[newerow1989]

function f(x:real):real;
begin
   Result:=Tan(2*x)+CoTan(2*x);
end;

var a,b,dx,s,x:real;
begin
   Write('a = '); Readln(a);
   Write('b = '); Readln(b);
   Write('dx = '); Readln(dx);
   s:=0;
   x:=a-dx;
   Repeat
      x:=x+dx;
      s:=s+f(x);
   Until x>=b;
   Writeln(s:5:5);
   Readln;
end.