VPF::[VBA] Оценка остаточного члена в формуле Симпсона

STRIKER221

Дата 14.3.2010, 18:14 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 84
Регистрация: 2.11.2006
Где: Belarus

Репутация: нет
Всего: нет

Необходимо вычислить интеграл по формуле Симпсона определяя шаг интегрирования по оценке остаточного члена.
Метод Симпсона я реализовал:

Код


Private Sub SIMP() 'Метод Симпсона
    n = (b - a) / h
    
    S1 = 0
    For i = 1 To 2 * n - 1 Step 2
        S1 = S1 + f(a + (b - a) * i / (2 * n))
    Next i
    
    S2 = 0
    For i = 2 To 2 * n - 2 Step 2
        S2 = S2 + f(a + (b - a) * i / (2 * n))
    Next i
    
    s = ((b - a) / (6 * n)) * ((f(a) + f(b) + 4 * S1 + 2 * S2))
    ListBox1.AddItem "По формулам Симпсона"
    ListBox1.AddItem "Интеграл = " & Round(s, Len(eps))
    ListBox1.AddItem "--------------------------------------"
End Sub

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

Это сообщение отредактировал(а) STRIKER221 - 15.3.2010, 00:22

Sanaff

Дата 14.3.2010, 22:00 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 151
Регистрация: 15.9.2009
Где: г. Северодвинск

Репутация: -1
Всего: 1

Кажется в Википедии есть. Там про Симпсона много написано.

--------------------

Программист - это локальный бог ©ICQ 373-628-456

STRIKER221

Дата 15.3.2010, 00:23 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 84
Регистрация: 2.11.2006
Где: Belarus

Репутация: нет
Всего: нет

Смотрел, про оценку остаточного члена ничего не нашел

STRIKER221

Дата 16.3.2010, 04:34 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 84
Регистрация: 2.11.2006
Где: Belarus

Репутация: нет
Всего: нет

Сам разобрался. Вот код:

Код


Private Sub SIMP() 'Метод Симпсона
    Dim i As Single
    Dim j As Single
    
    If q = 1 Then
        Exit Sub
    End If

    M = 0
    
    For j = a To b Step 0.1
        Z = Abs(df_4(j))
        If Z > M Then
            M = Z
        End If
    Next j
        
    h = (eps * 180) / ((b - a) * M)
    n = Fix((b - a) / (2 * h))
    
    s1 = 0
    For i = 1 To 2 * n - 1 Step 2
        s1 = s1 + f(a + (b - a) * i / (2 * n))
    Next i
    
    S2 = 0
    For i = 2 To 2 * n - 2 Step 2
        S2 = S2 + f(a + (b - a) * i / (2 * n))
    Next i
    
    s = ((b - a) / (6 * n)) * ((f(a) + f(b) + 4 * s1 + 2 * S2))
    ListBox1.AddItem "По формуле Симпсона"
    ListBox1.AddItem "Интеграл = " & Round(s, Len(eps))
    ListBox1.AddItem "--------------------------------------"
End Sub

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »