Форум программистов [Powered by Invision Power Board]

Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате

Форум программистов > Центр помощи > [VBA] Оценка остаточного члена в формуле Симпсона

Автор: STRIKER221 14.3.2010, 18:14

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

Код


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

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

Автор: Sanaff 14.3.2010, 22:00

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

Автор: STRIKER221 15.3.2010, 00:23

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

Автор: STRIKER221 16.3.2010, 04:34

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

Код


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