Форум программистов, компьютерный форум CyberForum.ru

метод Эйткена-Стеффенсона - C++

Восстановить пароль Регистрация
 
fits
3 / 3 / 1
Регистрация: 08.12.2012
Сообщений: 47
28.05.2013, 18:19     метод Эйткена-Стеффенсона #1
добрый день,помогите,пожалуйста ,разобраться где ошибка в коде,он должен давать ответ ОТВЕТ: x=1.171229 на интервале [a=1, b=2]

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
36
37
38
39
#include "stdafx.h"
#include <iostream>
#include <conio.h>
#include <math.h>
#include <stdio.h>
using namespace std; 
 
    
 
// this function returns value of f(x)
double f(double x)
{
  return x-sin(x)-0.3;
}
 
 
double Solve(double x0, double epsilon)
{
  double x1,x2,xn,znam,a,b,c;
  do
  {
    x0=0;
    xn=x0;
    x1=f(x0);
    x2=f(x1);
    znam=x0-2*x1+x2;
    x0=(x0*x2-x1*x1)/znam;
  } while (fabs(x0-xn)>epsilon && znam!=0);
  return x0;
}
 
int main()
{   
    double  epsilon=1e-5;
 
    printf("x-sin(x)-0.3, [a; b]=[1.0; 2.0]");
    printf("\n x0 = %f",Solve(1.0,epsilon));
    getch();
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
28.05.2013, 18:19     метод Эйткена-Стеффенсона
Посмотрите здесь:

Объясните метод Эйткена C++
Метод деления отрезка пополам для решения нелинейных уравнений (метод дихотомии) C++
C++ Метод медиан из трех элементов VS улучшенный быстрый метод сортировки(метод Бентли-Макилроя)
метод Эйткена-Стеффенсона C++
C++ Производный класс: метод возведения в произвольную степень, и метод для вычисления логарифма числа
Класс vector (поля: координаты, 2 конструктора, метод нахождения длины вектора и метод вывода координат на экран) C++
C++ СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя
Нахождения корней уравнения: метод половинного деления (бисекции) или метод хорд C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 18:27. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru