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

Проверка метода Гаусса - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Объяснить работу программы http://www.cyberforum.ru/cpp-beginners/thread253775.html
Во введенном слове подсчитать количество различных пар букв. Например, в слове babacabacd 5 различных пар букв, ba – 3 пары, ac – 2 пары. #pragma argsused int main(int argc, char* argv) { { ...
C++ Вычислить произведение чисел кратных 7 на промежутке от 1 до 100. Вычислить произведение чисел кратных 7 на промежутке от 1 до 100. http://www.cyberforum.ru/cpp-beginners/thread253763.html
Вычислить сумму чисел от 1 до заданного числа C++
Вычислить сумму чисел от 1 до заданного числа!
Вывести фамилию 10 раз (FOR) C++
Вывести на экран свою фамилию 10 раз (FOR) Пожалуйста посмотрите !!!!!! Спасибо)
C++ Кодировка в разных ОС http://www.cyberforum.ru/cpp-beginners/thread253750.html
ВОт кусок программы ( консольная игра): #include<iostream> #include<string> #include"City.h" #include"city/build.h" int main(int argc, char *argv) { string n; setlocale(LC_ALL,...
C++ Числа А как создать входной файл для программы на Паскале? К примеру, необходимо обработать 1000 чисел, неужели все вводить вручную? Если нет, повторяю: как же создать этот входной файл с необходимой... подробнее

Показать сообщение отдельно
Jupiter
Каратель
Эксперт С++
6556 / 3977 / 227
Регистрация: 26.03.2010
Сообщений: 9,273
Записей в блоге: 1
Завершенные тесты: 2
08.03.2011, 16:16
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
#include <iostream>
 
void ReadSys(double **A, double *B, size_t dim)
{
    for(size_t i = 0; i < dim; i++)
    {
    for(size_t j = 0; j < dim; j++)
    {
        std::cout << "A[" << i << "][" << j << "] = ";
        std::cin >> A[i][j];
    }
    std::cout << "B[" << i << "] = ";
    std::cin >> B[i];
    }
}
 
void swap(double &a, double &b)
{
    double tmp = a;
    a = b;
    b = tmp;
}
 
void Gauss(double **A, double *B, size_t dim)
{
    for (size_t k = 0; k < dim; ++k)
    {
        if (A[k][k] == 0.0)
        {
            for (size_t i = k + 1; i < dim; ++i)
            {
                if (A[i][k] != 0.0)
                {
                    for (size_t j = 0; j < dim; ++j)
                        swap(A[i][j], A[k][j]);
                    swap(B[i], B[k]);
                    break;
                }
            }
        }
        double t = A[k][k];
        for (size_t j = 0; j < dim; ++j)
            A[k][j] /= t;
        B[k] /= t;
        for (size_t i = k + 1; i < dim; ++i)
        {
            double tmp = A[i][k];
            for (size_t j = 0; j < dim; ++j)
                A[i][j] -= tmp * A[k][j];
            B[i] -= tmp * B[k];
        }
    }
    for (size_t k = dim - 1; k > 0; --k)
    {
        for (size_t i = k - 1; i + 1 > 0; --i)
        {
            double tmp = A[i][k];
            for (size_t j = 0; j < dim; ++j)
                A[i][j] -= tmp * A[k][j];
            B[i] -= tmp * B[k];
        }
    }
}
 
int main()
{
    const size_t N = 10;
    double ** A = new double* [N];
    for(size_t i = 0; i < N; i++)
        A[i] = new double [N];
    double * B = new double [N];
 
    ReadSys(A, B, N);
    Gauss(A, B, N);
 
    for(size_t i = 0; i < N; i++)
        std::cout << "X[" << i << "] = " << B[i] << std::endl;
 
    std::cin.get();
    std::cin.get();
    for (int i = 0; i < N; ++i)
        delete [] A[i];
    delete [] A;
    delete [] B;
    return 0;
}
1
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru