Алгоритм сделать явью

@vladrrom · Регистрация: 13.10.2016

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

Доброго времени суток. Снова пишу сюда.. Помогите ребята. Нужно составить код алгоритма, но не имею понятия каким образом сделать одну вещицу. Расскажу по порядку..
Алгоритм:

Алг. метода Ньютона:
Шаг 1. Задать функцию f(x), числа a и b отрезка [a,b], число ε - точность алгоритма, x₀ ∈ [a,b] - начальное приближение к искомой точке, k=0.
Шаг 2. Найти точку x_k+1 по формуле
x_k+1 = x_k - f'(x_k)/f''(x_k).
Перейти к шагу 3.
Шаг 3. Проверка на окончание поиска. Если |f'(x_k+1)| <= ε, то положить x*=x_k+1 и завершить алгоритм. Иначе, положить k=k+1 и перейти к шагу 2.

Алгоритм я понимаю как строить, но как задать переменную x_k и к чему её приравнять для меня остается загадкой. Сначала думал произвести это через массив, но так до реализации и не дошло..
Код алгоритма:

C++

#include <iostream>
#include <stdio.h>
using namespace std;
 
double fs (double x) //первая производная
{
    return (exp(x) - 1/x*x);
}
double fss (double x) //вторая производная
{
    return (exp(x) + 2/(x*x*x));
}
int main()
{
    double b, a, xs, eps; // xs - это x'
    eps = 0.1;
        int k=0;
    // объявление x0
    a = 0.5;
    b = 1.5;
    flag = true;
    while (flag)
    {
        // формула
        if (abs(fs(/* xk+1 */)) <= eps)
        {
            xs = /* xk+1 */;
            flag = false;
            break;
        }
        else
        {
            k=k+1;
            continue;
        }
    }
    cout << "Method Newton" << endl;
    cout << "x* = " << xs << "; \t" << "f* = " << f(xs) << endl;
    system("pause");
    return 0;
}

	31.03.2018, 14:06