Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <iostream>
- #include<cmath>
- using namespace std;
- int main()
- {
- const int n = 15; //количество итераций
- const float g = 9.8, k = 0.3;
- float alpha; // угол в градусах
- float alfa; // угол в радианах
- float vzero, vx, vy;
- float x[n+1], y[n+1];
- float tmax; // время полета
- float h; // интервалы
- float eps[n + 1] = {0}; // погрешности
- float yeps;
- cout<<"v0 = ";
- cin>>vzero;
- cout<<"alpha = ";
- cin>>alpha;
- alfa = alpha*M_PI/180;
- tmax = (2*vzero*sin(alfa))/g;
- h = tmax/n;
- vx = vzero*cos(alfa);
- vy = vzero*sin(alfa);
- //cout<<x[0]<<" "<<y[0]<<" "<<eps[0]<<endl;
- for(int i = 1; i < n + 1; i++)
- {
- x[i] = x[i - 1] + vx*h;
- y[i] = y[i - 1] + vy*h - g*h*h/2;
- yeps = y[i - 1] + vy*(h/2) - g*(h/2)*(h/2)/2;
- eps[i] = abs(yeps - y[i])/abs(yeps);
- vy -= g*h;
- //cout<<x[i]<<" "<<y[i]<<" "<<eps[i]<<endl;
- }
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment