Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- double ObGaussa(double a, double b, int w, double TabT[], double TabA[])
- {
- double wyn = 0;
- for(int i = 0; i < w; i++)
- {
- wyn += TabA[i]*(1/(((b-a)/2)*TabT[i]+((b+a)/2)));
- }
- wyn *= ((b-a)/2);
- return wyn;
- }
- double W1Gaussa(double a, double b)
- {
- double TabT[1];
- double TabA[1];
- TabA[0] = 2;
- TabT[0] = 0;
- return ObGaussa(a, b, 1, TabT, TabA);
- }
- double W2Gaussa(double a, double b)
- {
- double TabT[2];
- double TabA[2];
- TabA[0] = 1;
- TabA[1] = 1;
- TabT[0] = -0.57735;
- TabT[1] = 0.57735;
- return ObGaussa(a, b, 2, TabT, TabA);
- }
- double W3Gaussa(double a, double b)
- {
- double TabT[3];
- double TabA[3];
- TabA[0] = 0.55556;
- TabA[1] = 0.55556;
- TabA[2] = 0.88889;
- TabT[0] = -0.7746;
- TabT[1] = 0.7746;
- TabT[2] = 0;
- return ObGaussa(a, b, 3, TabT, TabA);
- }
- double W4Gaussa(double a, double b)
- {
- double TabT[4];
- double TabA[4];
- TabA[0] = 0.34785;
- TabA[1] = 0.34785;
- TabA[2] = 0.65215;
- TabA[3] = 0.65215;
- TabT[0] = -0.86114;
- TabT[1] = 0.86114;
- TabT[2] = -0.33998;
- TabT[3] = 0.33998;
- return ObGaussa(a, b, 4, TabT, TabA);
- }
- double W5Gaussa(double a, double b)
- {
- double TabT[5];
- double TabA[5];
- TabA[0] = 0.23693;
- TabA[1] = 0.23693;
- TabA[2] = 0.47863;
- TabA[3] = 0.47863;
- TabA[4] = 0.56889;
- TabT[0] = -0.90618;
- TabT[1] = 0.90618;
- TabT[2] = -0.53847;
- TabT[3] = 0.53847;
- TabT[4] = 0;
- return ObGaussa(a, b, 4, TabT, TabA);
- }
- void Gaussa()
- {
- int w;
- double a, b, wyn;
- cout << "Podaj a: ";
- cin >> a;
- cout << "Podaj b: ";
- cin >> b;
- cout << "Podaj w (1 - 5): ";
- cin >> w;
- switch (w)
- {
- case 1 : wyn = W1Gaussa(a, b); break;
- case 2 : wyn = W2Gaussa(a, b); break;
- case 3 : wyn = W3Gaussa(a, b); break;
- case 4 : wyn = W4Gaussa(a, b); break;
- case 5 : wyn = W5Gaussa(a, b); break;
- default : cout << "Brak implementacji powyzej 5 wezlow!" << endl; return; break;
- }
- cout << "Pole: " << wyn << endl;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement