daily pastebin goal
80%
SHARE
TWEET

Untitled

a guest Apr 17th, 2018 56 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. #include <iostream>
  2. #include <cmath>
  3. #include <conio.h>
  4.  
  5. using namespace std;
  6.  
  7. int main()
  8. {
  9.         int stop;
  10.         double a1, a2, a3;
  11.         double y0[3];
  12.         double tablica[1000][3];
  13.         double macierz[3][3];
  14.         double lambdaMax;
  15.         double wektor[3];
  16.  
  17.         y0[0] = 1; y0[1] = 0; y0[2] = 0;
  18.         cout << "Wprowadz wspolczynnik:" << endl;
  19.         for (int i = 0; i < 3; i++)
  20.         {
  21.             for (int j = 0; j < 3; j++)
  22.             {
  23.                 cout << "a" << i + 1 << j + 1 << " = ";
  24.                 cin >> macierz[i][j];
  25.             }
  26.  
  27.         }
  28.  
  29.         cout << "Podaj liczbe iteracji: ";
  30.         cin >> stop;
  31.        
  32.        
  33.  
  34.         tablica[0][0]=0.0;
  35.  
  36.         for(int i=0;i<=stop;i++)
  37.         {
  38.             if(i==0)
  39.             {
  40.                 for(int j=0;j<3;j++)
  41.                 {
  42.                     tablica[i][j]=y0[j];
  43.                 }
  44.             }
  45.             else
  46.             {
  47.                     for(int j=0;j<3;j++)
  48.                     for(int w=0;w<3;w++)
  49.                     {
  50.                     tablica[i][j] = tablica[i][j] + macierz[w][j]*tablica[i-1][w];
  51.                     }
  52.             }
  53.         }
  54.  
  55.  
  56.         a1 = tablica[stop][0]/tablica[stop-1][0];
  57.         a2 = tablica[stop][1]/tablica[stop-1][1];
  58.         a3 = tablica[stop][2]/tablica[stop-1][2];
  59.        
  60.         lambdaMax = (a1+a2+a3)/3;
  61.  
  62.         double pom = wektor[0];
  63.  
  64.         for(int i=0;i<3;i++)
  65.         {
  66.             wektor[i] = tablica[stop][i] / pom;
  67.         }
  68.  
  69.         cout << "Dominujaca wartosc wlasna zadanej macierzy wynosi: " << lambdaMax << endl;
  70.         cout << "Wektor wlasny (zwiazany z dominujaca wartoscia):" << endl;
  71.  
  72.         for(int i=0;i<3;i++)
  73.         {
  74.             cout << wektor[i] << endl;
  75.         }
  76.  
  77.     return 0;
  78. }
RAW Paste Data
Top