NonPeriodic Tilings With Nfold Symmetry 
A tiling consists of an arrangement of shapes covering the plane without overlap or gaps. Usually we restrict ourselves to a finite number of shapes and fixed sizes (although we can also consider tiles of geometrically decreasing and increasing sizes). The simplest and, in some sense, most natural tilings are periodic, in the sense that they consist of a metatiling of duplicated regions under translation. Three common periodic tilings, one consisting of triangles, one of squares, and one of hexagons, are illustrated below. 
These tilings have rotational periodicities as well, and since they are rotationally periodic about one point, and they are translationally periodic, it follows that they are rotationally periodic about infinitely many points, as distinct from, say, a sliced pie, which has rotational symmetry only about the central point. It's trivial to define patterns with rotational symmetry about a single point, but if we also require translational symmetry, it can be shown that the only possible rotational symmetries are of order 2, 3, 4, and 6. This is called the crystallographic restriction. 
It's also possible to tile the plane nonperiodically. The bestknown nonperiodic tilings are the Penrose tilings based on the pentagon and the golden ratio. (These tiles also have the property that they can tile the plane only nonperiodically, not periodically. Such tilings are called "aperiodic". Typically this is achieved by imposing edgematching conditions on a set of tiles to prevent them from being able to tile periodically.) It's interesting to consider other nonperiodic tilings. The figure below illustrates a tiling based on the septagon (7sided regular polygon), although it actually produces a 14way symmetry (at least about the central point). 
The 14sided polygon is tiled with just the three types of tiles shown below: 
The key to this tiling is the fact that these three basic tiles can be combined to produce scaledup copies of themselves. The scale factor is a/b = 2.24697... Alternatively, these tiles can be subdivided into scaleddown versions of themselves. Consequently, beginning with the 14sided tiled polygon shown above, we can "inflate" the figure by a factor of 2.24697... This is done by scaling up the entire figure by this factor, and then subdividing each expanded tile into tiles of the original size. Obviously we can continue this process indefinitely, so this proves that we can tile the entire plane with this 14way symmetry using just these three tiles. The result of the first inflation step is shown below. 
We will call the tiles with edge lengths (a,a,b) the "G" tiles (for Green), and the tiles with edge lengths (a,b,c) will be called "B" tiles (for Blue), and the tiles with edge lengths (a,a,c) will be called Y tiles (for Yellow). We can start with an initial 14gon tiling consisting of just 14 G tiles, but since the overall figure has perfect 14way symmetry, we need consider only one of these slices, consisting of just a single G tile. Thereafter, with each inflation step, each G tile is partitioned into 2 G tiles, 2 B tiles, and 1 Y tile. Likewise, each B tile is partitioned into 1 G tile, 2 B tiles, and 1 Y tile. Also, each Y tile is partitioned into 2 G tiles, 3 B tiles, and 2 Y tiles. Thus if we let G_{n}, B_{n}, and Y_{n} denote the numbers of tiles of the respective types contained in the original slice after the nth inflation step, we have the relations 
The determinant of the coefficient array is 1, and the characteristic polynomial is 
which has the roots 
We also have the recurrence relation for the number of G tiles after the nth inflation 
and likewise for the numbers of B and Y tiles. We can also give the number of the different kinds of tiles explicitly in terms of the characteristic roots as follows. 
where 
Interestingly, the coefficients of the expressions for B_{n} and Y_{n} are just permutations of each other (with opposite signs). The values of the characteristic roots are l_{1} = 0.307978..., l_{2} = 0.643104..., and l_{3} = 5.048917..., so as n increases the terms involving l_{1} and l_{2} go to zero and only the term involving l_{3} is significant. As a result, the numbers of tiles of the three types are asymptotic to 
and the total number of tiles goes to 
From this we can determine the asymptotic proportions of the three kinds of tiles. We have 
Notice that B_{n}/T_{n} equals the ratio of edge lengths b/a, which equals 2sin(p/14). We can also see that G_{n}/T_{n} equals 1 + a/b = l_{1}, and Y_{n}/T_{n} equals a/b  2 = a/c  1. In addition, we have Y_{n}/G_{n} = c/a. Since the characteristic polynomial is irreducible, we know the roots are irrational, so the asymptotic proportions of the tiles are also irrational, which implies that the tiling is nonperiodic. (For a periodic tiling, the entire plane is covered with copies of some finite region, in which the proportions of different tile types is obviously rational, and this is the asymptotic proportion as well.) 
Incidentally, although the characteristic polynomial (1) is irreducible over the rationals, if we put l = s^{2} we have the factorization 
In terms of the edge lengths a,b,c of our tiles, we notice that the roots of the right hand factor are c/a, a/b, and b/a  1, and the roots of the left hand factor are a/b, c/a, and c/a + 1. is a root of the left hand factor, and c/a is a root of the right hand factor. Thus it's not surprising that there are numerous relations between the edge lengths and the characteristic roots, such as 
We also have several relations between the edge lengths themselves, such as 
The left hand relation implies that c/a is the sum of the infinite geometric series in b/a, whereas the rightmost relation signifies that b is half the harmonic mean of a and c. 
The tiling is selfsimilar at each inflation stage, in the sense that the pattern on the existing region is replicated once the tiling has been expanded and the larger tiles have been partitioned back into tiles of the original size. Therefore, it is possible in principle to construct the entire tiling incrementally by adding tiles of a fixed size to the outer edge, one at a time. However, the correct pattern that must be followed in order to continuing replicating the tiling on larger and larger scales becomes increasingly complex. 
One notable feature of this tiling is that the vertices of the adjoining tiles don't appear to always be aligned. In other words, it appears that the vertex of one tile contacts another tile at some intermediate point on its edge. However, the distinction between edge points and vertices is somewhat conventional, because any point on an edge can be regarded as a vertex between two edges making an angle of p. 
The relations exhibited by this nonperiodic tiling have analogs for tilings with other rotational symmetries. In general, to find a nonperiodic tiling with nway symmetry, we need an irreducible polynomial P(x) of degree (n1)/2 with integer coefficients such that P(x^{2}) factors into two polynomials of degree (n1)/2. This is motivated by the fact that we need to be able to combine areas as well as edge lengths rationally, so we need polynomials whose roots are the squares of the roots of another related polynomial. A suitable family of polynomials is given by diagonal slices through Pascal's triangle 

