Нахождение корней уравнений. Метод Ньютона

@NaiMI · Регистрация: 13.03.2013

Студворк — интернет-сервис помощи студентам

Код взят из книги Numerical Recipes in C The Art of Scientific Computing Second Edition авторов William H. Press(Harvard-Smithsonian Center for Astrophysics),Saul A. Teukolsky(Department of Physics, Cornell University),William T. Vetterling(Polaroid Corporation),Brian P. Flannery(EXXON Research and Engineering Company)

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
#include <cmath>
#include <stdio.h>
#include <iostream>
 
#define JMAX 20
 
float rtnewt(void (*funcd)(float, float *, float *), float x1, float x2, float xacc)
{
    int j;
    float df,dx,f,rtn;
    rtn=0.5*(x1+x2); 
    for (j=1;j<=JMAX;j++)
    {
        (*funcd)(rtn,&f,&df);
        dx=f/df;
        rtn -= dx;
        if ((x1-rtn)*(rtn-x2) < 0.0)
            printf("Jumped out of brackets in rtnewt \n");
        if (fabs(dx) < xacc) return rtn;
    }
    printf("Maximum number of iterations exceeded in rtnewt \n");
    return 0.0;
}
 
void myf(float x, float* f, float* df) 
 {
    *f = exp(-pow(x,2))-0.01;
    *df = -2*x*exp(-pow(x,2));
 }
int main()
 {
    printf("%.10f \n", rtnewt(myf, 1.0f, 3.0f, 0.0001f));
    system("pause");
    return 0;
 }

@volvo · 29.10.2014, 19:26

Прекрасно. И что не так с этим кодом? И зачем он в разделе по Билдеру?

Новые блоги и статьи Все статьи Все блоги /
Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем	10 пpимет, которые всегда сбываются Maks 31.03.2026 1. Чтобы, наконец, пришла маршрутка, надо закурить. Если сигарета последняя, маршрутка придет еще до второй затяжки даже вопреки расписанию. 2. Нaдоели зима и снег? Не надо переезжать. Достаточно. . .	Перемещение выделенных строк ТЧ из одного документа в другой Maks 31.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .
Functional First Web Framework Suave DevAlt 30.03.2026 Sauve. IO Апнулись до NET10. Из зависимостей один пакет, работает одинаково хорошо как в режиме проекта так и в интерактивном режиме. из сложностей - чисто функциональный подход. Решил. . .	Автоматическое создание документа при проведении другого документа Maks 29.03.2026 Реализация из решения ниже выполнена на нетиповых документах, разработанных в конфигурации КА2. Есть нетиповой документ "ЗаявкаНаРемонтСпецтехники" и нетиповой документ "ПланированиеСпецтехники". В. . .	Настройка движения справочника по регистру сведений Maks 29.03.2026 Решение ниже реализовано на примере нетипового справочника "ТарифыМобильнойСвязи" разработанного в конфигурации КА2, с целью учета корпоративной мобильной связи в коммерческом предприятии. . . .	Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .

@volvo Супер-модератор 33403 / 21513 / 8236 Регистрация: 22.10.2011 Сообщений: 36,914 Записей в блоге: 12
	29.10.2014, 19:26
	Прекрасно. И что не так с этим кодом? И зачем он в разделе по Билдеру? 0