Нахождение интегралов в среде 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.