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

Решение СЛАУ методом простых итераций - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ изобразить меню функциональная клавиша в виде строки внизу экрана и окно 2 с рамкой http://www.cyberforum.ru/cpp-beginners/thread418448.html
прошу помочь, т.к во6ще не понимаю этого:wall:
C++ потоковый вывод. язык С. ребята помогите пожалуйста, разобраться с задачей, завтра здавать, я основную часть сделала оставалось только потоковый уровень, проблема в том что нужно написать не структурный массив, а считать с... http://www.cyberforum.ru/cpp-beginners/thread418439.html
C++ Найти минимальное время прохождения трассы.
Трасса для соревнований задана в виде n -угольника , в одной из вершин которого находится место старта, а одна из сторон - линия финиша (место старта - не на линии финиша!). Путь по трассе...
перекомпиляция из VS 2008 в VS 2010 C++
помогите пожалуйста
C++ Вывести на печать все пифагоровы числа, не превышающие 55 http://www.cyberforum.ru/cpp-beginners/thread418426.html
Программы с++ Прошу помощь в написании этих программ ( (в Паскале,по возможности) и С++), с комментариями : 1) 1)Вывести на печать все пифагоровы числа,не превышающие 55.(Пиф цифры это тройки...
C++ Чтение двоичного файла Есть файл, без типа,просто файл. Содержит набор чисел (размером по 4 б - int ) нужно его прочитать. Пишу: FILE *fp; if ((fp = fopen("ITI","rb"))==NULL) //ITI - название этого файла,хранится... подробнее

Показать сообщение отдельно
case
7 / 7 / 2
Регистрация: 29.09.2010
Сообщений: 32
27.12.2011, 01:38
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
#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
#define n 4
 
void del(double **A, double *X, double *Xk,   double *b);
 
int main()
 
{
 
    int i,j;
    double **A;     //matrix equations
    double *X;      //communications solution and iteration (k+1)
    double *Xk;     //iteration (k)
    double *b;      //free to members of
    double sumd = 0;//additional variables for computation
    double eps;
 
//--------------dynamic memory allocation------------------------------//
    A  = new double* [n];
    X  =  new double [n];
    Xk =  new double [n];
    b  =  new double [n];
    for(i = 0; i < n; i++)
        A[i] = new double [n];
//-------------input system of equations-------------------------------//
    printf("Type elements in the variables in the system: \n");
    for(i = 0; i < n; i++)
       for(j = 0; j < n; j++)
        scanf("%lf",&A[i][j]);
 
 
    printf("Members enter free system: \n");
    for(i = 0; i < n; i++)
        scanf("%lf",&b[i]);
    printf("enter the calculation accuracy: \n");
        scanf(" %lf",&eps);
//-------------test system on the diagonal elements--------------------//
    for(i = 0; i < n; i++)
       if (A[i][i] == 0) {del( A, X, Xk, b); printf("!On the main diagonal of the system is zero"); exit(0);}
//-------------checking for convergence of-----------------------------//
    for(i = 0; i < n; i++)
      {
       for(j = 0; j < n; j++)
          if (i != j)       sumd += fabs(A[i][j]);
       if (fabs(A[i][i]) < sumd)  {del( A, X, Xk, b);  printf("!Diverging methods for this system"); exit(0);}
       sumd = 0;
      }
//----- approximate solutions of the system-//
    for(i = 0; i < n; i++)
         X[i]=0.0;
    
//-------------calculation system of linear algebraic equations---------//
    int  count = 0, flag = 1;
    double x, zh = 0;
    for(i = 0; i < n; i++)
    Xk[i] = X[i];
    do
      {
       count++;
       for(i = 0; i < n; i++)
     {
       x=0;
       for(j = 0; j < n; j++)
          {
         if (i != j)  x += Xk[j] * A[i][j];
         if (i == j) zh = A[i][j];
          }
       x = (b[i] - x) / zh;
       X[i] = x;
       if ((fabs(X[i]-Xk[i]))<=eps) flag=0;
     }
     for(i = 0; i < n; i++)
        Xk[i] = X[i];
      }while(flag);
    printf("solution call method of iterations: \n");
    for(i = 0; i < n; i++)
       printf("x[%d]=%0.5lf\n",i+1,X[i]);
       printf("count=%d\n",count);
    
       del( A, X, Xk, b);
       return 0;
}
 
 
void del(double **A, double *X, double *Xk,   double *b)
  {
    for(int i = 0; i < n; i++)
      delete[]A[i];
     delete[]A;
     delete[]X;
     delete[]Xk;
     delete[]b;
  }
2
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru