Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <bits/stdc++.h>
- using namespace std;
- #define lol long long
- #define fi first
- #define se second
- #define pb push_back
- #define sz(s) (lol)s.size()
- #define must ios_base::sync_with_stdio(0)
- #define inp(s) freopen(s, "r", stdin)
- #define out(s) freopen(s, "w", stdout)
- #define TABLE22
- #ifdef TABLE21
- double m1 = 51.1, m2 = 49.8;
- double v10[] = {0.45, 0.47, 0.49, 0.48, 0.48};
- double v1[] = {0.22, 0.22, 0.24, 0.23, 0.23};
- #else
- double m1 = 51.1, m2 = 104;
- double v10[] = {0.49, 0.48, 0.47, 0.50, 0.48};
- double v1[] = {0.14, 0.14, 0.13, 0.13, 0.13};
- #endif
- double p10[5];
- double P[5];
- double p[5];
- double wE[5];
- double wT;
- double sp;
- double ssp;
- double swE;
- double sswE;
- double talpha = 2.055;
- int main() {
- for(int i = 0; i < 5; i++) {
- p10[i] = v10[i] * m1;
- P[i] = (m1 + m2) * v1[i];
- p[i] = P[i] / p10[i] - 1;
- wE[i] = ((m1 + m2) * v1[i] * v1[i]) / (m1 * v10[i] * v10[i]) - 1;
- sp += p[i];
- swE += wE[i];
- cout << p10[i] << ' ' << P[i] << ' ' << p[i] << ' ' << wE[i] << '\n';
- }
- wT = -m2 / (m1 + m2);
- sp /= 5.;
- swE /= 5.;
- cout << sp << ' ' << swE << "\n";
- for(int i = 0; i < 5; i++) {
- ssp += (p[i] - sp) * (p[i] - sp);
- sswE += (wE[i] - swE) * (wE[i] - swE);
- }
- ssp = talpha * sqrt(ssp / 20.);
- sswE = talpha * sqrt(sswE / 20.);
- cout << ssp << ' ' << sswE << "\n";
- cout << "wT=" << wT << "\n";
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement