Вычисление интеграла методом прямоугольников

[mustard]

Никак не могу разобраться и понять, как написать программу в паскале, уже кучу раз писал и переписывал и ничего не выходит
Вот программа, которую написал:

var
a,b,x,y,h,S: real;
i,n: integer;
begin
readln (a,b,n);
h:=(b-a)/n;
S:=0;
for i:=1 to n do
begin
x:=a+i*h;
y:=((1-cos(2*x))/(1+cos(2*x)))/(1+sin(x)/cos(x));
S:=S+y;
x:=x+h;
end;
S:=h*S;
writeln (S);
end.

[type_Oleg]

Зачем в цикле

x:=a+i*h;
// ...
x:=x+h;  // это лишнее

И почему функцию Tan нельзя использовать? Если что, uses Math;

[type_Oleg]

Вот, проверял, работает

program Inte;
uses math;  // если надо
var
 a,b,x,y,h,S,t: real;
 i,n: integer;
begin
 readln (a,b,n);
 h:=(b-a)/n;
 S:=0;
 for i:=1 to n do
  begin
   x:=a+i*h;
   t:=Tan(x);  // лучше так, чтобы быстрее
   y:=(1+t*t)/(1+t);
   S:=S+y;
  end;
 S:=h*S;
 writeln (S);
end.

Число пи/4 правилно вводил? В смысле точки или запятой.

[mustard]

type_Oleg, спасибо огромное за помощь, всё сходится :3