|
|
ru.algorithms- RU.ALGORITHMS ---------------------------------------------------------------- From : Evgenij Masherov 2:5020/175.2 14 Feb 2003 10:49:50 To : Alexander Shevchenko Subject : интегpал -------------------------------------------------------------------------------- Thu Feb 13 2003 16:00, Alexander Shevchenko wrote to Evgenij Masherov: AS>>> Hе да. То что помню, то и сделал :) EM>>>> Метод трапеций несколько точнее, Симпсона еше лучше. AS>>> Им и сделал. Спасибо. EM>>>> В методе трапеций крайние две точки берутся с половинным весом, EM>>>> в методе Симпсона EM>>>> S(f)=1/6*h*SUMi(f(x[i])+4*f((x[i]+x[i+1])/2)+f(x[i+1]) AS>>> Хм... у меня, мякого говоря, другой алгоритм... AS>>> h:=b; AS>>> s:=f(a)+f(b); AS>>> s2:=s; AS>>> repeat AS>>> s3:=s2; AS>>> h:=h/2; AS>>> s1:=a; AS>>> x:=h; AS>>> repeat AS>>> s1:=s1+2*(f(x)); AS>>> x:=x+2*h; AS>>> until (not(x<t)); AS>>> s:=s+s1; AS>>> s2:=(s+s1)*h/3; AS>>> x:=abs(s3-s2)/15; AS>>> until (not(x > ep)); AS>>> *Result*:=s2; EM>> Тут несколько по разному понимается h. AS> А как у тебя? Как правильно? И так, и так правильно. Просто приведенное мной выражение описывает расчет для отрезка длиной h. А в программе делаются шаги по 2*h. Евгений Машеров АКА СанитарЖеня --- ifmail v.2.15dev5 * Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2) Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
Архивное /ru.algorithms/3300932463bf.html, оценка из 5, голосов 10
|