Предыдущая страница реферата | 1  2  3  4 | Следующая страница реферата

Пологая для краткости y0+yn=Yкр (крайние), y1+y2+…+yn-1=Yпром (промежуточные), получим

	ò ydx » ((b-a)/2* n)*(Yкр+2*Yпром)

                                                                                                                

Эту формулу можно записать в другом виде

ò f(x)dx » (h/2)*[f(a)+f(b)+2åf(xi)]

(где h – длина одного из n равных отрезков, xi=a+i*h). Эта приближенная формула и  называется формулой трапеций. Она оказывается тем более точной, чем больше взятое нами число n. Погрешность одного шага вычисляется по формуле: -(h^3)/12.

        Задача. Пусть нужно проинтегрировать функцию f(x) = x³ +2x²-3x-8 на отрезке [0, 6]. На этом отрезке функция непрерывна.

Для выполнения поставленной задачи составлена нижеописанная программа,  приближенно вычисляющая определенный интеграл с помощью метода трапеций. Программа состоит из трех функций main, f и trap. Функция main позволяет ввести интервалы интегрирования и задать точность вычисления интеграла, а также вызывает функцию trap для вычисления интеграла и распечатывает на экране результат. Функция f принимает аргумент x типа float и возвращает значение интегрируемой функции в этой точке. Trap – основная функция программы: она выполняет все вычисления, связанные с нахождением определенного интеграла. Trap принимает четыре параметра: пределы интегрирования типа float (a и b), допустимую относительную ошибку типа float и указатель на интегрируемую функцию. Вычисления выполняются до тех пор, пока относительная ошибка, вычисляемая по формуле | S-Sn |,  не будет меньше или равна требуемой. Функция реализована с экономией вычислений, т. е. учитывается, что S0 постоянная и S1=S1+f(a+(2*i+1)*h), поэтому эти значения вычисляются единожды. Метод трапеций обладает высокой скоростью вычисления, но меньшей точностью, чем метод Симпсона, поэтому его применение удобно там, где не требуется очень высокая точность.

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

2. Блок-схема программы:

eps" >                                                                                                          

                                                                                                                      ДА

                                                                  НЕТ

                                                                       

                                                                                                                             i=1

S1=S1+f(a+(2*i+1)*h)

 

                                                                                                                         

                                                                                                                             i=n/2

                                                                         

3. Листинг:

Рекомендуем скачать другие рефераты по теме: атлетика реферат, сообщение об открытии счета.





Предыдущая страница реферата | 1  2  3  4 | Следующая страница реферата