נוסחת סטירלינג

עבור <math>x</math> גדול, <math>\ln(x!)</math> מתקרב ל-<math>x\ln(x)-x</math>
נוסחת סטירלינג (על שם המתמטיקאי הסקוטי ג'יימס סטירלינג) היא קירוב מתמטי לערך של <math>n!</math> (במילים: <math>n</math> עצרת) עבור ערכים גדולים של <math>n</math>.
<math display=block>n!\sim\sqrt{2\pi n}\left(\frac{n}{e}\right)^{n}</math>
זוהי נוסחה אסימפטוטית בשימוש בסימון אסימפטוטי, ופירושה שבגבול היחס שואף ל-1:
:<math display=block>\lim_{n\to\infty}\frac{n!}{\sqrt{2\pi n}\left(\frac{n}{e}\right)^n}=1</math>
כתוצאה מכך (כפי שיפורט להלן) <math>\ln(n!)\approx n\ln(n)-n</math>.

בגרסה כללית יותר, הנוסחה נותנת הערכה לפונקציית גמא המהווה הרחבה של פונקציית העצרת: <math>\Gamma(n+1)=n!</math>.

משפט: קיימת פונקציה ממשית <math>\eta:(0,\infty)\to\R</math> המקיימת: <math>\Gamma(x)=\sqrt{2\pi}x^{x-0.5}e^{-x}e^{\eta(x)}\qquad0<\eta(x)<\frac{1}{12x}</math>

עבור <math>n</math>-ים גדולים <math>e^{\eta(n)}\approx1</math> ולכן <math>\Gamma(n)=(n-1)!\approx\sqrt{2\pi}n^{n-1/2}e^{-n}</math>. כפל שני האגפים ב-<math>n</math> יתן את הנוסחה ל-<math>n!</math>.

פיתוח הנוסחה מתבסס על פיתוח אסימפטוטי לטור של האינטגרל המגדיר את פונקציית גמא והפיכתו לאינטגרל של גאוסיאן כפול תיקונים מסדרים שונים.

==הוכחה אסימפטוטית==
הצדקה אסימפטוטית לקירוב סטירלינג עוברת דרך ההוכחה כי <math>\ln(n!)\approx n\ln(n)-n</math>.

נפתח את הביטוי <math>\ln(n!)</math>:
:<math>\ln(n!)=\ln\left(\prod_{k\,=\,1}^nk\right)=\sum_{k\,=\,1}^n\ln(k)\approx\int\limits_1^n\ln(x)dx=n\ln(n)-n+1</math>
כאשר את המעבר מצד שמאל של המשוואה עשינו בעזרת הגדרת העצרת. את השוויון השני קיבלנו בעזרת חוקי הלוגריתמים. את הקירוב קיבלנו בעזרת סכום דארבו עליון. כל המעברים האלה נכונים מבחינה מתמטית גם עבור <math>\lim_{n\to\infty}\ln(n!)</math> ולכן <math>\ln(n!)\in\Theta(n\ln(n))</math>.

לוגריתם העצרת אינו שווה במדויק לערך האינטגרל שחושב קודם, וכדי להגיע לצורה המדויקת של נוסחת סטירלינג יש להעריך את האינטגרל בעזרת קירוב מסדר גבוה יותר. את האינטגרל ניתן לקרב בעזרת כלל הטרפז הנובע מנוסחת אוילר-מקלורן – נחלק את השטח מתחת לגרף של <math>\ln(x)</math> בתחום <math>[1,n]</math> ל-<math>n-1</math> טרפזים בעלי רוחב 1 וגבהים <math>\ln(k),\ln(k+1)</math>, וכך נוכל לקרב את האינטגרל על ידי הסכום:
:<math>\int\limits_1^n\ln(x)dx\approx\sum_{k\,=\,1}^{n-1}\frac{1}{2}(\ln(k)+\ln(k+1))=\sum_{k\,=\,1}^{n-1}\ln(k)+\frac{1}{2}\ln(n) = \ln(n!)-\frac{1}{2}\ln(n)</math>
שכן כל האיברים בסכום פרט לאיבר הראשון והאחרון מופיעים פעמיים (בעוד האיבר הראשון הוא <math>\ln(1)=0</math>). מכאן הדרך לקירוב סטירלינג ישירה:
:<math>\ln(n!)-\frac{1}{2}\ln(n)\approx\int\limits_1^n\ln(x)dx=n\ln(n)-n+1\implies\ln(n!)\approx(n+0.5)\ln(n)-n=\ln(n^{n+0.5})-\ln(e^n)=\ln\bigl(\sqrt{n}\bigl(\tfrac{n}{e}\bigr)^n\bigr)</math> 
כאשר השמטנו את הקבוע 1 שהופיע בערך המדויק של האינטגרל המסוים משום שהוא לא תורם להתנהגות האסימפטוטית של <math>n!</math>, אלא נכנס למעשה בחישוב הקבוע <math>\sqrt{2\pi}</math> (שאינו מובא כאן), כך שמקבלים בסופו של דבר:
:<math>n!\sim\sqrt{2\pi n}\left(\frac{n}{e}\right)^n</math>

==היסטוריה==
הנוסחה התגלתה לראשונה על ידי אברהם דה מואבר בצורה:
:<math display=block>n!\sim[\text{constant}]\cdot n^{n+0.5}e^{-n}</math>
דה מואבר נתן קירוב רציונלי לערך הלוגריתם הטבעי של הקבוע. תרומתו של סטירלינג הייתה להראות שהקבוע שווה במדויק ל-<math>\sqrt{2\pi}</math>.

==קישורים חיצוניים==

*
*

קטגוריה:משפטים באנליזה
קטגוריה:משפטים בתורת המספרים