Approximating Pi using Liu Hui's Cyclotomic Method
2020-06-25 13:13
This example demonstrates the cyclotomic method by ancient Chinese mathematician Liu Hui to approximate Pi. By increasing the sides of an inscribed regular polygon, its area approaches the circle's area. Using GeoGebra commands, it dynamically illustrates the limit concept as the area ratio converges to Pi.
<?xml version="1.0" encoding="utf-8"?>
<geogebra format="5.0" version="5.0.590.0" app="classic" platform="d" id="91842247-2de3-4620-ba9d-c90d55d6bc8b" xsi:noNamespaceSchemaLocation="http://www.geogebra.org/apps/xsd/ggb.xsd" xmlns="" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" >
<gui>
<window width="2490" height="1385" />
<perspectives>
<perspective id="tmp">
<panes>
<pane location="" divider="0.4775100401606426" orientation="1" />
</panes>
<views>
<view id="512" toolbar="0 | 1 501 5 19 , 67 | 2 15 45 18 , 7 37 | 514 3 9 , 13 44 , 47 | 16 | 551 550 11 , 20 22 21 23 , 55 56 57 , 12 | 69 | 510 511 , 512 513 | 533 531 , 534 532 , 522 523 , 537 536 , 535 | 521 520 | 36 , 38 49 560 | 571 30 29 570 31 33 | 17 | 540 40 41 42 , 27 28 35 , 6 , 502" visible="false" inframe="false" stylebar="true" location="1,1,1" size="500" window="100,100,600,400" />
<view id="2" visible="false" inframe="false" stylebar="false" location="3,3" size="475" window="1001,345,640,1141" />
<view id="4097" visible="false" inframe="true" stylebar="true" location="1,1" size="370" window="283,611,1012,946" />
<view id="1" visible="true" inframe="false" stylebar="true" location="3" size="1189" window="100,100,600,400" />
<view id="4" toolbar="0 || 2020 , 2021 , 2022 , 66 || 2001 , 2003 , 2002 , 2004 , 2005 || 2040 , 2041 , 2042 , 2044 , 2043" visible="false" inframe="false" stylebar="true" location="1" size="559" window="1210,577,643,669" />
<view id="8" toolbar="1001 | 1002 | 1003 || 1005 | 1004 || 1006 | 1007 | 1010 || 1008 1009 || 66 68 || 6" visible="false" inframe="true" stylebar="true" location="3" size="300" window="1068,712,1127,400" />
<view id="16" visible="true" inframe="false" stylebar="true" location="1" size="1270" window="1122,676,500,500" />
<view id="32" visible="false" inframe="true" stylebar="true" location="1" size="419" window="1115,597,500,785" />
<view id="64" toolbar="0" visible="false" inframe="true" stylebar="true" location="1" size="150" window="50,50,500,500" />
<view id="70" toolbar="0 || 2020 || 2021 || 2022" visible="false" inframe="true" stylebar="true" location="1" size="150" window="50,50,500,500" />
</views>
<toolbar show="true" items="0 39 73 62 | 1 501 67 , 5 19 , 72 75 76 | 2 15 45 , 18 65 , 7 37 | 4 3 8 9 , 13 44 , 58 , 47 | 16 51 64 , 70 | 10 34 53 11 , 24 20 22 , 21 23 | 55 56 57 , 12 | 36 46 , 38 49 50 , 71 14 68 | 30 29 54 32 31 33 | 25 17 26 60 52 61 | 40 41 42 , 27 28 35 , 6 | 100003 100004" position="1" help="true" />
<input show="true" cmd="true" top="false" />
<dockBar show="false" east="true" />
</perspective>
</perspectives>
<labelingStyle val="2"/>
<font size="28"/>
</gui>
<euclidianView>
<viewNumber viewNo="1"/>
<size width="1189" height="1095"/>
<coordSystem xZero="584.5577758147786" yZero="451.5939177919264" scale="40.59632169145178" yscale="40.59632169145177"/>
<evSettings axes="false" grid="false" gridIsBold="false" pointCapturing="3" rightAngleStyle="1" checkboxSize="26" gridType="3" lockedAxesRatio="1.0"/>
<bgColor r="255" g="255" b="255"/>
<axesColor r="0" g="0" b="0"/>
<gridColor r="192" g="192" b="192"/>
<lineStyle axes="1" grid="0"/>
<axis id="0" show="false" label="" unitLabel="" tickStyle="1" showNumbers="true"/>
<axis id="1" show="false" label="" unitLabel="" tickStyle="1" showNumbers="true"/>
</euclidianView>
<euclidianView>
<viewNumber viewNo="2"/>
<size width="1270" height="1095"/>
<coordSystem xZero="711.4030014471587" yZero="466.8297710337339" scale="38.116155809355725" yscale="38.1161558093553"/>
<evSettings axes="false" grid="false" gridIsBold="false" pointCapturing="3" rightAngleStyle="1" checkboxSize="26" gridType="3"/>
<bgColor r="255" g="255" b="255"/>
<axesColor r="0" g="0" b="0"/>
<gridColor r="192" g="192" b="192"/>
<lineStyle axes="1" grid="0"/>
<axis id="0" show="false" label="" unitLabel="" tickStyle="1" showNumbers="true"/>
<axis id="1" show="false" label="" unitLabel="" tickStyle="1" showNumbers="true"/>
</euclidianView>
<algebraView>
<mode val="1"/>
</algebraView>
<kernel>
<continuous val="false"/>
<usePathAndRegionParameters val="true"/>
<decimals val="5"/>
<angleUnit val="degree"/>
<algebraStyle val="0" spreadsheet="0"/>
<coordStyle val="0"/>
</kernel>
<tableview min="-2.0" max="2.0" step="1.0"/>
<scripting blocked="false" disabled="false"/>
<construction title="" author="" date="">
<command name="Circle">
<input a0="(0, 0)" a1="10"/>
<output a0="c"/>
</command>
<element type="conic" label="c">
<show object="true" label="false" ev="6"/>
<objColor r="0" g="0" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1" opacity="178"/>
<eigenvectors x0="1.0" y0="0.0" z0="1.0" x1="-0.0" y1="1.0" z1="1.0"/>
<matrix A0="1.0" A1="1.0" A2="-100.0" A3="0.0" A4="-0.0" A5="-0.0"/>
<eqnStyle style="specific"/>
</element>
<element type="numeric" label="多边形边数">
<value val="1.0"/>
<show object="true" label="false"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
<layer val="0"/>
<labelMode val="0"/>
<slider min="1" max="6" absoluteScreenLocation="true" width="200.0" x="59.0" y="51.0" fixed="false" horizontal="true" showAlgebra="true"/>
<lineStyle thickness="10" type="0" typeHidden="1"/>
<animation step="1" speed="1" type="0" playing="false"/>
</element>
<expression label="m" exp="(6 * 2^(多边形边数 - 1))" />
<element type="numeric" label="m">
<value val="6.0"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
</element>
<command name="Intersect">
<input a0="c" a1="yAxis" a2="2"/>
<output a0="A"/>
</command>
<element type="point" label="A">
<show object="false" label="false"/>
<objColor r="68" g="68" b="68" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<coords x="0.0" y="10.0" z="1.0"/>
<pointSize val="4"/>
<pointStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Rotate[A, (360° * i / m), (0, 0)]" a1="i" a2="0" a3="m"/>
<output a0="l1"/>
</command>
<element type="list" label="l1">
<show object="false" label="false" ev="6"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Polygon">
<input a0="l1"/>
<output a0="poly1"/>
</command>
<element type="polygon" label="poly1">
<lineStyle thickness="5" type="0" typeHidden="1" opacity="178"/>
<show object="true" label="false" ev="4"/>
<objColor r="0" g="102" b="153" alpha="0.10000000149011612"/>
<layer val="0"/>
<labelMode val="0"/>
</element>
<expression label="text1" exp=""正多边形边数：" + (LaTeX[m]) + """ />
<element type="text" label="text1">
<show object="true" label="false" ev="40"/>
<objColor r="0" g="0" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<isLaTeX val="true"/>
<font serif="true" sizeM="1.0" size="0" style="1"/>
<absoluteScreenLocation x="43" y="86"/>
</element>
<command name="Polygon">
<input a0="l1"/>
<output a0="poly2"/>
</command>
<element type="polygon" label="poly2">
<lineStyle thickness="3" type="0" typeHidden="1" opacity="178"/>
<show object="true" label="false" ev="11"/>
<objColor r="204" g="0" b="0" alpha="0.10000000149011612"/>
<layer val="0"/>
<labelMode val="0"/>
</element>
<command name="Sequence">
<input a0="Rotate[A, (180° * i / m), (0, 0)]" a1="i" a2="0" a3="(2 * m)"/>
<output a0="l3"/>
</command>
<element type="list" label="l3">
<show object="false" label="false" ev="7"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Segment[(0, 0), Element[l1, j]]" a1="j" a2="1" a3="m"/>
<output a0="l2"/>
</command>
<element type="list" label="l2">
<show object="true" label="false" ev="7"/>
<objColor r="204" g="0" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="3" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Segment[Element[l1, i], Element[l1, i + 1]]" a1="i" a2="1" a3="m"/>
<output a0="l5"/>
</command>
<element type="list" label="l5">
<show object="false" label="false" ev="7"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Midpoint[Element[l5, i]]" a1="i" a2="1" a3="m"/>
<output a0="l4"/>
</command>
<element type="list" label="l4">
<show object="false" label="false" ev="7"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Area">
<input a0="poly1"/>
<output a0="b"/>
</command>
<element type="numeric" label="b">
<value val="259.8076211353316"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
</element>
<expression label="text2" exp=""蓝色正多边形面积：" + (LaTeX[b]) + """ />
<element type="text" label="text2">
<show object="true" label="false" ev="40"/>
<objColor r="0" g="51" b="153" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<fixed val="true"/>
<isLaTeX val="true"/>
<font serif="true" sizeM="1.0" size="0" style="1"/>
<absoluteScreenLocation x="589" y="875"/>
</element>
<expression label="e" exp="b / 100" />
<element type="numeric" label="e">
<value val="2.598076211353316"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
</element>
<command name="Sequence">
<input a0="Vector[Element[l4, i], Element[l3, (2 * i)]]" a1="i" a2="1" a3="m"/>
<output a0="l6"/>
</command>
<element type="list" label="l6">
<show object="false" label="false" ev="7"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Translate[Element[l1, i], Vector[Element[l6, i]]]" a1="i" a2="1" a3="m"/>
<output a0="l7"/>
</command>
<element type="list" label="l7">
<show object="false" label="false" ev="11"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Translate[Element[l1, i + 1], Vector[Element[l6, i]]]" a1="i" a2="1" a3="m"/>
<output a0="l8"/>
</command>
<element type="list" label="l8">
<show object="false" label="false" ev="7"/>
<objColor r="0" g="100" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="5" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sequence">
<input a0="Polygon[Element[l1, i], Element[l1, i + 1], Element[l8, i], Element[l7, i]]" a1="i" a2="1" a3="m"/>
<output a0="l9"/>
</command>
<element type="list" label="l9">
<show object="true" label="false" ev="7"/>
<objColor r="204" g="0" b="0" alpha="0.10000000149011612"/>
<layer val="0"/>
<labelMode val="0"/>
<lineStyle thickness="3" type="0" typeHidden="1"/>
<pointSize val="5"/>
<angleStyle val="0"/>
</element>
<command name="Sum">
<input a0="l9"/>
<output a0="a"/>
</command>
<element type="numeric" label="a">
<value val="80.38475772933683"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
</element>
<expression label="d" exp="a + b" />
<element type="numeric" label="d">
<value val="340.1923788646684"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
</element>
<expression label="f" exp="d / 100" />
<element type="numeric" label="f">
<value val="3.401923788646684"/>
<objColor r="0" g="0" b="0" alpha="0.10000000149011612"/>
</element>
<expression label="text3" exp=""圆的面积范围：" + (LaTeX[b]) + "<S<" + (LaTeX[d]) + """ />
<element type="text" label="text3">
<show object="true" label="false" ev="40"/>
<objColor r="0" g="102" b="51" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<fixed val="true"/>
<isLaTeX val="true"/>
<font serif="true" sizeM="1.0" size="0" style="1"/>
<absoluteScreenLocation x="397" y="959"/>
</element>
<expression label="text4" exp=""红色区域面积：" + (LaTeX[d]) + """ />
<element type="text" label="text4">
<show object="true" label="false" ev="43"/>
<objColor r="204" g="0" b="0" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<fixed val="true"/>
<isLaTeX val="true"/>
<font serif="true" sizeM="1.0" size="0" style="1"/>
<absoluteScreenLocation x="112" y="865"/>
</element>
<expression label="text5" exp=""圆周率估计值：" + (LaTeX[e]) + "<\pi<" + (LaTeX[f]) + """ />
<element type="text" label="text5">
<show object="true" label="false" ev="43"/>
<objColor r="0" g="102" b="51" alpha="0.0"/>
<layer val="0"/>
<labelMode val="0"/>
<fixed val="true"/>
<isLaTeX val="true"/>
<font serif="true" sizeM="1.0" size="0" style="1"/>
<absoluteScreenLocation x="110" y="955"/>
</element>
</construction>
</geogebra>