View difference between Paste ID: a0SRm6W9 and 5yJJKyYt

 // método de euler para EDOs:
 // f(x+dx) = f(x) + dxF(f,x) (+ condição inicial)
 
-// particula livre 1D
+// EDO Colisão unidimensional
-// colisão elástica (v=-v) na parede em x=5;
+
 #include <stdio.h>
-#define x0 0
+
-#define v0 1
+#define dt 0.00001
-#define dt 0.01
+#define v01 30
 #define v02 10
 #define m1 5
 #define m2 8
 
 int main(void)
-  int i=0;
+
-  double x,v,t;
+  int i;
   double v1,v2,t,x1,x2;
-  // v(t+dt)=v(t) cte
+
-  // x(t+dt)=x(t)+dt*v(t)
+  x1=-10;
   x2=5;
-  v=v0;
+
   for (i=0;i<2e5;i++)
-  for (i=0;i<10000;i++)
+
       v1=v01;
       v2=v02;
-      x=x+dt*v;
+
       t=++i*dt;
-      if (x>=5)
+
       x1+=v1*dt;
-	  v=-v;
+      x2+=v2*dt;
 
       if (x1==x2 || (x1-x2)<=1)
-      printf("%lf %lf\n",x,t);
+
 	  v1=((m1-m2)/(m1+m2))*v01+ ((2*m2)/(m1+m2))*v02;
-
+	  v2=((2*m1)/(m1+m2))*v01-((m1-m2)/(m1+m2))*v02;
 	}
       printf("%lf\t %lf\t %lf\t %lf\t %lf\n",t,x1,x2,v1,v2);
     }
 
   return 0;
 }