Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Практикум на ЭВМ.doc
Скачиваний:
111
Добавлен:
15.02.2016
Размер:
2.73 Mб
Скачать

Лабораторная работа №5 Вычисление определенных интегралов и табулирование первообразных функций

Часто в научно-технических задачах возникает необходимость вычисления определенного интеграла или значений первообразной функции. Умея вычислить первообразную функцию

мы можем вычислить определенный интеграл

и наоборот. Но, как правило, выразить первообразную функцию через элементарные функции не удается. Поэтому приходится прибегать к приближенному интегрированию. Для решения этой задачи существует много численных методов, из которых мы рассмотрим два: метод трапеций и метод Сипсона.

Метод трапеций.

Разобьем отрезок интегрирования [a,b] на n равных частей длиной h=(b-a)/n. В точках разбиения x0=a, x1=a+h,... xi=a+ih,..., xn=b. Вычислим ординаты y0=f(x0), y1=f(x1),..., yi=f(x1).

Тогда приближенные значения интеграла методом трапеций вычисляется по формуле

Метод Сипсона.

Разобьем отрезок интегрирования [a,b] на 2n равных частей длиной h=b-a/2n. Пусть точкам разбиения x0=a, x1=a+h, x2=a+2h,..., xi=a+ih,...,

x2n-1=a+(2n-1)h, x2n=b. соответствуют значения подынтегральной функции

y0=f(x0)=f(a), y1=f(x1),..., yi=f(xi),..., y2n-1=f(x2n-1), y2n=f(x2n)=f(b).

Тогда формула Cипсона имеет вид

Требования к работе:

1. Записать расчетные формулы для решения задачи.

2. Составить блок-схему алгоритма и программу для вычисления определенного интеграла указанным в варианте методом (табл.2), разбивая отрезок интегрирования [a,в] на указанное число (n) частей. Предусмотреть в программе вычисление точного значения определенного интеграла через первообразную.

3. Составить блок - схему и программу вычисления значений первообразной. Предусмотреть в программе печать точных значений первообразной с тем же шагом и вычерчиванием графиков точного и приближенного значений первообразной.

4. Отладить обе программы.

Лабораторные задания

Таблица 2

вар.

Подынтегральная

функция

Промежуток интегр.

[a;b]

Метод численного решения определ. интегр.

Кол-во частей разб.

Шаг. вычисл. знач. перво-обр.

h

Точность вычисл. знач. первообр

Точное значение первообразной

1

[1; 3,5]

Симпсона

30

0,25

0,001

2

lg2x+ctg2x

трапеций

54

/36

0,001

3

[2; 3]

Симпсона

36

0,2

0,001

2,3026 (ln ln x - ln ln 2)

4

[1; 4]

трапеций

52

0,5

0,001

1/3 ln3x

5

[ 0; ln2]

Симпсона

104

0,001

6

xex sin x

[0; 1]

трапеций

48

0,2

0,001

7

x sh x

[0; 2]

Симпсона

48

0,4

0,001

8

[0; 2]

трапеций

208

0,25

0,001

9

[1; 2,5]

Симпсона

44

0,3

0,001

10

x arctg x

[0, ]

трапеций

48

/ 8

0,001

11

[0; 3]

Симпсона

36

0,5

0,001

12

xx(1+ln x)

[1; 3]

трапеций

40

0,2

0,001

xx-1

13

[0; 1]

Симпсона

44

0,2

0,001

14

1, 2]

трапеций

160

0,001

15

23x

[0; 1]

Симпсона

240

0,2

0,001

16

[0; 1]

трапеций

22

1/8

0,001

17

[0; 2]

Симпсона

48

0,25

0,001

18

sin2 x

трапеций

22

0,001

19

[0; 1,9999]

Симпсона

96

0,25

0,001

20

ex cos2 x

[0; ]

трапеций

60

/6

0,001

21

(x ln x)2

[1; e]

Симпсона

52

0,001

22

[0; 3]

трапеций

176

0,6

0,001

23

[0; 1]

Симпсона

36

0,25

0,001

24

sinx ln(tgx)

[1; 1,5]

трапеций

52

0,1

0,001

25

[0; 1,5]

Симпсона

132

0,3

0,001

26

[0;3/4]

трапеций

40

3/20

0,001

27

[0; 1]

Симпсона

78

0,2

0,001

28

[1; 2]

трапеций

40

0,2

0,001

29

[1; 2]

Симпсона

72

0,125

0,001

30

[1; 2]

трапеций

36

0,25

0,001