Численное интегрирование Опции

[trminator]

Численное интегрирование методом прямоугольников:

Код

program integral;
 {        Вычисляет приближенное значение       }
 { интеграла функции F методом прямоугольников  }
 { /b                                           }
 { |  f(x) dx;интервал разбивается на n частей  }
 { /a                                           }

var
 n,i : integer;
 a,b,shag,sum,itog : real;

{======================================}
{Введите сюда нужную функцию}
function F(x:real):real;
 begin
   F:=x*x*x*x+2*x*x+4      
 end;
{======================================}

BEGIN
 write('Начало интегрирования a = '); readln(a);
 write('Конец  интегрирования b = '); readln(b);
 write('Количество разбиений интервала n = '); readln(n);

 shag:=(b-a)/n;
 sum:=0;
 for i:=1 to n-1 do
   sum := sum + F(shag*i+a);
 sum := sum + (F(a)+F(b))/2;

 itog:=(b-a)/n * sum;
 writeLn('Интеграл = ', itog:0:5)
END.

[Altair]

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

UNIT C_Rect;

INTERFACE

TYPE Func=Function(x:real):real;

FUNCTION  Rectangles(a,b:Real; f:Func; n:word):Real;
FUNCTION RectanglesRunge(a,b:Real; f:Func; var n: word;eps:Real):Real;

IMPLEMENTATION

FUNCTION  Rectangles(a,b:Real; f:Func; n:word):Real;

VAR
  h,f1,sum,x,y:Real;
  i:Integer;

BEGIN

  h := (b-a) / n;
  sum := 0;
  x := a;
  for i := 1 to n  do
     begin
       y := x;
       x := x + h;
       f1 := f((x + y)/2);
       sum := sum + f1;
     end;
  Rectangles:= sum*h;
END;

FUNCTION RectanglesRunge(a,b:Real; f:Func; var n: word; eps:Real):Real;

VAR
  I1, I2: real;
BEGIN
   I1 := Rectangles(a,b,f,n);
   n := n+n;
   I2 := Rectangles(a,b,f,n);
   while abs(I1-I2) > 2*eps do
     begin
       if n >=  16383 then break;
       I1 := I2;
       n := n+n;
       I2 := Rectangles(a,b,f,n);
     end;
  RectanglesRunge := I2;
END;

BEGIN
END.

ОПИСАНИЕ МОДУЛЯ.  C_Rect.doc ( 59 килобайт ) Кол-во скачиваний: 3447