#include #include using namespace std; void main() { double pi = 3.141591; double g = 9.81; double m = 0.00054; double Ro1 = 11000; double Ro2 = 920; double mu = 0.95; double v = 0.6; double r = pow(3 * m / (4 * pi * Ro1), 1.0 / 3); double V = 4 * pi * pow(r, 3) / 3; double k1 = 6 * pi * mu * r; double Fa = Ro2 * g * V; double Ft = m * g; double h = 0.04; double dt = 0.001; double t = 0; while (h >= 0) { double Fc = k1 * v; double Fnet = Ft - Fc - Fa; double a = Fnet / m; v += a * dt; h -= v * dt; t += dt; // cout << "t = " << t << ", x = " << h << ", v = " << v << ", a = " << a << endl; cout << h << endl; } }