Наближене обчислення визначених інтегралів, що не беруться через елементарні функції - Курсова робота

.

Якщо проміжок розділити на рівних частин, то – для формули Сімпсона (10) – отримаємо залишковий член у вигляді

(16).

При зростанні цей вираз зменьшується приблизно як ; таким чином, формула Симпсона дійсно більш вигідна, ніж попередні дві формули.

Додаток 1.

Текст программи для автоматичного обчислення інтегралів на мові програмування QBASIC:

'Тут описуються сталі

e = 2.718281828459045#

pi = 3.141592653589793#

'Тут задається від під інтегральної функції

DEF fny# (x#) = e^x# ^2

DEF fncoef# (i#) = (i# MOD 2) * 2 + 2

DEF fnxi# (i#) = a# + i# * h#

DEF fnxis# (i#) = a# + i# * h# / 2

DEF fnxic# (i#) = a# + i# * h# + h# / 2

DEF fnxir# (i#) = a# + i# * h# + h# / 2

CLS

'Тут вводяться межі інтегрування та

'кількість проміжків

INPUT «Введіть нижню межу інтегрування » a#

INPUT «Введіть верхню межу інтегрування » b#

INPUT «Введіть кількість проміжків » n#

'Тут обчислюється крок

h# = (b# - a#) / n#

'Тут обчислюється наближене значення

'інтеграла за методом Сімпсона

integ# = 0

FOR i# = 1 TO ((2 * n#) - 1)

integ# = integ# + fncoef#(i#) * fny#(fnxis#(i#))

NEXT

integ# = integ# + fny#(a#) + fny#(b#)

integ# = integ# * (h# / 6)

PRINT "Simpson = "; integ#

'Тут обчислюється наближене значення

'інтеграла за методом трапецій

integ# = 0

FOR i# = 1 TO (n# - 1)

integ# = integ# + fny#(fnxi#(i#))

NEXT

integ# = integ# + (fny#(a#) + fny#(b#)) / 2

integ# = integ# * h#

PRINT "Trapeze = "; integ#

'Тут обчислюється наближене значення

'інтеграла за методом лівих прямокутників

integ# = 0

FOR i# = 0 TO (n# - 1)

integ# = integ# + fny#(fnxi#(i#))

NEXT

integ# = integ# * h#

PRINT "L Rectangle = "; integ#

'Тут обчислюється наближене значення

'інтеграла за методом центральних прямокутників

integ# = 0

FOR i# = 0 TO n#

integ# = integ# + fny#(fnxic#(i#))

NEXT

integ# = integ# * h#

PRINT "C Rectangle = "; integ#

'Тут обчислюється наближене значення

'інтеграла за методом правих прямокутників

integ# = 0

FOR i# = 1 TO n#

integ# = integ# + fny#(fnxir#(i#))

NEXT

integ# = integ# * h#

PRINT "R Rectangle = "; integ#

Додаток 2.

Далі подані результати роботи програми, яка викладена в додатку 1.

1) в межах від 0 до

n=1000

Метод Сімпсона -8.742278155181581D-08

Метод трапецій -8.742270585611512D-08

Метод лівих прямокутників 3.141505318306509D-03

Метод центральних прямокутників -3.14167628761223D-03

Метод правих прямокутників -6.283265152840917D-03

2) в межах від 0 до

n=1000

Метод Сімпсона 2.000000000000067

Метод трапецій 1.999998355065565

Метод лівих прямокутників 1.999998355202888

Метод центральних прямокутників 1.999995887392223

Метод правих прямокутників 1.999990952591778

3) в межах від 0 до 1

4) в межах від 0 до 1

n=1000

Метод Сімпсона .7468241385662959

Метод трапецій .7468240772530558

Метод лівих прямокутників .7471401375268841

Метод центральних прямокутників .7471916808878213

Метод правих прямокутників .7461916811378212

5) в межах від 0 до

n=1000

Метод Сімпсона .8323745796964475

Метод трапецій .8323723082182791

Метод лівих прямокутників .8325874590746988

Метод центральних прямокутників .8319367429487694

Метод правих прямокутників .8319318081462942

Висновки.

У данній роботі було розглянуто методи наближених обчислень визначених інтегралів, були виведині формули обчислень, формули додаткових членів. Результати, які наведені в додатку 2 наочно показують, що найбільш вигідним є використання формули Сімпсона.

Література.