To illustrate how this enables us to find nonperiodic tilings with nfold symmetry for a given n, let's take the simple case n = 5. Rather than just conducting an exhaustive search of all the possibilities for a tiling with 5fold symmetry, it's illuminating to follow a definite procedure which, according to our method, is based on the polynomial of degree (51)/2 = 2 with coefficients from a diagonal of Pascal's triangle, namely, x^{2}  3x + 1, which factors as (u^{2}  u  1)(u^{2} + u  1) where x = u^{2}. This will lead to two triangles that can be "inflated" rationally. The polynomial has the roots 2.618... and 0.381... and, as in the case of the 7fold tiling, the inflation factor will be the square root of the largest root, which gives 1.618..., i.e., the golden ratio f. In addition, the square root is 0.618..., i.e., 1/f, which will equal a ratio of edge lengths of the two tiles. 
By construction we know two copies of the larger (green) tile combined with one copy of the smaller (red) tile produce an inflated copy of the larger tile. In fact, this can be done in two distinct ways, as shown in the figure above. By inflation we can expand this pattern indefinitely. 
Incidentally, if we carry out this procedure for a 3fold symmetric tiling, we would be partitioning a slice of a hexagon, i.e., and equilateral triangle, into one part. The characteristic polynomial is x  1, so the inflation scale factor is just 1. This is a degenerate case. 
where x = u^{2}. From this we immediately know our inflation scale factor is where l is the largest characteristic value. We seek (71)/2 = 3 triangular tiles that partition a slice of a 14gon, one of which is similar to this slice, such that some combination of these parts produces inflated versions of the same parts. The area equations are lA_{1} = n_{11} A_{1} + n_{12} A_{2} + n_{13} A_{3}, and similarly for the other two areas. The characteristic equation of the coefficient array for this set of equations must be our Pascal polynomial, so we have the trace condition n_{11} + n_{22} + n_{33} = 6 in positive integers. The only possibilities are 1+1+4, 1+2+3, and 2+2+2. The most uniform is 2+2+2, so we take this option, and then the conditions on the remaining coefficients are 
For the second condition we can exclude the partitions 1+7 and 2+5, because they would violate the first condition. Thus the only possibility is 3+4, which implies one of the coefficients is 3, so the first partition of 7 must be 2+2+3. Setting n_{12}=3, we must have n_{21} = n_{31} = n_{23} = 1 and n_{13} = n_{32} = 2, so our area equation is 

Incidentally, such solutions are not unique. As an example, I've just learned that Ludwig Danzer is credited with a nonperiodic tiling with 7fold symmetry based on the same tile shapes as for the 7fold tiling described above, but they are arranged differently. It appears that the condition of "vertex to vertex" mating of the edges was imposed as an additional constraint. In any case, the area equation of the Danzer tiling is 

Notice that if we replace u with u + 1 in the first factor and with u  1 in the second factor we get the two factors of our previous characteristic polynomial. The opposite substitutions give u^{3} ± 7u^{2} + 14u ± 7. 
In general we can find a nonperiodic tiling corresponding to any given polynomial. This is particularly easy if we use rectangular tiles. To illustrate, consider the polynomial f(x) = x^{3}  x^{2}  x  1, and let s = 1.839286755... denote the largest real root. Thus we have the relations s^{3} = s^{2} + s + 1. We can now define a rectangular grid as shown on the left below, based on the powers of s. 
There are six distinct tile shapes, denoted by A through F, and the inflated versions of each tile (scaled up by a linear factor of s) can be composed of copies of the original tiles as indicated by the right hand figures above. The column and row arrangements are chosen so that no tile shares a common edge with another tile of the same type. Letting each tile symbol denote the area of the tile, the area equation is 
so the characteristic polynomial is 
Setting l = s^{2}, the first factor on the right splits into the two cubics 
consistent with the fact that the linear scale factor s is a root of x^{3} = x^{2} + x + 1. 
Clearly it's more challenging to find nonrectangular tilings. We saw previously how to determine tilings with nfold rotational symmetry for n=5 or n=7. These were both based on slices of a 2ngon, placing a tile similar to the overall slice at the end of the slice as well as at the hub. The results of carrying out this same arrangement for values of n from 3 to 8 are shown below. 
The cases of n = 3 (trivially), n=5, and n=7 exhibit special resonances that lead to successful tilings under this scheme. 
We can go on to consider the case of 9fold symmetry, based on the characteristic quartic x^{4}  10x^{3} + 15x^{2}  7x + 1. One area equation that maps to this quartic is 

and the roots of the second "u" cubic are the negations of these. We can then select a trace for the area matrix such as 7+1+1 and then determine the remaining coefficients of the matrix, leading to the area equation 

A set of tiles that satisfies these conditions is illustrated below. This figure shows four slices of a regular 18sided polygon. The first slice shows one way of partitioning such a slice into copies of three smaller shapes. The three shapes are colored yellow, green, and red, and their individual areas A_{1}, A_{2}, and A_{3} comprise an eigenvector of the area equation. The long edge of the overall slice is s times the short edge, so the inflated version of the yellow tile, labeled as Y', is simply an entire slice of the 18gon, and we can verify that it consists of 7 yellow tiles, 1 green tile, and 4 red tiles. The inflated version of the green tile is labeled as G', and is composed of 5 yellow, 1 green, and 3 red tiles. Lastly, the inflated red tile, labeled as R', is composed of 1 yellow and one red tile. 
Thus we can create an infinite nonperiodic tiling with 9fold symmetry from just the three tiles shown below. 

The integers inside the tiles signify the interior vertex angles in multiples of p/9. For a more uniform set of tiles, we can partition the yellow and green tiles, and arrive at the 9fold symmetrical tiling as indicated in the figure below. 
The slice of a regular 18gon on the left shows the basic partitioning into 3 blue, 4 yellow, and 1 green tile, and the three figures on the right indicate the inflated versions of those tiles. Naturally the area equation for this nonperiodic tiling 
still has the correct characteristic polynomial, i.e., x^{3}  9x^{2} + 6x  1. Arranging 18 slices such that all tiles meet edgetoedge and vertextovertex gives the pattern shown below. 
Notice the fundamental 9step symmetry of the dominant edge loops, as emphasized by the doublecovering cycle of the type illustrated in the figure below. 
If we carry out the first inflation step, replacing each expanded tile with its composition tiles, again requiring that tiles meet at vertices only, we get the pattern indicated by oneninth of the 18gon shown below. 
Clearly the requirement for vertextovertex joining implies that the blue tiles will always appear in the sixtile cluster, so we could replace the blue tiles with a single tile in the shape of this cluster. It's also worth noting that the yellow tiles join to other yellow tiles only along it's "end" or else as slices with a common hub, never in the opposite direction. With this restriction, it cannot tile the plane periodically. The green tiles sometimes meet parallel and sometimes antiparallel, so they do form a parallelegram, but they never meet "end to end", so they too can be prevented from tiling the plane periodically. Moreover, we can always rearrange a green and yello tile so that the green tiles meet only in the parallel sense. The consolidation of basic tiles of an nfold pattern such as this into metatiles that can tile the plane only nonperiodically (not periodically) is discussed in Aperiodic Tilings. 
The full 18gon following the first inflation step is shown in the figure below. 
All the tiles meet edgetoedge and vertextovertex at this stage (with suitable choices in the construction), but notice that this tiling involves some occurrences of two green tiles adjoined to each other "head to toe". At the next inflation and partition stage this will lead to some vertices lying on the edges of other tiles. It's possible to reverse one of the green tiles, but then the inconsistency appears at the neighboring boundary between the yellow and blue tiles. Therefore, at the next stage we necessarily have some vertextoedge meetings, as shown in the figure below. We also see the first occurrence of two yellow tiles meeting head to toe, and of a yellow tile with it's apex at the center of a face of a blue cluster. 
The red lines midway along the sides highlight the borders that do not meet vertextovertex. It's interesting that the region in the lower right corner outlined in red is identical to the original overall slice on the previous stage, although the corresponding region on the left side of this figure is not the same, which implies that the pattern of growth is not cumulative. The tiling needs to begin differently in order to cover larger and larger areas. This raises the interesting question of whether there exists a set of tiles that can tile the plane only asymptotically, in the sense that any given sequence of decisions covers only a finite region, but there is no limit to the size of the finite regions that can be covered. This is conceivable, because for a tiling based on inflation and partition, it's possible that on each stage the pattern in the original core region changes (as in the above example). The beginnings don't "converge" on any stable arrangement for any sequence of arrangements of arbitrarily increasing size. 