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

            int n = 4, i; /* Начальное число разбиений 4 */

            float s_ab = f(a) + f(b); /* Сумма значений фун – ции в a и b */

float h = (b – a) / n; /* Вычисляем шаг */

            float s_even = 0,  s_odd;

            float s_res = 0, s_pres;

            /* Сумма значений фун – ции в нечетных точках */

            for (i = 2; i < n; i += 2) {

                       s_even += f(a + i * h);

}

            do {

                       s_odd = 0;

                       s_pres = s_res;

                       

/* Сумма значений фун – ции в четных точках */

            for (i = 1; i < n; i += 2) {

                                   s_odd += f(a + i * h);

}

            /* Подсчет результата */

                       s_res = h / 3 * (s_ab + 2 * s_even + 4 * s_odd);

/* Избегаем деления на ноль */

                       if (s_res == 0) s_res = e;

                       s_even += s_odd;

                       n *= 2;

                       h /= 2;

} while (fabs((s_pres - s_res) / s_res) > e);/* Выполнять до тех  пор, пока результат не будет удовлетворять допустимой ошибке */

Рекомендуем скачать другие рефераты по теме: реферат связь, рим реферат.





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