КАТЕГОРИИ: Архитектура-(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) |
Вычисление определенного интегралаСледующий пример иллюстрирует применение стека для решения задачи, в которой часть работы выполняется по тому же алгоритму, что и вся работа. Каждое значение интеграла будем вычислять по формуле прямоугольников (хотя применяемый метод безразличен и можно было бы использовать любой). Вычисляем два очередных приближения интеграла – - как площадь одного прямоугольника: - как площадь двух прямоугольников: Если два приближения отличаются друг от друга более чем на назначенную меру погрешности eps, делим интервал интегрирования (a,b) пополам, и левую половину сразу обработаем тем же алгоритмом, а обработку правой отложим на потом – в стек. В противном случае, I2 – окончательный ответ. Прибавим его к уже вычисленной части результата и поинтересуемся, не осталось ли отложенной на потом работы в стеке. Если ДА, то возьмемся за ее выполнение, если НЕТ – решение получено. Текст функции, использующей ручное ведение стека для реализации рекурсии, приведен ниже:
#include <math.h> #include <values.h> struct INTERVAL{ double a,b; }; const int MAXSTACK=100; //---------------------------------------- double Integral(double a, double b, double (*f)(double), double eps){ // вычисление определенного интеграла от a до b с точн. eps // от ф-ии f по формуле прямоугольников double I1,I2; // два последовательных приближения double d; double Result=0; int v=-1; // указатель стека INTERVAL Stack[MAXSTACK]; Again: d=b-a; I1=(*f)(a)*d; // счет на одном интервале d=(b-a)/2; // ширина интервала I2=((*f)(a)+(*f)((a+b)/2))*d; // счет на 2х интервалах if(fabs(I1-I2)<eps){ Result+=I2; if(v>=0){ // если стек не пуст - взять из стека a=Stack[v].a; b=Stack[v].b; v--; goto Again; } else { // стек пуст - решение получено return Result; } } else { // правую часть - в стек // а левой займемся немедленно v++; if(v>=MAXSTACK){ MessageBox(0,"Переполнение стека","Ошибка",MB_OK); return MAXDOUBLE; // см. values.h } else { Stack[v].a=(a+b)/2; Stack[v].b=b; // a остается без изменений b=(a+b)/2; goto Again; } } }
Дата добавления: 2014-12-08; Просмотров: 420; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |