Koch Snowflake-An Iterated Function System

by nvanoort, published

Koch Snowflake-An Iterated Function System by nvanoort Apr 5, 2016
0 Share
Download All Files

Thing Apps Enabled

Order This Printed View All Apps



Apps Used

MakerBot Desktop

Desktop App

Liked By

View All

Give a Shout Out

If you print this Thing and display it in public proudly give attribution by printing and displaying this tag.

Print Thing Tag

Thing Statistics

257Views 96Downloads


Nicole Van Oort
George Mason University
Math 493-Mathematics Through 3D Printing
Shows the progression of an iterated function system from a hexagon to a Koch Snowflake.

The Koch Snowflake also known as the Koch curve is a fractal curve, and an iterated function system. Fractals are images, or mathematical sets that display a repeating pattern at every scale, also known as self-symmetry. In other words at any scale you look at the image, the image appears to copy or repeat itself. One can also describe a fractal as a shape made up of smaller images of itself. An iterated function system is a method of creating a fractal, repetition of a function creates a union of several copies of the same fractal varying in size, rotation and translation. The Koch Snowflake is based of iterations of the Koch curve, which was first described by Helge von Koch, a Swedish mathematician in 1904. This was one of the earliest fractal curves to be described and he described it in his paper “On a continuous curve without tangents constructible from elementary geometry”. The Koch curve follows the same iteration pattern as the Koch Snowflake, but starts out as a line and then emerges into equilateral triangles, rather than starting out as a complete shape for the first iteration like the Koch Snowflake does. The box containing fractal dimension of the Koch Curve is log 4/log 3 ≈ 1.26186.
The Koch Snowflake has multiple ways it can be created, although not all of them are true iterated function systems (adding iterations in different places as the iteration continues). One example of this version is starting the Koch Snowflake with an equilateral triangle and iterating smaller equilateral triangles off of each side. The issue with the version is that after the first iteration, equilateral triangles are only iterated off of two sides of the pervious triangles to produce the accurate Koch Snowflake. Therefore the iteration doesn’t stay the same throughout and it is not a true iterated functions system. This is the most common representative version of the Koch Snowflake seen in the mathematical world. The version of the Koch Snowflake I utilized is a version using a true iterated function system. Instead of using triangles, it utilizes hexagons and iterates 6 smaller hexagons around its edges, while translating and rotating them different amounts. The iterated function system of equations used to produce this Koch Snowflake are as follows where the scalex and scaley represent the scaling factors, theta represents rotation, and e and f represent translations in the x and y directions respectively. A more detailed version of this system can be seen in the photos provided.
scale1x = 1/sqrt(3);
scale2x = 0.333333;
scale3x = 0.333333;
scale4x = 0.333333;
scale5x = 0.333333;
scale6x = 0.333333;
scale7x = 0.333333;
scale1y = 1/sqrt(3);
scale2y = 0.333333;
scale3y = 0.333333;
scale4y = 0.333333;
scale5y = 0.333333;
scale6y = 0.333333;
scale7y = 0.333333;
scalez = 2;
theta1 = 30;
theta2 = 0;
theta3 = 0;
theta4 = 0;
theta5 = 0;
theta6 = 0;
theta7 = 0;
e1 = 0;
e2 = 0.57735027;
e3 = 0;
e4 = -0.57735027;
e5 = -0.57735027;
e6 = 0;
e7 = 0.57735027;
f1 = 0;
f2 = 0.33333333;
f3 = 0.66666666;
f4 = 0.33333333;
f5 = -0.3333333;
f6 = -0.6666666;
f7 = -0.3333333;

The way one can obtain a Koch Snowflake from the basic hexagon in 3D vary, but I utilized an OpenSCAD code by Torsten Paul and modified it to fit the parameters and function qualifications for the Koch curve and the Koch Snowflake. One of the first modifications his code needed was to increase the number of smaller images produced during each iteration. His specified that each iteration would multiply the image 3 times, where I needed 7 multiplications. This fix was easy by following the same pattern he used for his three copies and I just added more so there was a total of seven. Further in detail, I needed to provide 7 x scaling factors, 7 y scaling factors, 7 rotation degrees, and 7 x and y translation amounts. Each smaller copy created during an iteration has a certain scaling degree in the x and y direction, a specific amount it is rotated as well as a certain degree of translation. The specific code can be seen in the provided images. The importance of this code is to make sure each smaller image that is iterated is placed in the accurate spot, otherwise after multiple iterations the Koch Snowflake will not be the image that gets produced. Each iteration must be in the proper location so that the correct image on the Koch Snowflake gets iterated in the end. For scaling, the scaling in the x and y direction was the same for each iterated image. The first copy was scaled by 1 over the square root of three while the reaming 6 copies were all scaled by one third. This scaling factor is placed as both the x and y factor for each of the copies. The only copy to have any rotation was the first copy, rotated 30 degrees. The remaining copies all had zero degrees of rotation so zeroes were placed in those spots of the code. Lastly was translation factors. This can be seen as the second matrix in the iterated function system for each copy. The top number is the x translation while the bottom is the y translation. These varied for each copy.
The next step in producing the proper OpenSCAD code for the Koch Snowflake was to create the proper starting polygon, making sure that the dimensions of the vertices formed an equilateral hexagon. While the original code provided the spot to place the starting polygon dimensions, it did not have a hexagon as its starting shape so therefore I needed to add vertices and change them to the dimensions of an equilateral hexagon. These new vertices also required me to change the length and width of the starting image defined at the beginning of the code. These dimensions are used to scale the overall iteration, but if I did not change them to the new dimensions of the hexagon it would force the iterations to be off centered, since length and width is used later on in the code. Therefore this was also an important step in the process. Lastly, was to create iterations for 7 copies. The code I had obtained again only produced 3 copies each iteration, when I needed seven. This was a simple fix though, copying the iteration code for the first copy I placed it 4 more times below the previous iteration code and changed the values and scaling factors to 4, 5, 6, and 7 respectively so that it would call those numbers that we defined above for that specific copy. This code can also be seen in the images attached. To give it more detail though, the code for one image copy can be seen below:
ifuns( length, width, count - 1, m mt(e1, f1) mr(theta1)*ms(scale1x,scale1y));.
This calls into account the scaling factor in both the x and y direction for the first copy (scale1x and scal1y) the rotation degree for the 1st copy (theta1) and the translation in the x (e1) and y(f1) directions for the first copy. This code was repeated 6 more times changing the number 1 to 2, 3, 4, 5, 6 and 7 respectively each time to make sure that each copy needed was going to be iterated. This left my code ready to be rendered to create this 3D Koch Snowflake.
A few issues did arise while using this code to produce the Koch Snowflake. At first, I attempted to use the equilateral triangle method not realizing that it wasn’t a true iterated function. Not only was this difficult because equations for the iterated function system weren’t provided, it was also a guessing game as to where and how far to translate and rotate each copy of the triangle. Finally after manipulating the copies into the right spots and creating multiple iterations I realized that I was getting copies that aimed back toward the middle of my first triangle and that the corresponding image didn’t appear to resemble the Koch Snowflake, and thus I decided a new method must be used. The method of using hexagons provided a method that was much easier to follow not only because it was a true iterated functions system, but that a system of equations was actually provided with it making the code that much easier to modify and produce. The only issue I had with the hexagon method was that the scaling and translating factors were in relation to a unit size of 1, in other words they were numbers under the assumption that the first hexagon was of size 1. This provided problematic when finding the vertices to define this initial hexagon. At first making the first hexagon too large, made the whole iteration all off. Thus I scaled down the original hexagon to have a length and width of 1. This provided a more accurate iteration and produced image of the Koch Snowflake, yet it was very small. This was a much easier fix though because all I had to do was create a scale command right before the module for the iteration was called into command, that way the whole iteration system would be scaled all together. This created a proper Koch Snowflake that was larger and ready to print. The last step was to render it, save it as an .stl file and import it into MakerBot to specify settings and get the final touches on it to print.
The code and equations for this project were fairly straightforward thanks to the ease of already having a code that modeled the same idea to go off of and the ability to alter to fit my iterated system function needs. The links below also provided the useful information such as the actual iterated function system of equations used to produce the Koch Snowflake, as well as helpful background information if more information is needed. The STL file, as well as the OpenSCAD file of the Koch Snowflake are found under thing files.

Print Settings

Printer Brand:








This print was initially tried on a 2X MakerBot printer, but after a malfunction was successfully created on a 5th generation MakerBot printer with a total print time of 3 hours and 36 minutes.

All Apps

Auto-magically prepare your 3D models for 3D printing. A cloud based 3D models Preparing and Healing solution for 3D Printing, MakePrintable provides features for model repairing, wall thickness...

App Info Launch App

Kiri:Moto is an integrated cloud-based slicer and tool-path generator for 3D Printing, CAM / CNC and Laser cutting. *** 3D printing mode provides model slicing and GCode output using built-in...

App Info Launch App
KiriMoto Thing App

With 3D Slash, you can edit 3d models like a stonecutter. A unique interface: as fun as a building game! The perfect tool for non-designers and children to create in 3D.

App Info Launch App

Print through a distributed network of 3D printing enthusiasts from across the US, at a fraction of the cost of the competitors. We want to change the world for the better through technology, an...

App Info Launch App

Quickly Scale, Mirror or Cut your 3D Models

App Info Launch App

Treatstock is an online platform that offers decentralized manufacturing services such as 3D printing and CNC machining for clients all over the world. We offer free and instant access to comparati...

App Info Launch App

3D print your favourite design with NinjaPrototype, a professional 3D manufacture with consistent quality and speed.

App Info Launch App