Visual Basic FAQ

Вычисление определенного интеграла методом Симпсона

Private Sub Command1_Click()

A = 0 'нижний предел интегирования
B = 1 'верхний предел
h = 0.01 'шаг ( в твоей задаче это Е)
N = (B - A) / h 'кол-во точек, на которые разбивается диапазон

'вычисление первой суммы (см. формулу)
S1 = 0
For i = 1 To N - 1
  S1 = S1 + f(h * i + A)
Next i

'Вычисление второй суммы (см. формулу). Я преобразовал немного - 'так короче.
S2 = 0
For i = 1 To N
  S2 = S2 + f(h * i - 0.5*h  + A)
Next i

'Итоговая формула
S = h / 3 * (0.5 * f(A) + S1 + 2 * S2 + 0.5 * f(B))

MsgBox "Интеграл = " & S

End Sub
Function f(x) As Single
'Сама функция
  f = x * x + 1
End Function

Запись опубликована on 21.02.2010 at 15:42 и размещена в рубрике: Арифметика. Вы можете читать комментарии, используя RSS-ленту. Вы можете оставить комментарий, или отправить трекбек с Вашего сайта.