Advertisement
Guest User

Untitled

a guest
Oct 8th, 2015
93
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.45 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2.  
  3. using namespace std;
  4.  
  5. #define lol long long
  6. #define fi first
  7. #define se second
  8. #define pb push_back
  9. #define sz(s) (lol)s.size()
  10. #define must ios_base::sync_with_stdio(0)
  11.  
  12. #define inp(s) freopen(s, "r", stdin)
  13. #define out(s) freopen(s, "w", stdout)
  14.  
  15. #define TABLE22
  16.  
  17. #ifdef TABLE21
  18. double m1 = 51.1, m2 = 49.8;
  19. double v10[] = {0.45, 0.47, 0.49, 0.48, 0.48};
  20. double v1[] = {0.22, 0.22, 0.24, 0.23, 0.23};
  21. #else
  22. double m1 = 51.1, m2 = 104;
  23. double v10[] = {0.49, 0.48, 0.47, 0.50, 0.48};
  24. double v1[] = {0.14, 0.14, 0.13, 0.13, 0.13};
  25. #endif
  26. double p10[5];
  27. double P[5];
  28. double p[5];
  29. double wE[5];
  30. double wT;
  31. double sp;
  32. double ssp;
  33. double swE;
  34. double sswE;
  35. double talpha = 2.055;
  36.  
  37. int main() {
  38.     for(int i = 0; i < 5; i++) {
  39.         p10[i] = v10[i] * m1;
  40.         P[i] = (m1 + m2) * v1[i];
  41.         p[i] = P[i] / p10[i] - 1;
  42.         wE[i] = ((m1 + m2) * v1[i] * v1[i]) / (m1 * v10[i] * v10[i]) - 1;
  43.         sp += p[i];
  44.         swE += wE[i];
  45.         cout << p10[i] << ' ' << P[i] << ' ' << p[i] << ' ' << wE[i] << '\n';
  46.     }
  47.     wT = -m2 / (m1 + m2);
  48.     sp /= 5.;
  49.     swE /= 5.;
  50.     cout << sp << ' ' << swE << "\n";
  51.     for(int i = 0; i < 5; i++) {
  52.         ssp += (p[i] - sp) * (p[i] - sp);
  53.         sswE += (wE[i] - swE) * (wE[i] - swE);
  54.     }
  55.     ssp = talpha * sqrt(ssp / 20.);
  56.     sswE = talpha * sqrt(sswE / 20.);
  57.     cout << ssp << ' ' << sswE << "\n";
  58.     cout << "wT=" << wT << "\n";
  59.     return 0;
  60. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement