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

LU - разложение (выдает NaN) - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Примерчики ООП простые http://www.cyberforum.ru/cpp-beginners/thread747231.html
Всем добрый вечер Вот задачка Необходимо создать 3 класса Авто, Запчасти и Калькулятор в Классе Авто должны быть массивом запчасти и имя авто в Классе Запчасти Название запчасти и стоимость
C++ Массивы. Составить функцию для построения вектора по правилу Задали матрицу А = {aij} размерности n x n. Составить функцию для построения вектора B = {bi} по следующему правилу: если aii < 0, то bi взять равным сумме элементов i - й строки, иначе - сумму абсолютных величин i - го столбца. Затем в матрице А расположить элементы каждой строки в порядке возрастания. http://www.cyberforum.ru/cpp-beginners/thread747197.html
C++ Блок схема. С++
Помогите с блок схемой к этому коду #include <iostream.h> #include <conio.h> #include <math.h> #include <stdlib.h> void main() { clrscr(); int n,i=0,j=0,ch=0,z=0,k=0, buf;
Сложная обработка элементов файлов. Найти сумму положительных компонент файла, расположенных до максимальной компоненты... C++
Помогите пожалуйста! Создать числовой файл F . Найти сумму положительных компонент файла, расположенных до максимальной компоненты. Найденной суммой заменить отрицательные компоненты, строящие после максимума. Файл и сумму распечатать (как до, так и после преобразования). Добавлено через 3 часа 0 минут Пожалуйста Добавлено через 3 часа 11 минут Очень нужно, а я по файлам не шарю..
C++ Сортировка массива структур по одному из полей http://www.cyberforum.ru/cpp-beginners/thread747162.html
Добрый вечер. Есть такое задание: Организуйте сортировку массива студентов по номеру года обучения. Используйте стандартную функцию qsort. Указание: при сравнении понадобится разыменовывать указатель, а затем обращаться к полю структуры студента. Это могло бы выглядеть так: (*p).year. Однако обращение к полям структуры по указателю на структуру настолько частое действие в языке C, что для этого...
C++ Запись и считывание с файла с проверкой Написать программу, которая считывает из файла и выводит на экран только строки, не содержащие двузначныx чисел. Вот я написал программу, все работает. #include <iostream> #include <fstream> #include <conio.h> using namespace std; int main () { подробнее

Показать сообщение отдельно
ISunRiseI
0 / 0 / 0
Регистрация: 24.12.2012
Сообщений: 4
26.12.2012, 21:16     LU - разложение (выдает NaN)
Здравствуйте. Помогите исправить программу. Проблема заключается в том, что некоторые значение получаются NaN/
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
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
#include<iostream.h>
#include<conio.h>
#include<iomanip.h>
int main()
{
    int i,j,k,n;
    float a[10][10],b[10],x[10],l[10][10],z[10],u[10][10],sum;
    cout<<"Enter number of unknowns:";
    cin>>n;
    cout<<"Enter co-efficients "<<endl;
    for(i=0;i<n;i++)
    for(j=0;j<=n;j++)
    {
        if(j<n) cin>>a[i][j];
        else cin>>b[i];
    }
    for(i=0;i<n;i++)
    for(j=0;j<n;j++)
        if(i>j) u[i][j]=0;
        else if(i==j) l[i][j]=1;
        else l[i][j]=0;
    //--------------------------------------------------------------------
    //AX=B
    //Decomposition of A into L and U   
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            sum=0;
            if(i<=j)
            {
                for(k=0;k<n;k++) 
                    if(k!=i)
                        sum=sum+l[i][k]*u[k][j];
                u[i][j]=a[i][j]-sum;
            }
            else 
            {
                for(k=0;k<n;k++)
                    if(k!=j)
                        sum=sum+l[i][k]*u[k][j];
                    l[i][j]=(a[i][j]-sum)/u[j][j];
            }
        }
    }
    cout.precision(3);
    cout<<"\n\n  L= ";
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
             cout<<setw(6)<<l[i][j];
        cout<<endl;
        cout<<"     ";
    }
    cout<<"\n\n  U= ";
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
             cout<<setw(6)<<u[i][j];
        cout<<endl;
        cout<<"     ";
    }
    //--------------------------------------------------------------------
    //Now LUX=B
    //i.e LZ=B
    //Obtaining Z by forward subtitiution 
    z[0]=b[0]/l[0][0];
    for(i=1;i<n;i++)
    {
        sum=0;
        for(j=0;j<i;j++)
        {
            sum=sum+z[j]*l[i][j];
            
        }
        z[i]=b[i]-sum;
    }
    
    cout<<"\n\n  Z= ";
    for(i=0;i<n;i++)
    {
        cout<<setw(6)<<z[i];
        cout<<endl;
        cout<<"     ";
    }
    //--------------------------------------------------------------------
    //Now UX=Z
    //Finding X by backward substitution
    x[n-1]=z[n-1]/u[n-1][n-1];
    for(i=n-2;i>=0;i--)
    {
        sum=0;
        for(j=n-1;j>i;j--)
            sum=sum+x[j]*u[i][j];
        x[i]=(z[i]-sum)/u[i][i];
 
    }
    cout.precision(6);
    cout<<"\nThe solution is :"<<endl;
    for(i=0;i<n;i++)
        cout<<"\n x"<<i+1<<" = "<<x[i];
    cout<<endl;
    return 0;
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 01:28. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru