diff --git a/week09/improvement/trees.py b/week09/improvement/trees.py
index dfa03f3..6fab445 100644
--- a/week09/improvement/trees.py
+++ b/week09/improvement/trees.py
@@ -1,17 +1,64 @@
 """This code produces a tree-like plot."""
 
+import argparse
 from math import sin, cos
 from matplotlib import pyplot as plt
-s=1
-d=[[0,1,0]]
-plt.plot([0,0],[0,1])
-for i in range(5):
-    n=[]
-    for j in range(len(d)):
-        n.append([d[j][0]+s*sin(d[j][2]-0.1), d[j][1]+s*cos(d[j][2]-0.1), d[j][2]-0.1])
-        n.append([d[j][0]+s*sin(d[j][2]+0.1), d[j][1]+s*cos(d[j][2]+0.1), d[j][2]+0.1])
-        plt.plot([d[j][0], n[-2][0]],[d[j][1], n[-2][1]])
-        plt.plot([d[j][0], n[-1][0]],[d[j][1], n[-1][1]])
-    d=n
-    s*=0.6
-plt.savefig('tree.png')
\ No newline at end of file
+
+
+# Make a class to hold each branch
+class Branch():
+    def __init__(self, x_coord, y_coord, base_angle):
+        self.x = x_coord
+        self.y = y_coord
+        self.base_angle = base_angle
+
+def make_tree(filename, num_levels = 5):
+    # Set the length of the first branch (from 0,0 to 0,trunk_length)
+    trunk_length = 1
+
+    # Set up a list holding the end of the first branch/trunk
+    # First element is the x coordinate of the centre of the tree
+    # The second element is the length of the intial trunk (the height where the branches start)
+    # Third element is angle of the branch (radians, relative to vertical)
+    tree=[Branch(0, trunk_length, 0)]
+
+    # This scales the successive branches of the tree down in length
+    scale_factor = 0.6
+
+    # This is the number of times that the tree branches after the initial trunk
+    num_branches = num_levels
+
+    # Add a plot for the "trunk" branch
+    plt.plot([0,0],[0,trunk_length])
+
+    # Make the successive branches
+    for _ in range(num_branches):
+        # Store all the branches from the current "end" here
+        branches = []
+        # Loop over each branch in the tree and add the child nodes
+        for j in range(len(tree)):
+            # For two sub branches, we go at -0.1 radians and +0.1 radians from the parent/lower branch
+            for subbranch_angle in [-0.1, 0.1]:
+                # Create the new branch and store
+                new_branch = Branch(
+                    tree[j].x + trunk_length * sin(tree[j].base_angle + subbranch_angle),
+                    tree[j].y + trunk_length * cos(tree[j].base_angle + subbranch_angle),
+                    tree[j].base_angle + subbranch_angle
+                )
+                branches.append(new_branch)
+                # Plot the branch that we just created
+                plt.plot([tree[j].x, new_branch.x], [tree[j].y, new_branch.y])
+        # Store the branches for the next iteration
+        tree=branches
+        # Reduce the trunk length for the next branches
+        trunk_length *= scale_factor
+    # Save the figure
+    plt.savefig(filename)
+
+if __name__ == "__main__":
+    parser = argparse.ArgumentParser()
+    parser.add_argument("filename")
+    parser.add_argument("--num_branches")
+    args = parser.parse_args()
+    make_tree(args.filename)
+    
\ No newline at end of file