EDO Partícula livre 1D

// método de euler para EDOs:
// f(x+dx) = f(x) + dxF(f,x) (+ condição inicial)

// particula livre 1D
// colisão elástica (v=-v) na parede em x=5;

#define x0 0
#define v0 1
#define dt 0.01

#include <stdio.h>

int main(void)
{
  int i=0;
  double x,v,t;

  // v(t+dt)=v(t) cte
  // x(t+dt)=x(t)+dt*v(t)

  v=v0;

  for (i=0;i<10000;i++)
    {
      t=++i*dt;
      x=x+dt*v;

      if (x>=5)
    {
      v=-v;
    }

      printf("%lf %lf\n",x,t);
    }

  return 0;
}