Вычисление интеграла

program integral;
uses crt;
label 1,2;
var h,s,c,fa,fb:real;
i:integer;
v1,v2:char;
const a=1;
const b=4;
const n=60;
begin
1:clrscr;
writeln('Pramoyg ili trapezia');
writeln('Pramoyg#1');
writeln('Trapezia#2');
v1:=readkey;
case v1 of
'1':begin
clrscr;
h:=(b-a)/n;
writeln('h=',h:4:2);
s:=0;
c:=a-h/2;
for i:=1 to n do
begin
c:=c+h;
s:=s+(ln(c)*ln(c))/c;
writeln(i,'':4,'c=',c:4:2,'':2,'s=',s:4:2);
end;
s:=s*h;
writeln('i=',s:4:2);
readln;
end;
'2':begin
clrscr;
h:=(b-a)/n;
writeln('h=',h:4:2);
fa:=sqr(a);
fb:=sqr(b);
s:=(fa+fb)/2;
c:=a;
for i:=1 to n-1 do
begin
c:=c+h;
s:=s+(ln(c)*ln(c))/c;
writeln(i,'':4,'c=',c:4:2,'':2,'s=',s:4:2);
end;
s:=s*h;
writeln('i=',s:4:2);
readln;
end;
else writeln('Error');
end;
writeln('Выйти Y - да/N - нет?');
v2:=readkey;
case v2 of
'Y','y':goto 1;
'N','n':goto 2;
end;
2:end.

program integral;
uses crt;
const a=1;
const b=4;
const n=60;
function f(x:real):real;
begin
f:=sqr(ln(x))/x;
end;
procedure rect;
var h,x,s:real;
    i:integer;
begin
clrscr;
h:=(b-a)/n;
x:=a+h/2;
s:=0.0;
for i:=1 to n do
 begin
  s:=s+f(x);
  x:=x+h;
 end;
s:=s*h;
writeln('Интеграл методом средних прямоугольников=',s:0:4);
readln;
end;
procedure trap;
var h,s:real;
    i:integer;
begin
clrscr;
h:=(b-a)/n;
s:=(f(a)+f(b))/2;
for i:=1 to n-1 do
s:=s+f(a+i*h);
s:=s*h;
writeln('Интеграл методом трапеций=',s:0:4);
readln;
end;
var v1:char;
begin
clrscr;
repeat
clrscr;
writeln('Выберите действие');
writeln('1-интеграл методом прямоугольников');
writeln('2-интеграл методом трапеций');
writeln('другое выход');
readln(v1);
case v1 of
'1':rect;
'2':trap;
else exit;
end;
until not(v1 in ['1','2']);
readln
end.