6inline double formula_d(
double d,
double t,
double da,
double ta) {
7 return d + da * exp((ta - t) / 200);
10inline double formula_p(
double s,
double u,
double t,
double d) {
11 const double kM = 1 / (1 - exp(-0.005));
12 double m = s - (s - u) * pow((1 - exp(-t / 10000)), 10);
13 return (d < 20) ? m + (0.5 - m) * (d / kM)
14 : m + (1 - m) * (pow(4, (d / kM)) - 1) / 3;
double formula_p(double s, double u, double t, double d)
Definition dynamics.h:10
double formula_d(double d, double t, double da, double ta)
Definition dynamics.h:6