Організація реклами

Дата добавления: 2015-09-15; просмотров: 540

#include "stdafx.h"

#include <cstdio>

#include <cstring>

#include <cstdlib>

#include <cmath>

#include <iostream>

#include <iomanip>

using namespace std;

long double xI[11] = {1.415, 1.420, 1.425, 1.430, 1.435, 1.440, 1.445, 1.450, 1.455, 1.460, 1.465};

long double yI[11] = {0.888551, 0.889599, 0.890637, 0.891667, 0.892687, 0.893698, 0.894700, 0.895693, 0.896677, 0.897653, 0.898619};

long double delta_1_yI[10] = {0.001048, 0.001038, 0.00103, 0.00102, 0.001011, 0.001002, 0.000993, 0.000984, 0.000976, 0.000966};

long double delta_2_yI[9] = {-0.0001, -0.000008, -0.00001, -0.000009, -0.000009, -0.000009, -0.000009, -0.000008, -0.00001};

long double delta_3_yI[8] = {0.000092, 0.000018, 0.000001, 0.0, 0.0, 0.0, 0.000001, -0.000002};

long double formula_newton(long double x, unsigned N)

{

N--;

long double q = (x - xI[N]) / 0.005;

return (yI[N] + q*delta_1_yI[N] + q*(q-1)*delta_2_yI[N]/2 + q*(q-1)*(q-2)*delta_3_yI[N]/6);

}

int main(int argc, char** argv)

{

cout << "x1\t" << fixed << setprecision(16) << formula_newton(1.4161, 1) << "\t" << "1.4146\t" << "1" << endl;

cout << "x2\t" << fixed << setprecision(16) << formula_newton(1.4625, 10) << "\t" << "1.4625\t" << "10" << endl;

cout << "x3\t" << fixed << setprecision(16) << formula_newton(1.4135, 1) << "\t" << "1.4135\t" << "1" << endl;

cout << "x4\t" << fixed << setprecision(16) << formula_newton(1.470, 11) << "\t" << "1.470\t" << "11" << endl;

cout << "\n\n\n";

system("pause");

return 0;

}

Вывод программы

Третий столбец – значение переменной

Второй столбец – соответствующее значение функции

Четвертый столбец – номер из таблицы функции (выбирался приближенный по значению)