КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Формула Симпсона
Введите число точек - 100 Искомый определенный интеграл равен 1.999836
Введите левую и правую границы отрезка - 0 3.1415 Введите число точек - 1000 Искомый определенный интеграл равен 1.999998
Введите левую и правую границы отрезка - 0 3.1415 Введите число точек - 10000 Искомый определенный интеграл равен 2.000000
Найдем коэффициенты формулы , где , i=0,1,…,n при n= 2. При i = 0 При i = 1 При i = 2 Формула на отрезке [ x0, x2 ] примет вид: При n=2m применив формулу к каждой паре частичных отрезков [ x2i-2, x2i ] (i=1,2,…,m) получим формулу Симпсона: Рассмотрим погрешность: На отрезке [ a, b ] или , где При вычислении по методу повторного счета можно использовать формулу: . Если при вычислении интеграла требуемая точность не достигнута (т.е. ), предусматривается повторный счет с шагом, уменьшенным вдвое. Программа вычисления по формуле Симпсона методом повторного счета: program lab4_2; var n: integer; S,a,b,e,h,x,I_n,I_n2,M: real; function f(x: real):real; begin {записать, функцию в виде f:=[математическое выражение]} f:=sin(x); end; begin write('Введите концы отрезка интегрирования: '); readln(a,b); write('Введите погрешность e: '); readln(e); I_n:=0; n:=4; Repeat h:=abs(b-a)/n; s:=(f(a)-f(b))/2; x:=a+h; repeat s:=s+2*f(x)+f(x+h); x:=x+2*h; until x>=b; I_n2:=2*h*s/3; n:=n*2; M:=abs(I_n - I_n2); I_n:=I_n2; Until M<=15*e; writeln('Интеграл I=',I_n2:12:7); readln; end.
Дата добавления: 2014-01-06; Просмотров: 242; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |