PIB

/**
 * UAM - Uniform Accelerated Motion
 */

#include <cstdio>
#include <cstdlib>
#include <vector>
#include <boost/numeric/odeint.hpp>

struct UAM {
  typedef std::vector<double> State;
  State a; // a_{x,y}
  UAM(const State &_a) : a(_a) {}
  void operator()(const State& x, State& dxdt, const double t) {
    dxdt[0] = x[2]; // v_x(t)
    dxdt[1] = x[3]; // v_y(t)
    dxdt[2] = a[0]; // a_x(t)
    dxdt[3] = a[1]; // a_y(t)
  }
};

int main(int argc, char *argv[])
{
  UAM::State x(4); // position_{x,y}(t), v_{x,y}(t)
  UAM::State a(2); // a_{x,y}(t)
  UAM::State t(2); // t_{0,1}
  double dt;
  if (argc != 10) {
    fprintf(stderr, "Usage: %s <x(t0)> <y(t0)> <vx(t0)> <vy(t0)> <ax(t)> <ay(t)> <t0> <t1> <dt>\n", argv[0]);
    return EXIT_FAILURE;
  }
  x[0] = atof(argv[1]); // x(t0)
  x[1] = atof(argv[2]); // y(t0)
  x[2] = atof(argv[3]); // v_x(t0)
  x[3] = atof(argv[4]); // v_y(t0)
  a[0] = atof(argv[5]); // a_x(t)
  a[1] = atof(argv[6]); // a_y(t)
  t[0] = atof(argv[7]); // t_0
  t[1] = atof(argv[8]); // t_1
  dt   = atof(argv[9]);
  
  boost::numeric::odeint::runge_kutta4<UAM::State> stepper;
  UAM uam(a);
  size_t steps = boost::numeric::odeint::integrate_const(stepper, uam, x, t[0], t[1], dt,
    [&](const UAM::State& x, const double t) {
      printf("%f %f %f %f %f\n", t, x[0], x[1], x[2], x[3]);
    }
  );

  return EXIT_SUCCESS;
}