Разработка структур данных, необходимых для решения задачи

Переменные целого типа:

i — переменная цикла,.

n — количество интервалов разбиения Переменные вещественного типа:

a — нижняя граница интегрирования;

b — верхняя граница интегрирования;

R — параметр функции;

e — точность вычислений;

I, I2 — вычисленные значения интеграла;

x — переменная функции;

h — величина шага между двумя соседними точками разбиения интервала интегрирования.

Разработка структур алгоритма решения задачи

Ниже приведена схема алгоритм вычисления значение интеграла (рис. 3, рис. 4, рис. 5, рис.6).

Рисунок 3 — Ввод и проверка переменных.

Рисунок 4 — Ввод и проверка переменных (продолжение).

Рисунок 5 — Вычисление интеграла функции и вывод результата.

Текст программы на языке программирования Turbo Pascal 7.0

Далее текст программы «LinkinPark» на языке программирования IDE Borland Turbo Pascal 7.0, выполняющей вычисления таблицы значений функции, заданной в виде разложения в ряд.

program LinkinPark;

Uses crt;

var pr, j, n:integer;

a, b, R, e, I, I2, x, h, fi, fi2, pi:real;

str1:string;

begin.

clrscr;

I2:=0;

repeat.

writeln ('Введите нижнюю границу интегрирования. a=');

readln (str1);

val (str1,a, pr);

if pr0 then.

writeln ('Ошибка, введен символ. Введите цифру.');

until pr=0;

repeat.

writeln ('Введите верхнюю границу интегрирования b=');

readln (str1);

val (str1,b, pr);

if pr0 then.

writeln ('Ошибка. Введен символ. Введите цифру.');

if (a=b) or (a>b) then.

writeln ('Ошибка. аb и b>a').

if b=0 then.

writeln (`Ошибка. b0');

until ((pr=0)and (b>a)and (b0));

repeat.

writeln ('Введите параметр функции R');

readln (str1);

val (str1,R, pr);

if pr0 then.

writeln ('Ошибка. Введен символ. Введите цифру.');

until pr=0;

repeat.

writeln ('Введите точность вычислений');

readln (str1);

val (str1,e, pr);

if pr0 then.

writeln ('Ошибка. Введен символ. Введите цифру.');

until pr=0;

I:=0;

pi:=3.14 159 265;

n:=10;

while abs (I2-I)<=e do.

begin.

h:=(b-a)/n;

I2:=I;

I:=0;

for j:=1 to n do.

begin.

x:=a+h*(j-1);

if x<-R then.

fi:=-(x+R)*(sin (pi/3)/cos (pi/3)).

else.

begin.

if x.

fi:=cos (x).

else.

begin.

if x<3*R then.

fi:=sqrt (R*R-((x-2*R)*(x-2*R))).

else.

fi:=(x-3*R)*(sin (pi/6)/cos (pi/6));

end;

end;

I:=I+fi*h;

n:=n*2;

end;

end;

writeln (I:5:7);

end.