close

Вход

Забыли?

вход по аккаунту

?

Integrirovanie

код для вставкиСкачать
Московский Технический Университет Связи и Информатики
Кафедра Информатики
Лабораторная работа №6.3 по дисциплине информатика (спецглавы)
за 3 семестр
"Численное интегрирование"
Студента 2 курса
Митрофанова Егора
Группа БИН1207
Преподаватель: Гловацкая А. П.
Задание:
1. Выбрать индивидуальное задание для численного интегрирования: • f(x) - подынтегральную функцию; • a, b- пределы интегрирования; • метод интегрирования на компьютере -по формуле Симпсона
• метод интегрирования вручную - по формуле трапеций
• начальный шаг интегрирования h0. 2. Составить схему алгоритма и написать программу по выбранному методу численного интегрирования (или по указанному преподавателем), провести контрольное тестирование 3. Вычислить интеграл I=∫_a^b▒f(x)dx с точностью E= 0.1*〖10〗^(1-i), i=1, 2, 3, 4... 4. Вычислить "ручным расчетом" интеграл методом, определяемым значением столбца m из таблицы 1.4-1, с шагом h и h/2 (I_h0 и I_(h0/2)) и оценить погрешность по правилу Рунге. Вариант задания:
18cos(2x + 1) ln (2 / x) + 3
a=1, b=3
На компьютере: Интегрирование по формуле Симпсона
Вручную: По формуле трапеций
Шаг интегрирования:
h=0.5
Схема алгоритма по формуле Симпсона
Результаты вычисления "расчета на ПК"
0.140.55,5664860.0180.255,5627080.00180.255,562708
Зависимости числа итераций от точности
E0.10.010.001n488
Код программы:
Imports System.Math
Public Class Form1
Function f(ByRef x As Single) As Single
Return Cos(2 * x + 1) * Log(2 / x) + 3
End Function
Sub resh(ByVal a As Single, ByVal b As Single, ByVal E As Single, ByRef s As Single, ByRef n As Integer, ByRef h As Single)
Dim s1 As Single
Dim c, x As Single
n = 2
h = (b - a) / n
s = (f(a) + 4 * f(a + b) / 2 + f(b)) * h / 3
Do
n = 2 * n
h = (b - a) / n
s1 = s
c = 4
x = a
s = f(a) + f(b)
For i = 1 To n - 1
x = x + h
s = s + c * f(x)
c = 6 - c
Next
s = s * h / 3
Loop Until Abs(s - s1) / 15 < E
End Sub
Private Sub Button1_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles Button1.Click
Dim a, b, s, h As Single
Dim n As Integer
Dim Viv As String
a = TextBox1.Text
b = TextBox2.Text
TextBox3.Text = " "
s = 0
h = 0
resh(a, b, 0.1, s, n, h)
Viv = Space(2) + CStr(s) + Space(20) + CStr(n) + Space(30) + CStr(h) + vbNewLine
resh(a, b, 0.01, s, n, h)
Viv = Viv + Space(2) + CStr(s) + Space(20) + CStr(n) + Space(30) + CStr(h) + vbNewLine
resh(a, b, 0.001, s, n, h)
Viv = Viv + Space(2) + CStr(s) + Space(20) + CStr(n) + Space(30) + CStr(h) + vbNewLine
TextBox3.Text = Viv
End Sub
End Class
1
Документ
Категория
Рефераты
Просмотров
7
Размер файла
148 Кб
Теги
integrirovanny
1/--страниц
Пожаловаться на содержимое документа