Организовать чтение из файла

program tr_s;
uses crt;
var
a,b:real; { Границы отрезка }
r,r2:real; { Предыдущее и текущее приближенные значения интеграла}
n:integer; { Счетчик }
{ Интегрируемая функция }
function f(x:real):real;
begin
f:=1/(x*ln(x)*0.43429);
end;
{ Метод трапеций }
function trap(a,b:real;n:integer):real;
var
s:real; { Полученная сумма }
h:real; { Шаг }
m:integer; { Счетчик }
begin
h:=(b-a)/(n-1); { Определяется шаг }
s:=(f(a)+f(b))/2; { Начальное значение суммы }
for m:=1 to n-2 do s:=s+f(a+m*h); { Суммиование остальных элементов}
trap:=s*h; { Возвращается значение интеграла }
end;
{ Метод Симпсона }
function simpson(a,b:real;n:integer):real;
var
s:real; { Сумма }
h:real; { Шаг }
m:integer; { Счетчик }
mn:integer; { Очередной множитель }
begin
h:=(b-a)/(n-1); { Рассчитывается шаг }
s:=f(a)+f(b); { Начальное значение шага }
mn:=4; { Первый мнодитель - 4 }
{ Суммирование остальных элементов }
for m:=1 to n-2 do begin
s:=s+mn*f(a+h*m);
if (mn=4) then mn:=2 else mn:=4;{ Именение мноителя 2<>4 }
end;
simpson:=s*h/3; { Возвращается вычисленное значение }
end;
{ Основная программа }
begin
{ Ввод границ отрезков }
clrscr;
write(' Введите A,B: ');
readln(a,b);
{ Вычисляется интеграл по методу трапеций }
n:=3;
r:=trap(a,b,n); { Начальное значение }
repeat
r2:=r; { Запоминается предыдущее значение }
n:=n+2; { Увеличивается количество шагов }
r:=trap(a,b,n); { Рассчитывается новое значение }
until (abs(r-r2)<0.001);{ Повторяется до достижения необходимой точности }
{ Вывод результатов }
writeln(' Резльтат по методу трапеций равен: ',r:6:3);
writeln(' для получения необходимой точности интервал был разбит на');
writeln(n,' отрезков');
{ Вычисляется интеграл по методу Симпсона }
n:=3;
r:=simpson(a,b,n); { Начальное значение }
repeat
r2:=r; { Запоминается предыдущее значение }
n:=n+2; { Увеличивается количество шагов }
r:=simpson(a,b,n); { Рассчитывается новое значение }
until (abs(r-r2)<0.001);{ Повторяется до достижения необходимой
точности }
{ Вывод результатов }
writeln;
writeln(' Резльтат по методу Симпсона равен: ',r:6:3);
writeln(' для получения необходимой точности интервал был разбит на ');
writeln(n,' отрезков');
end.

var
a,b : integer;
f: text;
begin
assign (f,' {адрес файла}' ) ;
reset (f);
read (f, a,b);
close(f);
end.