From 881b2dd9cb222635b2e8e7f5ec9067ea3ae79b8e Mon Sep 17 00:00:00 2001
From: Ben Hambrecht <ben@hambrecht.ch>
Date: Wed, 11 Apr 2018 17:03:28 +0200
Subject: [PATCH] got code from EoC working (for the 2nd area model example)

---
 active_projects/eop/chapter1.py | 48 +++++++++++++++++++++++++++++++++
 1 file changed, 48 insertions(+)

diff --git a/active_projects/eop/chapter1.py b/active_projects/eop/chapter1.py
index 87656e85..7cb6c469 100644
--- a/active_projects/eop/chapter1.py
+++ b/active_projects/eop/chapter1.py
@@ -755,6 +755,54 @@ class IllustrateAreaModel1(Scene):
 
 
 
+class IllustrateAreaModel2(AreaIsDerivative):
+
+    CONFIG = {
+        "y_max" : 4,
+        "y_min" : -4,
+        "num_iterations" : 7,
+        "y_axis_label" : "",
+        "num_rects" : 400,
+        "dT" : 0.25,
+        "variable_point_label" : "T",
+        "area_opacity" : 0.8,
+    }
+    def construct(self):
+
+        self.setup_axes()
+        self.introduce_variable_area()
+
+        graph, label = self.get_v_graph_and_label()
+
+        rect_list = self.get_riemann_rectangles_list(
+            graph, self.num_iterations
+        )
+        VGroup(*rect_list).set_fill(opacity = 0.8)
+        rects = rect_list[0]
+
+        self.play(ShowCreation(graph))
+        self.play(Write(rects))
+        for new_rects in rect_list[1:]:
+            rects.align_submobjects(new_rects)
+            for every_other_rect in rects[::2]:
+                every_other_rect.set_fill(opacity = 0)
+            self.play(Transform(
+                rects, new_rects,
+                run_time = 2,
+                submobject_mode = "lagged_start"
+            ))
+        self.wait()
+
+#        self.play(FadeOut(self.x_axis.numbers))
+        self.add_T_label(6)
+        self.change_area_bounds(
+            new_t_max = 4,
+            rate_func = there_and_back,
+            run_time = 2
+        )
+
+    def func(self, x):
+        return np.exp(-x**2/2)
 
 
 class AreaSplitting(Scene):