diff --git a/README.md b/README.md
index b16b571..b99ab66 100644
--- a/README.md
+++ b/README.md
@@ -44,7 +44,7 @@ A little of context for the follow videos:
 
 ### Before
 
-[before video](before.mp4)
+[before video](before.webm)
 
 ### After
 
diff --git a/before.webm b/before.webm
new file mode 100644
index 0000000..5926cc3
Binary files /dev/null and b/before.webm differ
diff --git a/game.h b/game.h
index d41456d..4c58c0c 100644
--- a/game.h
+++ b/game.h
@@ -13,8 +13,8 @@ extern double delta; /* simulation time delta in seconds */
 extern int width;  /* game canvas width */
 extern int height; /* game canvas height */
 
-#define DEFAULT_WIDTH 800
-#define DEFAULT_HEIGHT 800
+#define DEFAULT_WIDTH 1000
+#define DEFAULT_HEIGHT 1000
 
 extern GtkWidget* canvas; /* game canvas object */
 
diff --git a/polygons.cc b/polygons.cc
index 74f0f16..241926b 100644
--- a/polygons.cc
+++ b/polygons.cc
@@ -14,7 +14,7 @@
 #include <utility>
 
 #define ARROW_VAL_RATIO 1.7
-
+#define SUBSETPS 5
 
 bool draw_speed = true;
 
@@ -35,6 +35,21 @@ static double random_color_component() {
     return 1.0 * (rand() % 200 + 56) / 255;
 };
 
+static vec2d random_velocity() {
+    double r2;
+    vec2d v;
+    do {
+        v.x = v_min + rand() % (v_max + 1 - v_min);
+        v.y = v_min + rand() % (v_max + 1 - v_min);
+        r2 = vec2d::dot(v, v);
+    } while (r2 > v_max * v_max || r2 < v_min * v_min);
+    if (rand() % 2)
+        v.x = -v.x;
+    if (rand() % 2)
+        v.y = -v.y;
+    return v;
+}
+
 void polygons_init_state() {
     n_polygons = 20;
     polygons = new polygon[n_polygons];
@@ -68,7 +83,6 @@ void polygons_init_state() {
                         .set_center({width / 2., height - height / 17.})
                         .set_angle(45);
 
-
     // ---------- Shapes flying around start here ----------
     polygons[n++] = poly_generate::regular(100, 3)
                         .set_center({100, 400})
@@ -88,9 +102,14 @@ void polygons_init_state() {
                         .set_speed({10, 0});
 
     polygons[n++] = poly_generate::rectangle(100, 150).set_center({600, 200});
+
+    polygons[n++] = poly_generate::regular(72, 8)
+                        .set_center({500, 500})
+                        .set_speed({100, 0});
+
     polygons[n++] = poly_generate::regular(50, 5)
                         .set_center({150, 150})
-                        .set_speed({100, 0});
+                        .set_speed({150, 0});
 
     polygons[n++] = poly_generate::general({{40, 20},
                                             {40, 40},
@@ -232,19 +251,24 @@ static void check_collisions(polygon* current_p) {
 }
 
 void polygons_update_state() {
-    for (polygon* p = polygons; p != polygons + n_polygons; ++p) {
-        if (p->mass == INFINITY) // immovable objects don't need to be updated
-            continue;
-        check_collisions(p);
+    const static float sub_delta = delta / SUBSETPS;
 
-        p->rotate(delta * p->angular_speed);
-        p->angle = std::fmod(p->angle, 360);
+    for (uint i = SUBSETPS; i--;) {
+        for (polygon* p = polygons; p != polygons + n_polygons; ++p) {
+            // immovable objects don't need to be updated
+            if (p->mass == INFINITY)
+                continue;
+            check_collisions(p);
 
-        p->translate(delta * p->speed);
+            p->rotate(sub_delta * p->angular_speed);
+            p->angle = std::fmod(p->angle, 360);
 
-        vec2d g = gravity_vector(p);
-        p->translate(.5 * delta * delta * g);
-        p->speed += delta * g;
+            p->translate(sub_delta * p->speed);
+
+            vec2d g = gravity_vector(p);
+            p->translate(.5 * sub_delta * sub_delta * g);
+            p->speed += sub_delta * g;
+        }
     }
 }
 
@@ -308,7 +332,6 @@ void polygon::draw(cairo_t* cr) const {
     vec2d centroid = this->centroid();
     draw_circle(cr, centroid, 1);
 
-
     // draw speed
     if (draw_speed && this->mass != INFINITY) {
         hsv_t hsv = rgb2hsv(this->color);