Нахождение интегралов в среде Pascal

Нахождение интегралов в среде Pascal

Текст задания: найти интеграл функции sin (x+10)+x⁴=0 используя два метода: метод прямоугольников и метод трапеций; пользователь вводит значения и выбирает метод решения.

Текст программы:

program LW_5;

uses

Crt;

var

dx, x1, x2, e, i, c, s:real;

function Fx (x:real):real;

begin

Fx:=sin (x+10)+x*x*x*x; {function}

end;

procedure CountViaBar;

var

xx1, xx2:real;

c:longint;

begin

writeln ("–");

writeln ("–>MEthod PRYAMOUGOL""NIKOVVVV");

writeln ("vsego interaciy:", round (abs(x2-x1)/e));

i:=0;

for c:=1 to round (abs(x2-x1)/e) do begin

write ("integration ", c, chr(13));

xx1:=Fx (x1+c*e);

xx2:=Fx (x1+c*e+e);

i:=i+abs (xx1+xx2)/2*e;

end;

writeln ("–");

writeln ("INTEGRAL=", i);

end;

procedure CountViaTrap;

var

xx1, xx2, xx3:real;

c:longint;

begin

writeln ("–");

writeln ("–>MEthod TRAPECIY");

writeln ("vsego interaciy:", round (abs(x2-x1)/e));

i:=0;

for c:=1 to round (abs(x2-x1)/e) do begin

write ("Integration ", c, chr(13));

xx1:=Fx (x1+c*e);

xx2:=Fx (x1+c*e+e);

if xx2>xx1 then xx3:=xx1 else xx3:=xx2;

i:=i+abs (xx2-xx1)*e+abs(xx3)*e;

end;

writeln ("–");

writeln ("INTEGRAL=", i);

end;

begin

clrscr;

writeln ("vvedite ishodniye zna4eniya:");

write ("NA4AL""NOYE ZNE4ENIYE x (x1)="); Readln(x1);

write ("KONE4NOYE ZNA4ENIYE x (x2)="); Readln(x2);

write ("TO4NOST"" VbI4ISLENIY (e)="); Readln(e);

writeln ("–======vvedite metod====== – :");

writeln ("DLYA METODA PRYAMOUGOL""NIKOV VVEDITE 1, DLYA METODA TRAPECIY VVEDITE 2");

 //writeln ("DLYA METODA TRAPECIY VVEDITE 2");

readln (s);

if s=1 then begin

CountViaBar;

end;

if s=2 then begin

CountViaTrap;

end;

readln;

end.