Помогиет найти ошибку в коде.

[sinh]

Писал прогу по следующим алгоритмам, но результаты совершенно разные, хотя решают одну и туже задачу. На С++ пишу 1-й раз, подскажите, где допускаю ошибки.

#include "stdafx.h"
#include "math.h"
#include <iostream>
#include <stdio.h>
 
 
double signum(double a)
{
    if (a > 0)
        a = 1;
    if (a < 0)
        a = -1;
    if (a == 0)
        a = 0;
    return a;
}
 
double function(double x, double w)
{
    double f = sinh(w*x);
    return f;
}
 
double vegstein(double x_0, double w, double e)
{
    double x_1 = function(x_0, w);
    double x_00 = x_0;
    double x_11 = x_1;
    int k = 0;
    shag2:
    double x_2 = function(x_11, w);
    if (abs(x_2-x_11) > e)
        {
            double x_22 = (x_2 * x_00 - x_1 * x_11)/(x_2 + x_00 - x_1 - x_11);
            x_00 = x_11;
            x_1 = x_2;
            x_11 = x_22;
            k++;
            goto shag2;
        }
    else 
        {
            double ksi = x_2;
            std::cout << "ksi = " << ksi << "\n";
            std::cout << "k = " << k << "\n";
        }
    getchar();
    return 0;
}
 
double eitken(double x_0, double  w, double e)
{
    double x_1 = function(x_0, w);
    double x_2 = function(x_1, w);
    int k = 0;
    shag2:
    double x_22 = (x_0*x_2 - x_1*x_1)/(x_2-2*x_1+x_0);
    k++;
    double x_3 = function(x_22, w);
    if (abs(x_3-x_22) > e)
        {
            x_0 = x_22;
            x_1 = x_3;
            x_2 = function(x_1, w);
            goto shag2;
        }
    else 
        {
            double ksi = x_3;
            std::cout << "ksi = " << ksi << "\n";
            std::cout << "k = " << k << "\n";
        }
    getchar();
    return 0;
}
 
 
int main(int argc, _TCHAR* argv[])
{
    double x_0 = 1.0;
    double w = 0.5;
    double e = powl(10,-15);
    printf("Metod Vegsteina: \n");
    vegstein(x_0, w, e);
    printf("Metod Eitkena: \n");
    eitken(x_0, w, e);
    return 0;
}

[Sciv]

double signum(double a)
{
    if (a > 0)
        a = 1;
    if (a < 0)
        a = -1;
    if (a == 0)
        a = 0;
    return a;
}

можно заменить на:

double signum(double a)
{
    if (a!=  0)   a = a/fabs((float)a); 
    return a;
}

Но это так, к слову.

[Гром]

Если уж на то пошло, то сигнум лучше так записывать:

int signum (double a)
{
if (a > 0)
 return 1;
if (a < 0)
 return -1;
return 0;
}

Аналогично function можно в одну строчку, сразу ретурном без ввода лишней переменной.
Также goto по-хорошему бы заменить на цикл while. А для чисел с плавающей точкой надо использовать не abs, а fabs.
Навскидку несоответствий не нашел.

[Sciv]

А проблема у Вас, скорее всего, как раз в модулях чисел: функция abs применяется только с целочисленными аргументами, то есть abs(-5.5634)=abs(5)=5;

Нужно менять:

fabs((float)(x_2-x_11)) > e

и

fabs((float)(x_3-x_22))>e

UPD. Ну вот, уже и ответили

[Гром]

Приводить к float-у необязательно (даже не стоит):

     double fabs (      double x );
      float fabs (       float x );
long double fabs ( long double x );

Ну и вообще тогда. Раз используете цпп-шные потоки, то зачем вам принтф? Выкиньте его вместе с stdio.h (точнее, cstdio).

#include "math.h"

тоже замените на

#include <cmath>

Это все не ошибки, но признак не лучшего стиля и устаревшие (хотя и все еще поддерживаемые) конструкции.

[sinh]

cпасибо большое, про fabs не знал.

[Somebody]

Если на то пошло,

int signum(double a)
{
    return (a > 0) - (a < 0);
}

[sinh]

изменение abs на fabs ничего не поменяло, результаты все равно значительно отличаются. Мб до нужной точности сравнение не происходит? И еще в Turbo C++ и Visual Studio выдает разные результаты))

[Гром]

Ищите багу. Подставьте получающиеся результаты в выражение x = phi(x), чтобы узнать, считает ли хотя бы одна правильно. Выводите на экран все промежуточные значения, чтобы проверить вручную хотя бы первые несколько итераций. Замените функцию на какую-нибудь попроще - например, на phi(x) = 2*x (ответ - x = 0). И так далее.

[sinh]

Тему можно закрывать.