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

@pavelIST · Регистрация: 05.01.2014

Студворк — интернет-сервис помощи студентам

Написал программу в делфи для вычисления определенного интеграла методом трапеции, но ответ в делфи не сходится сответом в мадкаде. Интеграл в строке(нижний предел - 0,7, верхний - 2,1)
(Sqrt(1.7*x*x+0.5)/(1.4+sqrt(1.2*x+1.3))

Код

Delphi
1
2
3
4
5
6
7
8
9
10
11
12
13
var
  a,b,h,S,x,y:Real;
begin
    a:=StrToFloat(edt1.Text);
    b:=StrToFloat(edt2.Text);
    h:=StrToFloat(edt3.Text);
    x:=0;
    S:=0;
    y:=h/2*((Sqrt(1.7*x*x+0.5)/(1.4+sqrt(1.2*x+1.3))))+(Sqrt(1.7*(x+h)*(x+h)+0.5)/(1.4+sqrt(1.2*(x+h)+1.3)));
    S:=S+y*h;
    Mmo1.Lines.Add('S= '+FloatToStr(S));
    x:=x+h;
end;

Нужно найти ошибку, или что-то изменить.
Игрик(y) высчитывал по формуле h/2(f(x)+f(x+h)) Скорее всего проблемма в S, но не найду ее. / Спасибо/

@Puporev · 07.01.2014, 08:03

Так интеграл как сумму площадей трапеций нужно считать в цикле, например по числу разбиений

Delphi
1
2
3
4
5
6
7
8
s:=0;
h := (b-a)/n;//n-число разбиений
for i:=1 to n-1 do
 begin
  x:= a + i* h;
  s:= s + f(x);
 end;
y:= h*((f(a)+f(b))/2+s);//f()-ваша функция

Новые блоги и статьи Все статьи Все блоги /
Модель заражения группы наркоманов alhaos 17.04.2026 Условия задачи сформулированы тут Суть: - Группа наркоманов из 10 человек. - Только один инфицирован ВИЧ. - Колются одной иглой. - Колются раз в день. - Колются последовательно через. . .	Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .	My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .	Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.
Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .	Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .	Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .