<p>We present an interesting application of the solution to the simple harmonic oscillator (SHO) that can serve as a computation of <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq3.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\pi\)</EquationSource> <EquationSource Format="MATHML"><math> <mi>π</mi> </math></EquationSource> </InlineEquation>. We begin with a review of a compact teaching strategy for solving its equation of motion through integration in a general physics course, where many students face difficulties with conventional methods for solving differential equations. This integration approach leads to the arcsine function, the inverse of the sine function, ultimately providing the solution to the SHO. We investigate various series for approximating <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq3.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\pi\)</EquationSource> <EquationSource Format="MATHML"><math> <mi>π</mi> </math></EquationSource> </InlineEquation>, focusing on the arcsine series and their difference in convergence speed. We begin with Newton’s arcsine series for <InlineEquation ID="IEq5"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq5.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="96" /> </InlineMediaObject> <EquationSource Format="TEX">\(\pi =2\arcsin 1\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>π</mi> <mo>=</mo> <mn>2</mn> <mo>arcsin</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation>. We then explore a series based on powers of <InlineEquation ID="IEq6"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq6.gif" Format="GIF" Height="21" Rendition="HTML" Resolution="72" Type="Linedraw" Width="86" /> </InlineMediaObject> <EquationSource Format="TEX">\(\sin \frac{\pi }{2^{k+1}} \ll 1\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mo>sin</mo> <mfrac> <mi>π</mi> <msup> <mn>2</mn> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msup> </mfrac> <mo>≪</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation>, where <i>k</i> is a large positive integer and the sine term is computed using nested radicals through half-angle formulas, resembling Viète’s formula. The small sine term acts as a power-counting parameter, making the series better convergent to <InlineEquation ID="IEq7"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq3.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\pi\)</EquationSource> <EquationSource Format="MATHML"><math> <mi>π</mi> </math></EquationSource> </InlineEquation> with reliable error estimation. We extend this approach to a fractional-angle method, generalizing the factor from 1/2 to <InlineEquation ID="IEq8"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq8.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="31" /> </InlineMediaObject> <EquationSource Format="TEX">\(1/{p'}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>1</mn> <mo stretchy="false">/</mo> <msup> <mi>p</mi> <mo>′</mo> </msup> </mrow> </math></EquationSource> </InlineEquation> for a prime number <InlineEquation ID="IEq9"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq9.gif" Format="GIF" Height="18" Rendition="HTML" Resolution="72" Type="Linedraw" Width="16" /> </InlineMediaObject> <EquationSource Format="TEX">\(p'\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>p</mi> <mo>′</mo> </msup> </math></EquationSource> </InlineEquation>, by employing Chebyshev polynomials of the second kind, which commonly arise in physics problems. This leads to a series involving powers of <InlineEquation ID="IEq10"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq10.gif" Format="GIF" Height="22" Rendition="HTML" Resolution="72" Type="Linedraw" Width="33" /> </InlineMediaObject> <EquationSource Format="TEX">\(\sin \frac{\pi }{p}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mo>sin</mo> <mfrac> <mi>π</mi> <mi>p</mi> </mfrac> </mrow> </math></EquationSource> </InlineEquation>, where <i>p</i> is an arbitrary integer expressed as a product of prime factors, further enhancing convergence with a smaller power-counting parameter. The power counting allows us to identify significant terms in the Chebyshev polynomials and to truncate numerically insignificant contributions that optimize and simplify the computation of the sine term. Our novel strategies are pedagogical and suitable for advanced physics undergraduates, enabling them to approximate <InlineEquation ID="IEq11"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="40042_2025_1356_Article_IEq3.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\pi\)</EquationSource> <EquationSource Format="MATHML"><math> <mi>π</mi> </math></EquationSource> </InlineEquation> with high accuracy using techniques covered in physics courses.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Chebyshev polynomials boost the π-series convergence

  • Sungwoong Cho,
  • Daekyoung Kang,
  • U-Rae Kim,
  • Jungil Lee,
  • Jiawei Zhu

摘要

We present an interesting application of the solution to the simple harmonic oscillator (SHO) that can serve as a computation of \(\pi\) π . We begin with a review of a compact teaching strategy for solving its equation of motion through integration in a general physics course, where many students face difficulties with conventional methods for solving differential equations. This integration approach leads to the arcsine function, the inverse of the sine function, ultimately providing the solution to the SHO. We investigate various series for approximating \(\pi\) π , focusing on the arcsine series and their difference in convergence speed. We begin with Newton’s arcsine series for \(\pi =2\arcsin 1\) π = 2 arcsin 1 . We then explore a series based on powers of \(\sin \frac{\pi }{2^{k+1}} \ll 1\) sin π 2 k + 1 1 , where k is a large positive integer and the sine term is computed using nested radicals through half-angle formulas, resembling Viète’s formula. The small sine term acts as a power-counting parameter, making the series better convergent to \(\pi\) π with reliable error estimation. We extend this approach to a fractional-angle method, generalizing the factor from 1/2 to \(1/{p'}\) 1 / p for a prime number \(p'\) p , by employing Chebyshev polynomials of the second kind, which commonly arise in physics problems. This leads to a series involving powers of \(\sin \frac{\pi }{p}\) sin π p , where p is an arbitrary integer expressed as a product of prime factors, further enhancing convergence with a smaller power-counting parameter. The power counting allows us to identify significant terms in the Chebyshev polynomials and to truncate numerically insignificant contributions that optimize and simplify the computation of the sine term. Our novel strategies are pedagogical and suitable for advanced physics undergraduates, enabling them to approximate \(\pi\) π with high accuracy using techniques covered in physics courses.