Вычислить определённый интеграл заданным методом с заданной точностью.

[461R]

Помогите исправить программу, препод сказал что константа используется как шаг, а должна использоваться как точность.
Метод прямоугольников по серединам.

program Integral;
Uses Crt;
var
x1,x2,i:real;
xx1,xx2:real;
c:longint;
const
e=0.001;
function Fx(x:real):real;
begin
Fx:=ln(x)/(x*sqrt(1+ln(x)));
end;
begin
write('Введите значение a=');Readln(x1);
write('Введите значение b=');Readln(x2);
begin
for c:=1 to round(abs(x2-x1)/e) do begin
xx1:=Fx(x1+c*e);
xx2:=Fx(x1+c*e+e);
i:=i+abs(xx1+xx2)/2*e;
end;
writeln('Интеграл=',i:0:3);
end;
end.

[ViktorR]

Определи некотурую переменную n - целого типа. Эта переменная - число разбиений отрезка, на котором вычисляется интеграл. Пусть для начала

n := 10;

Получи значение интеграла. Измени n, например,

n := n*2;

Получи новое значение интеграла и разность с первым значением:

if abs(i1-i2) <= eps then ...

Получено значение интеграла с заданной точностью.



Как-то так, ...

[Smitt&Wesson]

Народ! Ну, я уже ржать устал! Имейте совесть!

[461R]

Не могу понять что куда подставлять в программе, можешь скинуть готовую программу?

[461R]

VictorR, не могу понять что куда подставлять, можешь скинуть готовый вариант программы?

[ViktorR]

Ну вот, ...

program Integral;
Uses Crt;
type fun = f(x:real) : real; {определим функциональный тип}
{$F+}
function F(x:real):real; {опишем нашу функцию}
begin
   F:=ln(x)/(x*sqrt(1+ln(x)));
end;
{$F-}

function Integr(f : fun; a, b : real; n : integer) : real; {Вычисление интеграла}
var sum, x, h : real;
   i : integer;
begin
   h := (b - a) / n;
   sum := 0;
   x := a + h/2;
   for i := 1 to n do
   begin
     sum := sum + f(x)
     x := x + h;
   end;
   integr := sum * h;
end;

var x1, x2, eps :real;
   si, sn : real;
      n : integer;
begin
   write('Введите левую границу интегрирования a= ');
   Readln(x1);
   write('Введите правую границу интегрирования b= ');
   Readln(x2);
   write('Введите требуемую точность eps= ');
   Readln(eps);
   write('Введите начальное число разбиений отрезка интегрирования n= ');
   Readln(n);
   sn := 0; {инициируем первое приближение для интеграла}
   si := integr(f, x1, x2, n); {Второе приближение}
while abs(si - sn) > eps do {вычисляем более точно}
begin
   sn := si;
   n := 2*n;
   si := integr(f, x1, x2, n);
   if n > 500 then
     break; {Прекратим вычисленя, если число разбиений велико}
end;
writeln('Интеграл=', si:8:3);
end.

Можно выполнить эту работу и для другой функции, которую следует описать между символами {$F+} и {$F+}

Как-то так, ...

PS: Так стало легче?

[461R]

намного) спасибо огромное)