VPF::Проблема с Методом Холецкого. - Форум программистов

slevin45

Дата 4.3.2014, 16:27 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 3
Регистрация: 4.3.2014

Репутация: нет
Всего: нет

Здравствуйте!
Помогите,пожалуйста,найти ошибку в коде.Выдаёт неправильные корни.Вот блок схема и формула,по которой находил элементы S и d.

Код

#include<iostream>
#include<math.h>
using namespace std;
int main (void)
{
    double a[6][6],b[6],x[6],y[6],delta,delta2=0,delta3=0,del2,del3,delta4=0,del4,s[6][6],d[6],H;
    int i,j,n,k;
    cout<<"vvedite poryadok matrici";
    cin>>n;
    n=n+1;
    a[1][1]=-4.88; a[1][2]=1;  a[1][3]=0;a[1][4]=0; a[1][5]=0;
    a[2][1]=1;     a[2][2]=-2; a[2][3]=1; a[2][4]=0;a[2][5]=0;
    a[3][1]=0;     a[3][2]=1;  a[3][3]=-2;a[3][4]=1; a[3][5]=0;
    a[4][1]=0;      a[4][2]=0;  a[4][3]=1; a[4][4]=-2;a[4][5]=1;
    a[5][1]=0;      a[5][2]=0;  a[5][3]=0; a[5][4]=1; a[5][5]=-4.88;
    b[1]=0;b[2]=-8;b[3]=-8;b[4]=-8;b[5]=0;
    for(k=1;k<n;k++)//нахождение элементов S и d;
    {
        for(i=1;i<=k-1;i++)
        {
            del2=d[i]*pow(fabs(s[i][k]),2);
            delta2=delta2+del2;
        }
        if((a[k][k]-delta2)>0)
            d[k]=1;
            if((a[k][k]-delta2)<0)
                d[k]=-1;
                s[k][k]=pow(fabs(a[k][k]-delta2),1/2);
                for(j=k+1;j<n;j++)
                {
                    for(i=1;i<=k-1;i++)
                    {
                        del4=d[i]*s[i][k]*s[i][j];
                        delta4=delta4+del4;
                    }
                    s[k][j]=(a[k][j]-delta4)/(s[k][k]*d[k]);
                }
    }
    for(k=1;k<n;k++)
    {
        delta=a[k][k];
        if(k!=1)
        {
            for(i=1;i<k-1;i++)
            {
                delta=delta-d[i]*pow(s[i][k],2);
            }
        }
        if(delta>0)
        d[k]=1;
        if(delta<0)
            d[k]=-1;
            s[k][k]=pow(fabs(delta),1/2);
                for(j=k+1;j<n;j++)
                {
                    delta=a[k][j];
                    if(k!=1)
                    {
                        for(i=1;i<k-1;i++)
                        {
                            delta=delta-d[i]*s[i][k]*s[i][j];
                        }
                    }
                    s[k][j]=delta/(s[k][k]*d[k]);
                }
    }
    y[1]=b[1]/(s[1][1]*d[1]);
    for(i=2;i<n;i++)
    {
        delta=b[i];
        for(k=1;k<=i-1;k++)
        {
            delta=delta-d[k]*y[k]*s[k][i];
        }
        y[i]=delta/(s[i][i]*d[i]);
    }
    x[n-1]=y[n-1]/s[n-1][n-1];
    for(i=n-2;i>=1;i--)
    {
        delta=y[i];
        for(k=i+1;k<n;k++)
        {
            delta=delta-x[k]*s[i][k];
        }
        x[i]=delta/s[i][i];
    }
    cout<<endl;
    for(i=1;i<n;i++)
    {
        cout<<x[i];
        cout<<endl;
    }
    cout<<"nevyazka";
     for(i=1;i<n;i++)
                  {
                      double z=0;
                        for(int j=1;j<n;j++)
                        {
                        z=z+(a[i][j]*x[j]);
                        }
                            H=fabs(b[i]-z);
                            cout<<endl;
                             cout<<"Nevyazka"<<H;
     }
    system("pause");
    return(0);
}

borisbn

Дата 5.3.2014, 09:34 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 4875
Регистрация: 6.2.2010
Где: Ростов-на-Дону

Репутация: 21
Всего: 135

Цитата(slevin45 @ 4.3.2014, 16:27

)

del2=d[i]*pow(fabs(s[i][k]),2);

первое, что бросилось в глаза - массив s используется, но нигде не заполняется.

Далее, в строках 24, 26, 51 (м.б. где-то ещё) какое-то странное форматирование.. Ты нигде не забыл поставить фигурные скобки ?

--------------------

Женщины отличаются от программистов тем, что у них чары состоят из стрингов

1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Для новичков \| Следующая тема »