Mathematics · الكتاب 4 · Bachelor Year 2

الرياضيات الجامعية — السنة 2

الرياضيات الجامعية — السنة 2 · Bachelor Year 2

6مقارنة الدوال

بدأ التحليل المقارب — وهو فنّ تعويض مقدار معقّد بمقدار بسيط مع خطأ محكوم — في مجلد السنة الأولى مع نشور تايلور. ويجعله هذا الفصل تخصّصًا قائمًا بذاته: النشور على سلالم عامة، ومقارنة المتسلسلة بالتكامل بكل قوتها المقاربة، وصيغة ستيرلينغ (مبرهَنةً كاملةً)، والدراسة المنهجية للمتتاليات المعرَّفة ضمنيًا. وهذه التقنيات هي الخبز اليومي للتحليل المقارب، وكل فصل لاحق يقدّر أي شيء — متسلسلات أو تكاملات أو احتمالات — يأكل من هذه المائدة.

6.1 علاقات المقارنة والسلالم

تعريف 6.1

قرب نقطة aa (مع aRa \in \R أو ±\pm\infty)، ومن أجل الدوال (أو المتتاليات، مع nn \to \infty): تُعرَّف f=o(g)f = o(g) وf=O(g)f = O(g) وfgf \sim g كما في مجلد السنة الأولى. وسلّم المقارنة عند aa عائلة من الدوال الموجبة، قابلة للمقارنة مثنى مثنى، ومرتَّبة ترتيبًا كليًا بالعلاقة o()o(\cdot) — والسلّم القياسي عند ++\infty هو

xα(lnx)β(α,βR),x^{\alpha} (\ln x)^{\beta} \qquad (\alpha, \beta \in \R),

مرتَّبًا معجميًا في (α,β)(\alpha, \beta)، ويُنقَّح عند الحاجة بالأسّيات eγx\eu^{\gamma x}.

تعريف 6.2 (النشر المقارب)

تقبل ff النشر المقارب

f=c1φ1+c2φ2++ckφk+o(φk)(φi+1=o(φi) في السلّم)f = c_1 \varphi_1 + c_2\varphi_2 + \dots + c_k \varphi_k + o(\varphi_k) \qquad (\varphi_{i+1} = o(\varphi_i) \text{ في السلّم})

إذا حققت البواقي المتتالية التقديرات المعروضة. وتكون المعاملات عندئذٍ وحيدة: c1=limf/φ1c_1 = \lim f/\varphi_1، ثم بالتراجع ci+1=lim(fjicjφj)/φi+1c_{i+1} = \lim\,(f - \sum_{j \leq i} c_j\varphi_j)/\varphi_{i+1}.

مثال 6.3

نشور تايلور نشورٌ مقاربة على السلّم (xa)k(x - a)^k عند aa. لكن المفهوم أوسع تمامًا: فعند ++\infty،

1xlnx=1x11lnxx=1x+lnxx2+o(lnxx2),\frac{1}{x - \ln x} = \frac1x \cdot \frac{1}{1 - \frac{\ln x}{x}} = \frac1x + \frac{\ln x}{x^2} + o\Bigl(\frac{\ln x}{x^2}\Bigr),

وهو نشر على سلّم مختلط — ولا تنطبق أي مبرهنة تايلور، بل المتسلسلة الهندسية وحساب المقادير oo فقط.

مثال 6.4 (السلّم القياسي مرتَّب فعلًا)

يحتاج الادعاء المعجمي في التعريف 6.1 سطرَ برهان لكل حالة. لنقارن xα(lnx)βx^{\alpha}(\ln x)^{\beta} وxα(lnx)βx^{\alpha'}(\ln x)^{\beta'} عند ++\infty. إذا كان α<α\alpha < \alpha': فالنسبة xαα(lnx)ββ0x^{\alpha - \alpha'}(\ln x)^{\beta - \beta'} \to 0، لأن قوة سالبة للمقدار xx تسحق أي قوة للمقدار lnx\ln x (ضع x=etx = \eu^t: e(αα)ttββ0\eu^{(\alpha - \alpha')t}\,t^{\beta - \beta'} \to 0 بنهاية “الأسّي يغلب كثير الحدود” في مجلد السنة الأولى). وإذا كان α=α\alpha = \alpha' وβ<β\beta < \beta': فالنسبة (lnx)ββ0(\ln x)^{\beta - \beta'} \to 0 مباشرةً. ومن ثم فإن الأزواج (α,β)(\alpha, \beta)، مرتَّبةً معجميًا، ترتّب السلّم بالعلاقة o()o(\cdot) — والتعويض x=etx = \eu^t هو الحيلة الصالحة لكل المقارنات المختلطة بين القوى واللوغاريتمات.

مثال 6.5 (ترتيب حديقة حيوان)

يجب أن تكون السلالم مرتَّبة؛ وهذا هو التمرين القياسي. عند ++\infty، نقارن n10n^{10} وelnnn\eu^{\sqrt{\ln n}\,\cdot\,\sqrt n} و2n2^n وnlnnn^{\ln n} بأخذ اللوغاريتمات:

10lnn    (lnn)2    nlnn    nln2,10\ln n \;\ll\; (\ln n)^2 \;\ll\; \sqrt{n\ln n} \;\ll\; n\ln 2 ,

حيث تعني anbna_n \ll b_n أن an=o(bn)a_n = o(b_n)؛ والمدخل الثاني هو ln(nlnn)\ln(n^{\ln n}). وتحافظ الأسّيات على هذه الفجوات التامة (فإذا كان lnunlnvn\ln u_n - \ln v_n \to -\infty فإن un/vn0u_n/v_n \to 0)، ومنه

n10=o(nlnn),nlnn=o(enlnn),enlnn=o(2n).n^{10} = o\bigl(n^{\ln n}\bigr), \qquad n^{\ln n} = o\bigl(\eu^{\sqrt{n\ln n}}\bigr), \qquad \eu^{\sqrt{n\ln n}} = o(2^n) .

والعبرة مزدوجة: قارِن دائمًا عبر اللوغاريتمات (بفروق اللوغاريتمات، لا بنسبها)، ولا تستنتج أبدًا unvnu_n \sim v_n من lnunlnvn\ln u_n \sim \ln v_n — فالزوج n10n^{10} وnlnnn^{\ln n} نسبة ln\ln فيه تؤول إلى \infty، لكن 2n2^n و4n4^n نسبة ln\ln فيهما تساوي 22 بالضبط وهما غير متكافئين إطلاقًا.

6.2 مقارنة المتسلسلة بالتكامل، مقاربيًا

مبرهنة 6.6

لتكن ff متصلة وموجبة ومتناقصة على [1,+)\intco{1}{+\infty}.

  1. إذا تقاربت 1f\int_1^{\infty} f، فإن البواقي تحقق

    n+1f    k>nf(k)    nf.\int_{n+1}^{\infty} f \;\leq\; \sum_{k > n} f(k) \;\leq\; \int_{n}^{\infty} f .
  2. وإذا تباعدت 1f\int_1^\infty f، فإن المجاميع الجزئية تحقق k=1nf(k)=1nf+C+o(1)\sum_{k=1}^{n} f(k) = \int_1^n f + C + o(1) من أجل ثابت CC ما: أي إن الفرق knf(k)1nf\sum_{k \leq n} f(k) - \int_1^n f يتقارب.

برهان. كان الحصر f(k+1)kk+1ff(k)f(k+1) \leq \int_k^{k+1} f \leq f(k) (بالتناقص) هو حيلة السنة الأولى؛ ويعطي جمعه على kn+1k \geq n+1 أو على knk \geq n البندَ (1). ومن أجل (2)، نضع uk=f(k)kk+1fu_k = f(k) - \int_k^{k+1} f: فبالحصر، 0ukf(k)f(k+1)0 \leq u_k \leq f(k) - f(k+1)، ومن ثم فإن المجاميع الجزئية للمتسلسلة uk\sum u_k محدودة بالمجموع التلسكوبي f(1)f(n+1)f(1)f(1) - f(n+1) \leq f(1): فتتقارب المتسلسلة. زيادةً على ذلك، تكون المتتالية (nn+1f)n\bigl(\int_n^{n+1} f\bigr)_n غير متزايدة (لأن ff متناقصة) وغير سالبة، ومن ثم فهي متقاربة. وبكتابة

k=1nf(k)1nf=k=1nuk+nn+1f,\sum_{k=1}^{n} f(k) - \int_1^n f = \sum_{k=1}^{n} u_k + \int_n^{n+1} f ,

يتقارب الطرف الأيمن حين nn \to \infty: أي إن الفرق يتقارب إلى ثابت CC، وهو البند (2).

مثال 6.7 (نشر المتسلسلة التوافقية)

من أجل f(t)=1tf(t) = \frac1t: Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1)، فنستعيد ثابت أويلر (مجلد السنة الأولى) ببرهان أنظف. وبدفعها رتبة أخرى (التمرين 6.3):

Hn=lnn+γ+12n+o(1n).H_n = \ln n + \gamma + \frac{1}{2n} + o\Bigl(\frac1n\Bigr).

وتُظهر الأعداد المكسب عند n=10n = 10: فلدينا H10=2.928968H_{10} = 2.928968\dots وln10=2.302585\ln 10 = 2.302585\dots، ومن ثم فالتقدير الخام للمقدار γ\gamma هو H10ln10=0.626383H_{10} - \ln 10 = 0.626383، بخطأ 0.0490.049؛ وطرح التصحيح 120\frac1{20} يعطي 0.5763830.576383، وهو يبعد عن γ=0.577216\gamma = 0.577216 بمقدار 8.31048.3\cdot10^{-4} فقط — وهو نفسه الحدّ التالي 112100\frac{1}{12\cdot100} في النشر، كما تبرهن مسألة نهاية الأسبوع (السؤال 8).

مثال 6.8 (تقدير خام للمقدار ln(n!)\ln(n!) بلا ستيرلينغ)

تحدّد حيلة الحصر وحدها موقع ln(n!)\ln(n!) بالفعل. فبما أن ln\ln متزايدة،

k1klnt ⁣dt    lnk    kk+1lnt ⁣dt,\int_{k-1}^{k}\ln t\,\dd t \;\leq\; \ln k \;\leq\; \int_{k}^{k+1}\ln t\,\dd t ,

وبالجمع على k=2,,nk = 2, \dots, n (مع 1nln=nlnnn+1\int_1^n\ln = n\ln n - n + 1):

nlnnn+1    ln(n!)    (n+1)ln(n+1)n.n\ln n - n + 1 \;\leq\; \ln(n!) \;\leq\; (n+1)\ln(n+1) - n .

والسياجان كلاهما nlnnn+O(lnn)n\ln n - n + O(\ln n): ومنه ln(n!)=nlnnn+O(lnn)\ln(n!) = n\ln n - n + O(\ln n)، وعلى الخصوص ln(n!)nlnn\ln(n!) \sim n\ln n. وما تضيفه صيغة ستيرلينغ هو الدرجتان التاليتان — المقدار 12lnn\frac12\ln n والثابت ln2π\ln\sqrt{2\pi} — وهما يكلّفان الحصر التلسكوبي الأدقّ في المبرهنة 6.13. ومعرفةُ أي دقة يشتريها كل أداة نصفُ صنعة التحليل المقارب.

مثال 6.9 (التلاشي يقتضي النشور)

لنحسب نهاية n2+nn\sqrt{n^2 + n} - n. فالحدّان كلاهما n\sim n، والكتابة “nn\sim n - n” لا معنى لها: إذ لا يمكن طرح المكافئات. فننشر بدل ذلك:

n2+nn=n(1+1n1)=n(12n18n2+O(1n3))=1218n+O(1n2):\sqrt{n^2 + n} - n = n\Bigl(\sqrt{1 + \tfrac1n} - 1\Bigr) = n\Bigl(\frac{1}{2n} - \frac{1}{8n^2} + O\Bigl(\frac1{n^3}\Bigr)\Bigr) = \frac12 - \frac{1}{8n} + O\Bigl(\frac1{n^2}\Bigr) :

فالنهاية 12\frac12، مع سرعة الاقتراب 18n\frac1{8n} مكافأةً. وتستحق الآلية اسمًا: فالفرق بين مقدارين كبيرين متكافئين يعيش كله في حدّيهما التاليين، ومن ثم يجب النشر إلى أول رتبة يختلف عندها الطرفان — مع حمل الباقي للتصديق على أن لا شيء آخر ينجو عند تلك الرتبة.

مثال 6.10 (مقارنة متباعدة، منفَّذة)

من أجل f(t)=1tlntf(t) = \frac{1}{t\ln t} على [2,+)\intco{2}{+\infty} (وهي متصلة وموجبة ومتناقصة): 2xf=lnlnxlnln2\int_2^x f = \ln\ln x - \ln\ln 2 \to \infty، ومن ثم حسب المبرهنة 6.6 (2)،

k=2n1klnk=lnlnn+C+o(1)\sum_{k=2}^{n}\frac{1}{k\ln k} = \ln\ln n + C + o(1)

من أجل ثابت CC ما. ودرسان. أولًا، التباعد حقيقي لكنه جليدي: فالمجموع الجزئي يتجاوز 44 أول مرة قرب nee4Cn \approx \eu^{\eu^{4 - C}}، وهو عدد فلكي. ثانيًا، أعطت الدالة الأصلية الشكل lnlnn\ln\ln n، ولم يُخمَّن: فمن أجل الحدود الرتيبة، يكون التكامل أداة الجمع القانونية، والثابت CC — مثل ثابت أويلر γ\gamma — هو ذاكرة الحدود الأولى.

6.3 صيغة ستيرلينغ

مبرهنة مساعدة 6.11 (تكاملات واليس، من جديد)

لتكن Wn=0π/2sinnt ⁣dtW_n = \int_0^{\pi/2} \sin^n t\,\dd t. عندئذٍ nWnWn1=π2nW_nW_{n-1} = \frac\pi2 من أجل n1n \geq 1، و(Wn)(W_n) متناقصة، وWnπ2nW_n \sim \sqrt{\dfrac{\pi}{2n}}.

برهان. تعطي المكاملة بالتجزئة nWn=(n1)Wn2nW_n = (n-1)W_{n-2} (مع n2n \geq 2)، ومن ثم فإن nWnWn1nW_nW_{n-1} ثابت في nn، ويساوي 1W1W0=π21 \cdot W_1 W_0 = \frac\pi2. التناقص: sinn+1sinn\sin^{n+1} \leq \sin^n على [0,π2]\intcc{0}{\frac\pi2}. والحصر، بالتفصيل: تعطي الرتابة Wn+1WnWn1W_{n+1} \leq W_n \leq W_{n-1}، وبالقسمة على Wn1>0W_{n-1} > 0،

nn+1=Wn+1Wn1WnWn11,\frac{n}{n+1} = \frac{W_{n+1}}{W_{n-1}} \leq \frac{W_n}{W_{n-1}} \leq 1 ,

والمتطابقة اليسرى من التراجع عند الدليل n+1n + 1. ويؤول الحدّان إلى 11: ومنه WnWn1W_n \sim W_{n-1}، ومن ثم

nWn2nWnWn1=π2Wnπ2n.nW_n^2 \sim nW_nW_{n-1} = \frac\pi2 \qquad\Longrightarrow\qquad W_n \sim \sqrt{\frac{\pi}{2n}} .

مثال 6.12 (تكاملات واليس الأولى)

من W0=π2W_0 = \frac\pi2 وW1=1W_1 = 1 والتراجع nWn=(n1)Wn2nW_n = (n-1)W_{n-2}:

W2=π4,W3=23,W4=3π16,W5=815,W6=5π32.W_2 = \frac\pi4, \qquad W_3 = \frac23, \qquad W_4 = \frac{3\pi}{16}, \qquad W_5 = \frac{8}{15}, \qquad W_6 = \frac{5\pi}{32}.

تحمل الأدلة الزوجية عاملًا π\pi، والفردية ناطقة — وهما الجداءان المتداخلان للصيغ المغلقة. وعدديًا W60.4909W_6 \approx 0.4909 في مقابل المقارب π/120.5116\sqrt{\pi/12} \approx 0.5116: فعند n=6n = 6 يكون المكافئ ضمن 5%5\% بالفعل، والمتطابقة الجدائية دقيقة عند كل nn: 6W6W5=65π32815=π26\,W_6W_5 = 6\cdot\frac{5\pi}{32}\cdot\frac8{15} = \frac\pi2. وجداول صغيرة كهذه أرخص طريقة لاصطياد زلة جبرية قبل أن تعدي حجة مقاربة.

مبرهنة 6.13 (ستيرلينغ)

n!    2πn(ne) ⁣n.n! \;\sim\; \sqrt{2\pi n}\, \Bigl(\frac{n}{\eu}\Bigr)^{\!n} .

برهان. الخطوة 1: n!Cn(n/e)nn! \sim C \sqrt n\, (n/\eu)^n من أجل ثابت C>0C > 0 ما. نضع

dn=ln(n!)(n+12)lnn+n.d_n = \ln(n!) - \Bigl(n + \frac12\Bigr)\ln n + n .

عندئذٍ

dndn+1=(n+12)lnn+1n1=(n+12)(1n12n2+13n3+o(n3))1=112n2+o(1n2),d_n - d_{n+1} = \Bigl(n + \frac12\Bigr) \ln\frac{n+1}{n} - 1 = \Bigl(n + \frac12\Bigr)\Bigl(\frac1n - \frac{1}{2n^2} + \frac{1}{3n^3} + o\bigl(n^{-3}\bigr)\Bigr) - 1 = \frac{1}{12n^2} + o\Bigl(\frac{1}{n^2}\Bigr),

بنشر تايلور للدالة ln(1+1n)\ln(1 + \frac1n). ومن ثم تتقارب المتسلسلة (dndn+1)\sum (d_n - d_{n+1}) بإطلاق (بالمقارنة مع n2\sum n^{-2})، فتتقارب (dn)(d_n)، ولتكن نهايتها dd؛ وبأخذ الأسّي، n!Cn(n/e)nn! \sim C\sqrt n\,(n/\eu)^n مع C=edC = \eu^{d}.

الخطوة 2: C=2πC = \sqrt{2\pi} عبر واليس. تتضافر الصيغة المغلقة W2p=(2p)!4p(p!)2π2W_{2p} = \frac{(2p)!}{4^p (p!)^2}\cdot\frac\pi2 (الآتية من التراجع، وهي حساب السنة الأولى معادًا في إطار المبرهنة المساعدة 6.11) مع الخطوة 1:

W2pC2p(2p/e)2p4p(Cp(p/e)p)2π2=2pCpπ2=πC12p.W_{2p} \sim \frac{C\sqrt{2p}\,(2p/\eu)^{2p}} {4^p\,\bigl(C\sqrt p\,(p/\eu)^p\bigr)^2}\cdot\frac{\pi}{2} = \frac{\sqrt{2p}}{C\,p}\cdot\frac{\pi}{2} = \frac{\pi}{C}\cdot\frac{1}{\sqrt{2p}} .

وبالمقارنة مع W2pπ4pW_{2p} \sim \sqrt{\frac{\pi}{4p}} (المبرهنة المساعدة 6.11): يفرض πC2p=π4p(1+o(1))\frac{\pi}{C\sqrt{2p}} = \sqrt{\frac{\pi}{4p}}\,(1 + o(1)) أن C=π4p2pπ=2πC = \pi \sqrt{\frac{4p}{2p\,\pi}} = \sqrt{2\pi}.

مثال 6.14 (المعامل الثنائي المركزي)

(2nn)=(2n)!(n!)24πn(2n/e)2n2πn(n/e)2n=4nπn:\binom{2n}{n} = \frac{(2n)!}{(n!)^2} \sim \frac{\sqrt{4\pi n}\,(2n/\eu)^{2n}}{2\pi n\,(n/\eu)^{2n}} = \frac{4^n}{\sqrt{\pi n}} :

واحتمال أن يعود مسير عشوائي متناظر إلى 00 عند الزمن 2n2n هو 1πn\sim \frac{1}{\sqrt{\pi n}} — وهو إعلان عن الفصل 22.

ملاحظة 6.15 (آفاق داخل هذا المجلد)

يتكلم كل فصل كمّي آتٍ لغةَ هذا الفصل. فالفصل الفصل 7 يصنّف المتسلسلات بمقارنة حدودها بالسلّم nα(lnn)βn^{-\alpha}(\ln n)^{-\beta} — وترسم مسألة نهاية الأسبوع فيه تلك الحدود كاملةً. ويفعل الفصل 9 الشيء نفسه من أجل التكاملات المعتلة، بالسلّم نفسه في المتغيّر المتصل. ويحسب الفصل 11 أنصاف أقطار التقارب من lim supan1/n\limsup\abs{a_n}^{1/n}، وهو تمرين في مكافئ الجذور من الرتبة nn يكون فيه ستيرلينغ المفتاحَ القياسي (n!nne\sqrt[n]{n!} \sim \frac n\eu، التمرين 6.4). وتصرف فصول الاحتمال ستيرلينغ مباشرةً: فالتقديرات المحلية في الفصل 22 للمعاملات الثنائية هي المثال 6.14 و المثال 6.21 حرفيًا. فالتحليل المقارب ليس فصلًا هنا؛ بل هو لكنة المجلد.

طريقة 6.16 (قائمة تدقيق الإقلاع المتدرّج)

قبل الوثوق بنشر متدرّج، دقّق أربع نقاط. (1) الوجود أولًا: يجب تثبيت الجذر أو المتتالية (بالرتابة والقيم الوسطى) قبل أي نشر — فالرموز بلا مرجع تنتشر بجمال ولا تعني شيئًا. (2) رتبة واحدة في كل مرور: لا يُوثق بكل تعويض إلا إلى رتبة التقدير المُدخَل؛ واستخراج حدّين جديدين من مرور واحد هو المصدر الكلاسيكي للمعاملات الخاطئة. (3) البواقي ترافقك: احمل المقدار o()o(\cdot) عبر كل خطوة جبرية ودَع الابتلاع (ابتلاع البواقي الأكبر للحدود الأصغر) يقع في النهاية، صراحةً. (4) التدقيق العددي: قيّم عند قيمة أمينة واحدة للمقدار nn؛ فخطأ المعامل ينجو من إعادة الاستنباط الجبري كثيرًا على نحو مفاجئ، ولا ينجو من الحساب أبدًا تقريبًا.

ملاحظة 6.17 (مزالق شائعة)

(1) تُجمَع المكافئات جمعًا سيئًا: فمن unn+lnnu_n \sim n + \ln n و vnnv_n \sim -n لا يجوز استنتاج un+vnlnnu_n + v_n \sim \ln n؛ فالتلاشيات تقتضي نشورًا ببواقٍ صريحة، لا مكافئات مجرّدة. (2) لا تأخذ أسّي تكافؤ أبدًا: فلدينا n+1nn + 1 \sim n لكن en+1≁en\eu^{n+1} \not\sim \eu^n؛ والاتجاه الآمن هو أخذ لوغاريتمات المكافئات المؤولة إلى ++\infty (مسألة نهاية الأسبوع في هذا الفصل، السؤال 24). (3) النشر المقارب مرتبط بسلّم: فكتابة f=1x+o(1x2)f = \frac1x + o\bigl(\frac1{x^2}\bigr) تدّعي أكثر مما تدّعيه f=1x+o(1x)f = \frac1x + o\bigl(\frac1x\bigr)، وخلط الاثنين يُبطل الجبر اللاحق. (4) في الإقلاع المتدرّج، عوّض النشر الجاري كله، بما فيه الباقي — فإسقاط مقدار o()o(\cdot) في منتصف المرور يُنتج معاملات معقولة لكنها خاطئة. (5) تحتاج مقارنة المتسلسلة بالتكامل إلى الرتابة: فهي تفشل تمامًا من أجل الحدود المتذبذبة (قارِن sinkk\sum\frac{\sin k}k، الفصل 7).

مثال 6.18 (ستيرلينغ بالأعداد)

عند n=10n = 10: تعطي الصيغة 20π(10/e)103598696\sqrt{20\pi}\,(10/\eu)^{10} \approx 3\,598\,696 في مقابل 10!=362880010! = 3\,628\,800: أي خطأ نسبي 8.31038.3\cdot10^{-3}، وهو مذهل من أجل قول “مقارب” عند n=10n = 10. وللخطأ بنية — وهي التنقيح الدقيق n!=2πn(n/e)n(1+112n+O(n2))n! = \sqrt{2\pi n}\,(n/\eu)^n\bigl(1 + \frac1{12n} + O(n^{-2})\bigr) — ويشرح تصحيحه الأول 11208.3103\frac1{120} \approx 8.3\cdot10^{-3} الفجوةَ المرصودة بدقة تقريبًا. وآلة أويلر–ماكلورين في مسألة نهاية الأسبوع هي بالضبط المصدر المنهجي لحدود التصحيح هذه.

ملاحظة 6.19 (أين يُستعمَل هذا الفصل)

المقارنة المقاربة هي نحو كل ما هو كمّي فيما بعد: محكّات التقارب وبانوراما برتران في الفصل 7، ومحكّات قابلية المكاملة في الفصل 9، وحسابات نصف قطر التقارب في الفصل 11، ومبرهنات النهاية في الفصل 22 (حيث يُشغّل ستيرلينغ تقديرات دي موافر–لابلاس). ويُصنّع مجلد السنة الثالثة الفكرةَ الواحدة التي نبرهن عليها هنا يدويًا — استخرج الحدّ الرئيسي وحُدّ الباقي — فتصير طريقة لابلاس والتقارب المهيمَن عليه.

مثال 6.20 (تكامل مقارَن بنفسه: 2x ⁣dtlnt\int_2^x \frac{\dd t}{\ln t})

تعمل صندوق أدوات المقارنة على التكاملات أيضًا. لتكن F(x)=2x ⁣dtlntF(x) = \int_2^x\frac{\dd t}{\ln t} (والدالة المكامَلة متصلة على [2,)\intco2\infty). نكامل بالتجزئة:

F(x)=[tlnt]2x+2x ⁣dt(lnt)2=xlnx+O(2x ⁣dt(lnt)2)+O(1),F(x) = \Bigl[\frac{t}{\ln t}\Bigr]_2^x + \int_2^x\frac{\dd t}{(\ln t)^2} = \frac{x}{\ln x} + O\Bigl(\int_2^x\frac{\dd t}{(\ln t)^2}\Bigr) + O(1),

والتكامل الباقي هو o(xlnx)o\bigl(\frac{x}{\ln x}\bigr): فنقسمه عند x\sqrt x، ونحدّه بالمقدارين

2x ⁣dt(lnt)2xوxx ⁣dt(lnt)2x(lnx)2=4x(lnx)2.\int_2^{\sqrt x}\frac{\dd t}{(\ln t)^2} \leq \sqrt x \qquad\text{و}\qquad \int_{\sqrt x}^{x}\frac{\dd t}{(\ln t)^2} \leq \frac{x}{(\ln\sqrt x)^2} = \frac{4x}{(\ln x)^2} .

ومنه F(x)xlnxF(x) \sim \frac{x}{\ln x}. وسيتعرّف القراء الذين لقوا مبرهنة الأعداد الأولية في مسألة نهاية الأسبوع في هذا الفصل على FF: فهو التكامل اللوغاريتمي، وهو المقدِّر الأفضل للمقدار π(x)\pi(x)، ويبيّن الحساب أنه يتفق مع xlnx\frac{x}{\ln x} إلى الرتبة الأولى.

مثال 6.21 (ستيرلينغ على معامل ثنائي غير متوازن)

يعطي روتين العواملات الثلاثة نفسه المستعمل من أجل المثال 6.14، ومن أجل (3nn)=(3n)!n!(2n)!\binom{3n}{n} = \frac{(3n)!}{n!\,(2n)!}:

(3nn)6πn(3n/e)3n2πn(n/e)n4πn(2n/e)2n=34πn(274) ⁣n.\binom{3n}{n} \sim \frac{\sqrt{6\pi n}\,(3n/\eu)^{3n}} {\sqrt{2\pi n}\,(n/\eu)^{n}\cdot\sqrt{4\pi n}\,(2n/\eu)^{2n}} = \sqrt{\frac{3}{4\pi n}}\, \Bigl(\frac{27}{4}\Bigr)^{\!n} .

والمعدل الأسّي 274=3322\frac{27}4 = \frac{3^3}{2^2} هو e3nH(1/3)\eu^{3n\,H(1/3)} بترميز الإنتروبيا في نظرية المعلومات: فالمعاملات الثنائية غير المتوازنة تنمو أبطأ تمامًا من المعامل المركزي 4n4^n لكل خطوتين — وهنا (27/4)1/31.89<2(27/4)^{1/3} \approx 1.89 < 2 لكل خطوة. وكل مقارب ثنائي في التوافقيات والاحتمال (الفصل 22) هو هذا الحساب الواحد بأوزان مختلفة.

6.4 المتتاليات المعرَّفة ضمنيًا

طريقة 6.22

لإيجاد مقاربات حلول xnx_n لمعادلة F(x,n)=0F(x, n) = 0:

  1. وطّن: برهن على وجود xnx_n ووحدانيته في فترة محدَّدة (بالرتابة ومبرهنة القيم الوسطى)، وجد سلوكه الخام (النهاية، ورتبة النموّ).
  2. أقلِع تدريجيًا: عوّض الصورة الخامة xn=(الحدّ الرئيسي)(1+εn)x_n = (\text{الحدّ الرئيسي})(1 + \varepsilon_n) في المعادلة وحُلّ من أجل الرتبة التالية للمقدار εn\varepsilon_n؛ وكرّر، فينقّح كل مرور رتبةً واحدة.

مثال 6.23

من أجل n1n \geq 1، للمعادلة tanx=x\tan x = x حلٌّ واحد بالضبط xnx_n في (nππ2,nπ+π2)\intoo{n\pi - \frac\pi2}{n\pi + \frac\pi2} (فالدالة tanxx\tan x - x متزايدة من -\infty إلى ++\infty هناك، ومشتقها tan2x0\tan^2 x \geq 0). الخام: xn=nπ+π2ynx_n = n\pi + \frac\pi2 - y_n مع yn(0,π)y_n \in \intoo{0}{\pi}؛ وبما أن xnx_n \to \infty وtanxn=xn+\tan x_n = x_n \to +\infty، فإن xnx_n تقترب من المقارب من اليسار: yn0y_n \to 0. الإقلاع المتدرّج: tanxn=cotyn=1tanyn1yn\tan x_n = \cot y_n = \frac{1}{\tan y_n} \sim \frac{1}{y_n}، وتعطي المعادلة cotyn=xnnπ\cot y_n = x_n \sim n\pi أن yn1nπy_n \sim \frac{1}{n\pi}. ومنه

xn=nπ+π21nπ+o(1n),x_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + o\Bigl(\frac1n\Bigr),

وتستمر العملية إلى أي رتبة (التمرين 6.6).

مثال 6.24 (تشغيلة ثانية للطريقة)

حُلّ x+lnx=nx + \ln x = n مقاربيًا. التوطين: xx+lnxx \mapsto x + \ln x متزايدة من -\infty إلى ++\infty على (0,+)\intoo{0}{+\infty}: فجذر وحيد xnx_n، وxnx_n \to \infty. الخام: يعطي lnxn=o(xn)\ln x_n = o(x_n) أن xnnx_n \sim n. الإقلاع المتدرّج: من xn=nlnxnx_n = n - \ln x_n ومن lnxn=lnn+o(1)\ln x_n = \ln n + o(1) (لوغاريتمات مكافئات، وكلا الطرفين \to \infty):

xn=nlnn+o(1);x_n = n - \ln n + o(1) ;

وبمرور آخر، مع lnxn=ln(nlnn+o(1))=lnnlnnn+o(lnnn)\ln x_n = \ln\bigl(n - \ln n + o(1)\bigr) = \ln n - \frac{\ln n}{n} + o\bigl(\frac{\ln n}n\bigr):

xn=nlnn+lnnn+o(lnnn).x_n = n - \ln n + \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr).

(وللتحقق عند n=100n = 100: الجذر هو x95.4415x \approx 95.4415؛ وتعطي صيغة الحدود الثلاثة 1004.6052+0.0461=95.4409100 - 4.6052 + 0.0461 = 95.4409، وصيغة الحدّين 95.394895.3948 — فيكسب كل مرور الرتبة المتوقعة.) الحلقة نفسها، ومنظر ثالث: فطريقة الطريقة 6.22 لا تعبأ بشكل المعادلة، بل بأن يعزل كل مرور المجهولَ المهيمن فحسب.

6.5 تمارين

تمرين 6.1

انشر عند ++\infty، حدّين بعد الحدّ الرئيسي:

x2+x+1,ln(x2+x)2lnx,x+sinxxlnx.\sqrt{x^2 + x + 1} , \qquad \ln(x^2 + x) - 2\ln x, \qquad \frac{x + \sin x}{x - \ln x} .
حل

حل التمرين 6.1.

x2+x+1=x1+1x+1x2=x+12+381x+o(1x)\sqrt{x^2 + x + 1} = x\sqrt{1 + \tfrac1x + \tfrac{1}{x^2}} = x + \frac12 + \frac38\cdot\frac1x + o\bigl(\frac1x\bigr) (بالنشر الثنائي: يعطي 12u18u2\frac12 u - \frac18 u^2 مع u=1x+1x2u = \frac1x + \frac{1}{x^2} أن 12x+12x218x2=12x+38x2\frac{1}{2x} + \frac{1}{2x^2} - \frac{1}{8x^2} = \frac{1}{2x} + \frac{3}{8x^2}، ثم نضرب في xx).

ln(x2+x)2lnx=ln(1+1x)=1x12x2+o(1x2)\ln(x^2 + x) - 2\ln x = \ln\bigl(1 + \tfrac1x\bigr) = \frac1x - \frac{1}{2x^2} + o\bigl(\frac{1}{x^2}\bigr).

الدالة الثالثة: ننشر كل عامل،

x+sinxxlnx=(1+sinxx)(1+lnxx+(lnx)2x2+O((lnx)3x3)).\frac{x + \sin x}{x - \ln x} = \Bigl(1 + \frac{\sin x}{x}\Bigr) \Bigl(1 + \frac{\ln x}{x} + \frac{(\ln x)^2}{x^2} + O\Bigl(\frac{(\ln x)^3}{x^3}\Bigr)\Bigr).

ونرتّب الإسهامات على السلّم عند ++\infty: lnxx1xsinxx(lnx)2x2\frac{\ln x}{x} \gg \frac{1}{x} \geq \bigl|\frac{\sin x}{x}\bigr| \gg \frac{(\ln x)^2}{x^2}. ومن ثم فالحدّان اللذان يليان الحدّ الرئيسي 11 هما lnxx\frac{\ln x}{x}، ثم حدّ التذبذب المحدود sinxx\frac{\sin x}{x}:

x+sinxxlnx=1+lnxx+sinxx+O((lnx)2x2).\frac{x + \sin x}{x - \ln x} = 1 + \frac{\ln x}{x} + \frac{\sin x}{x} + O\Bigl(\frac{(\ln x)^2}{x^2}\Bigr).

تمرين 6.2

أعطِ الطبيعة (تقارب أو تباعد) ومقاربات الحدّ الرئيسي عند التباعد للمقدار knkα\sum_{k \leq n} k^\alpha من أجل α>1\alpha > -1 وα=1\alpha = -1 وα<1\alpha < -1، عبر المبرهنة 6.6.

حل

حل التمرين 6.2.

f(t)=tαf(t) = t^\alpha (t1t \geq 1).

α>1\alpha > -1: تباعد، وحسب المبرهنة 6.6 (2)، knkα=nα+1α+1+C+o(1)\sum_{k\leq n} k^\alpha = \frac{n^{\alpha+1}}{\alpha+1} + C + o(1) إذا كان α<0\alpha < 0 (حيث تتناقص ff)؛ ومن أجل α0\alpha \geq 0 (حيث تتزايد ff) يعطي الحصر نفسه بمتراجحات معكوسة أن knkαnα+1α+1\sum_{k \leq n} k^\alpha \sim \frac{n^{\alpha + 1}}{\alpha + 1}.

α=1\alpha = -1: Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1) (المثال 6.7).

α<1\alpha < -1: تقارب، بباقٍ k>nkαnα+1(α+1)\sum_{k > n} k^\alpha \sim \frac{n^{\alpha+1}}{-(\alpha+1)} بالحصر (1) (فحدّا التكامل كلاهما مكافئ لتلك القيمة).

تمرين 6.3 ★★

برهن على أن Hn=lnn+γ+12n+o(1n)H_n = \ln n + \gamma + \frac{1}{2n} + o\bigl(\frac1n\bigr). (ادرس vn=Hnlnnγv_n = H_n - \ln n - \gamma: بيّن أن vnvn+1=12n2+O(n3)v_n - v_{n+1} = \frac{1}{2n^2} + O(n^{-3}) واجمع الذيل، بالمقارنة مع kn12k212n\sum_{k \geq n} \frac{1}{2k^2} \sim \frac{1}{2n}المبرهنة 6.6 (1).)

حل

حل التمرين 6.3.

لتكن vn=Hnlnnγ0v_n = H_n - \ln n - \gamma \to 0. عندئذٍ

vnvn+1=lnn+1n1n+1=(1n12n2)(1n1n2)+O(1n3)=12n2+O(1n3),v_n - v_{n+1} = \ln\frac{n+1}{n} - \frac{1}{n+1} = \Bigl(\frac1n - \frac{1}{2n^2}\Bigr) - \Bigl(\frac1n - \frac{1}{n^2}\Bigr) + O\Bigl(\frac{1}{n^3}\Bigr) = \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr),

باستعمال 1n+1=1n1n2+O(n3)\frac{1}{n+1} = \frac1n - \frac{1}{n^2} + O(n^{-3}). وبما أن vn0v_n \to 0، فبالتلسكوب على الذيل:

vn=kn(vkvk+1)=kn(12k2+O(k3))=12n+O(1n2),v_n = \sum_{k \geq n} (v_k - v_{k+1}) = \sum_{k\geq n} \Bigl(\frac{1}{2k^2} + O(k^{-3})\Bigr) = \frac{1}{2n} + O\Bigl(\frac{1}{n^2}\Bigr),

حسب المبرهنة 6.6 (1) مطبَّقة على t2t^{-2} (بباقٍ 1n\sim \frac1n، منصَّفًا) وعلى t3t^{-3}. ومنه Hn=lnn+γ+12n+o(1n)H_n = \ln n + \gamma + \frac{1}{2n} + o(\frac1n).

تمرين 6.4 ★★

باستعمال ستيرلينغ، جد مكافئات المقادير: (3n)!(n!)3\dfrac{(3n)!}{(n!)^3}؛ و  n!nn\;\dfrac{n!}{n^n}؛ و  n!n\;\sqrt[n]{n!} (حين ne(1+o(1))\frac n\eu(1 + o(1))، مصوغًا بدقة إلى حدّين).

حل

حل التمرين 6.4.

ستيرلينغ ثلاث مرات:

(3n)!(n!)36πn(3n/e)3n(2πn)3/2(n/e)3n=6  27n2πn12πn2πn  =327n2πn.\frac{(3n)!}{(n!)^3} \sim \frac{\sqrt{6\pi n}\,(3n/\eu)^{3n}} {(2\pi n)^{3/2}\,(n/\eu)^{3n}} = \frac{\sqrt{6}\; 27^{\,n}}{2\pi n} \cdot \frac{1}{\sqrt{2\pi n}}\cdot\sqrt{2\pi n}\; = \frac{\sqrt3\,27^n}{2\pi n} .

(وبانتباه: 6πn(2πn)3/2=6(2πn)2πnπn=32πn\frac{\sqrt{6\pi n}}{(2\pi n)^{3/2}} = \frac{\sqrt6}{(2\pi n)\sqrt{2\pi n}}\sqrt{\pi n} = \frac{\sqrt3}{2\pi n}.)

n!nn2πnen\dfrac{n!}{n^n} \sim \sqrt{2\pi n}\,\eu^{-n}.

n!n=exp(lnn!n)\sqrt[n]{n!} = \exp\bigl(\frac{\ln n!}{n}\bigr) مع lnn!=nlnnn+12ln(2πn)+o(1)\ln n! = n\ln n - n + \frac12\ln(2\pi n) + o(1):

n!n=exp(lnn1+ln(2πn)2n+o(lnnn))=ne(1+ln(2πn)2n+o(lnnn)).\sqrt[n]{n!} = \exp\Bigl(\ln n - 1 + \frac{\ln(2\pi n)}{2n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr) = \frac{n}{\eu}\Bigl(1 + \frac{\ln(2\pi n)}{2n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

تمرين 6.5 ★★

من أجل n2n \geq 2، برهن على أن للمعادلة xn+x=1x^n + x = 1 حلًا وحيدًا xn(0,1)x_n \in \intoo{0}{1}، وأن xn1x_n \to 1، وأثبت

xn=1lnnn+o(lnnn).x_n = 1 - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr).

(من xnn=1xnx_n^n = 1 - x_n: خذ اللوغاريتمات وأقلِع تدريجيًا مع xn=1εnx_n = 1 - \varepsilon_n.)

حل

حل التمرين 6.5.

g(x)=xn+x1g(x) = x^n + x - 1 متزايدة تمامًا على [0,1]\intcc{0}{1} من 1-1 إلى 11: فجذر وحيد xnx_n. وبما أن xnn=1xn(0,1)x_n^n = 1 - x_n \in \intoo{0}{1}: فلو كان xnc<1x_n \leq c < 1 على متتالية جزئية، لكان xnncn0x_n^n \leq c^n \to 0، ومنه 1xn01 - x_n \to 0: وهو تناقض مع xncx_n \leq c. ومنه xn1x_n \to 1.

نكتب xn=1εnx_n = 1 - \varepsilon_n، εn0+\varepsilon_n \to 0^+. وتُقرأ المعادلة (1εn)n=εn(1 - \varepsilon_n)^n = \varepsilon_n، أي

nln(1εn)=lnεnnεn(1+o(1))=lnεn.n\ln(1 - \varepsilon_n) = \ln \varepsilon_n \quad\Longrightarrow\quad -n\varepsilon_n\bigl(1 + o(1)\bigr) = \ln\varepsilon_n .

ومنه nεn=lnεn(1+o(1))+n\varepsilon_n = -\ln\varepsilon_n\,(1 + o(1)) \to +\infty، و بأخذ اللوغاريتمات من جديد: lnn+lnεn=ln(lnεn)+o(1)\ln n + \ln\varepsilon_n = \ln(-\ln\varepsilon_n) + o(1). وبما أن ln(lnεn)=o(ln(1/εn))\ln(-\ln \varepsilon_n) = o(\ln(1/\varepsilon_n))، يعطي هذا أن lnεnlnn\ln\varepsilon_n \sim -\ln n، ومنه εn=lnεnn(1+o(1))lnnn\varepsilon_n = \frac{-\ln\varepsilon_n}{n}(1 + o(1)) \sim \frac{\ln n}{n}:

xn=1lnnn+o(lnnn).x_n = 1 - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}{n}\Bigr) .

تمرين 6.6 ★★

ادفع المثال 6.23 رتبة أخرى:

xn=nπ+π21nπ+12n2π+o(1n2).x_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + \frac{1}{2n^2\pi} + o\Bigl(\frac{1}{n^2}\Bigr).

(اكتب cotyn=xn\cot y_n = x_n بالضبط، وانشر coty=1yy3+o(y)\cot y = \frac1y - \frac y3 + o(y) وxn=nπ(1+12n)x_n = n\pi(1 + \frac{1}{2n} - \dots)، ثم عيّن المعاملات.)

حل

حل التمرين 6.6.

العلاقة الدقيقة: cotyn=xn=nπ+π2yn\cot y_n = x_n = n\pi + \frac\pi2 - y_n، مع yn1nπy_n \sim \frac{1}{n\pi} (لأن المثال 6.23). وننشر coty=1yy3+O(y3)\cot y = \frac1y - \frac y3 + O(y^3):

1ynyn3+O(yn3)=nπ+π2yn1yn=nπ+π2+O(1n),\frac{1}{y_n} - \frac{y_n}{3} + O(y_n^3) = n\pi + \frac\pi2 - y_n \quad\Longrightarrow\quad \frac{1}{y_n} = n\pi + \frac\pi2 + O\Bigl(\frac1n\Bigr),

(فالحدّان yn-y_n وyn3-\frac{y_n}{3} هما O(1n)O(\frac1n)). ونقلب:

yn=1nπ11+12n+O(n2)=1nπ(112n+O(1n2))=1nπ12n2π+O(1n3).y_n = \frac{1}{n\pi}\cdot\frac{1}{1 + \frac{1}{2n} + O(n^{-2})} = \frac{1}{n\pi}\Bigl(1 - \frac{1}{2n} + O\Bigl(\frac{1}{n^2}\Bigr)\Bigr) = \frac{1}{n\pi} - \frac{1}{2n^2\pi} + O\Bigl(\frac{1}{n^3}\Bigr).

ومنه

xn=nπ+π2yn=nπ+π21nπ+12n2π+o(1n2).x_n = n\pi + \frac{\pi}{2} - y_n = n\pi + \frac\pi2 - \frac{1}{n\pi} + \frac{1}{2n^2\pi} + o\Bigl(\frac{1}{n^2}\Bigr).

تمرين 6.7 ★★

عيّن limn1n!k=0nk!\lim_{n\to\infty} \dfrac{1}{n!}\sum_{k=0}^{n} k! (حُدّ مجموع كل الحدود عدا الأخيرين)، واستنتج النشر المقارب knk!=n!(1+1n+O(n2))\sum_{k \leq n} k! = n!\bigl(1 + \frac1n + O(n^{-2})\bigr).

حل

حل التمرين 6.7.

نفصل الحدّين الأكبرين:

k=0nk!=n!+(n1)!+kn2k!,kn2k!(n1)(n2)!=(n1)!.\sum_{k=0}^{n} k! = n! + (n-1)! + \sum_{k \leq n-2} k! , \qquad \sum_{k\leq n-2} k! \leq (n-1)\,(n-2)! = (n-1)! .

ومنه 11n!k!1+2n1 \leq \frac{1}{n!}\sum k! \leq 1 + \frac{2}{n}: فالنهاية 11. وبالتنقيح: (n1)!n!=1n\frac{(n-1)!}{n!} = \frac1n، ويمكن شحذ الحدّ الخام kn2k!(n1)!\sum_{k \leq n-2}k! \leq (n-1)! بالطريقة نفسها: kn2k!=(n2)!(1+O(1n))=O(n!n2)\sum_{k\leq n-2} k! = (n-2)!\,(1 + O(\frac1n)) = O\bigl(\frac{n!}{n^2}\bigr). ومنه

k=0nk!=n!(1+1n+O(1n2)).\sum_{k=0}^{n} k! = n!\Bigl(1 + \frac1n + O\Bigl(\frac{1}{n^2}\Bigr)\Bigr).

تمرين 6.8 ★★★

لتكن u0>0u_0 > 0 وun+1=un+1unu_{n+1} = u_n + \dfrac{1}{u_n}. برهن على أن unu_n \to \infty، ثم على أن un2nu_n \sim \sqrt{2n} (ادرس un2u_n^2: فتزايداتها 2+un22 + u_n^{-2}؛ ثم اجمع)، ونقّح:

un=2n(1+lnn8n+o(lnnn)).u_n = \sqrt{2n}\Bigl(1 + \frac{\ln n}{8n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

(من un2=2n+k<nuk2+u02u_n^2 = 2n + \sum_{k<n} u_k^{-2} + u_0^2 ومن uk22ku_k^2 \sim 2k: يكون المجموع 12lnn\sim \frac12\ln n حسب المبرهنة 6.6.)

حل

حل التمرين 6.8.

(un)(u_n) متزايدة؛ ولو كانت محدودة لتقاربت إلى \ell مع =+1\ell = \ell + \frac1\ell: وهذا محال. ومنه unu_n \to \infty.

المربّعات: un+12=un2+2+un2u_{n+1}^2 = u_n^2 + 2 + u_n^{-2}، ومنه

un2=u02+2n+k=0n11uk2.u_n^2 = u_0^2 + 2n + \sum_{k=0}^{n-1} \frac{1}{u_k^2} .

والمجموع o(n)o(n) (فالحدود تؤول إلى 00، بتشيزارو)، ومنه un22nu_n^2 \sim 2n وun2nu_n \sim \sqrt{2n}.

التنقيح: 1uk212k\frac{1}{u_k^2} \sim \frac{1}{2k}، ومن ثم بالمقارنة (المبرهنة 6.6، أو بمكافئات المجاميع الجزئية للمتسلسلات الموجبة) k<nuk212lnn\sum_{k<n} u_k^{-2} \sim \frac12 \ln n. ومنه

un2=2n+lnn2(1+o(1))+O(1)un=2n1+lnn4n+o(lnnn)=2n(1+lnn8n+o(lnnn)).u_n^2 = 2n + \frac{\ln n}{2}\,(1 + o(1)) + O(1) \quad\Longrightarrow\quad u_n = \sqrt{2n}\sqrt{1 + \frac{\ln n}{4n} + o\Bigl(\frac{\ln n}{n}\Bigr)} = \sqrt{2n}\Bigl(1 + \frac{\ln n}{8n} + o\Bigl(\frac{\ln n}{n}\Bigr)\Bigr).

تمرين 6.9 ★★★

(مجموع ريمان بلمسة) عيّن السلوك المقارب للمقدار

Sn=k=1n1n+klnn.S_n = \sum_{k=1}^{n} \frac{1}{n + k\ln n} .

(أخرِج العامل nn: Sn=1nk(1+klnnn)1S_n = \frac1n\sum_k \bigl(1 + \frac{k\ln n}{n}\bigr)^{-1}؛ وتعرّف على مجموع من نمط ريمان بوسيط بطيء التغيّر t=lnnt = \ln n، واحسب 01 ⁣du1+tu=ln(1+t)t\int_0^1 \frac{\dd u}{1 + tu} = \frac{\ln(1+t)}{t}، واختم بأن SnlnlnnlnnS_n \sim \frac{\ln\ln n}{\ln n}.)

حل

حل التمرين 6.9.

نُخرِج العامل nn ونضع t=lnnt = \ln n:

Sn=1nk=1n11+tkn.S_n = \frac1n \sum_{k=1}^{n} \frac{1}{1 + t\,\frac kn} .

ومن أجل tt ثابت، يكون المجموع مجموعَ ريمان للدالة u11+tuu \mapsto \frac{1}{1 + tu} على [0,1]\intcc{0}{1}؛ والدالة رتيبة في uu، ومن ثم يُحصر مجموع ريمان بالتكامل مزاحًا بخطوة واحدة:

01 ⁣du1+tu1nSn01 ⁣du1+tu+1n\int_0^1 \frac{\dd u}{1 + tu} - \frac1n \leq S_n \leq \int_0^1 \frac{\dd u}{1 + tu} + \frac1n

(وهي مقارنة مجاميع ريمان لدالة رتيبة بتكاملها، وهي صالحة من أجل كل nn بالعدد t=lnnt = \ln n الخاص به). والآن 01 ⁣du1+tu=ln(1+t)t\int_0^1 \frac{\dd u}{1 + tu} = \frac{\ln(1 + t)}{t}، و 1n=o(lntt)\frac1n = o\bigl(\frac{\ln t}{t}\bigr): ومنه

Sn=ln(1+lnn)lnn+O(1n)    lnlnnlnn.S_n = \frac{\ln(1 + \ln n)}{\ln n} + O\Bigl(\frac 1n\Bigr) \;\sim\; \frac{\ln\ln n}{\ln n} .

تمرين 6.10

برهن على المتطابقة (lnn)lnn=nlnlnn(\ln n)^{\ln n} = n^{\ln\ln n}، ثم رتّب ما يلي ترتيبًا تصاعديًا بالعلاقة o()o(\cdot) عند اللانهاية، مع البراهين: n2n^2، (lnn)lnn(\ln n)^{\ln n}، 2n2^n، n!n!، nnn^n.

حل

حل التمرين 6.10.

المتطابقة: (lnn)lnn=elnnlnlnn=(elnn)lnlnn=nlnlnn(\ln n)^{\ln n} = \eu^{\ln n\,\ln\ln n} = \bigl(\eu^{\ln n}\bigr)^{\ln\ln n} = n^{\ln\ln n}. والترتيب: نقارن اللوغاريتمات. ln(n2)=2lnn\ln(n^2) = 2\ln n؛ ln((lnn)lnn)=lnnlnlnn\ln\bigl((\ln n)^{\ln n}\bigr) = \ln n\ln\ln n؛ ln(2n)=nln2\ln(2^n) = n\ln2؛ ln(n!)=nlnnn+O(lnn)\ln(n!) = n\ln n - n + O(\ln n) (بستيرلينغ، أو بالحصر الأخشن lnn!nlnn\ln n! \sim n\ln nln(nn)=nlnn\ln(n^n) = n\ln n. وبما أن 2lnn=o(lnnlnlnn)2\ln n = o(\ln n\ln\ln n)، lnnlnlnn=o(n)\ln n\ln\ln n = o(n)، nln2=o(nlnnn)n\ln 2 = o(n\ln n - n)، وnlnnnnlnnn \ln n - n \sim n\ln n لكن n!/nn0n! / n^n \to 0 (ففرق اللوغاريتمين n+O(lnn)-n + O(\ln n) \to -\infty):

n2=o((lnn)lnn),(lnn)lnn=o(2n),2n=o(n!),n!=o(nn).n^2 = o\bigl((\ln n)^{\ln n}\bigr),\quad (\ln n)^{\ln n} = o(2^n),\quad 2^n = o(n!),\quad n! = o(n^n).

(ففي كل خطوة: يؤول فرق اللوغاريتمين إلى ++\infty، ومن ثم تؤول النسبة إلى 00.)

تمرين 6.11 ★★

(ذيل المقدار 1/k2\sum 1/k^2، بحدّين) باستعمال التلسكوب الدقيق k>n1k(k+1)=1n+1\sum_{k > n} \frac{1}{k(k+1)} = \frac{1}{n+1} والتفكيك 1k2=1k(k+1)+1k2(k+1)\frac1{k^2} = \frac{1}{k(k+1)} + \frac{1}{k^2(k+1)}، برهن على أن

k>n1k2=1n12n2+O(1n3).\sum_{k > n} \frac{1}{k^2} = \frac1n - \frac{1}{2n^2} + O\Bigl(\frac{1}{n^3}\Bigr).
حل

حل التمرين 6.11.

نفكك 1k2=1k(k+1)+1k2(k+1)\frac1{k^2} = \frac1{k(k+1)} + \frac1{k^2(k+1)} ونجمع من أجل k>nk > n:

k>n1k2=1n+1+k>n1k2(k+1),\sum_{k>n}\frac1{k^2} = \frac1{n+1} + \sum_{k>n}\frac{1}{k^2(k+1)} ,

فالمجموع الأول تلسكوبي بالضبط (1k(k+1)=1k1k+1\frac1{k(k+1)} = \frac1k - \frac1{k+1}). وأما الثاني: فلدينا 1k2(k+1)=1k3+O(1k4)\frac{1}{k^2(k+1)} = \frac1{k^3} + O\bigl(\frac1{k^4}\bigr) (لأن 1k2(k+1)1k3=1k3(k+1)\frac{1}{k^2(k+1)} - \frac1{k^3} = \frac{-1}{k^3(k+1)})، وبمقارنة التكامل k>n1k3=12n2+O(1n3)\sum_{k>n}\frac1{k^3} = \frac1{2n^2} + O\bigl(\frac1{n^3}\bigr)، k>n1k4=O(1n3)\sum_{k>n}\frac1{k^4} = O\bigl(\frac1{n^3}\bigr). ومنه

k>n1k2=1n+1+12n2+O(1n3)=1n1n2+12n2+O(1n3)=1n12n2+O(1n3),\sum_{k>n}\frac1{k^2} = \frac1{n+1} + \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr) = \frac1n - \frac1{n^2} + \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr) = \frac1n - \frac{1}{2n^2} + O\Bigl(\frac1{n^3}\Bigr),

باستعمال 1n+1=1n1n2+O(1n3)\frac1{n+1} = \frac1n - \frac1{n^2} + O\bigl(\frac1{n^3}\bigr).

تمرين 6.12 ★★★

لتكن u0=12u_0 = \frac12 وun+1=un+eunu_{n+1} = u_n + \eu^{-u_n}. برهن على أن unu_n \to \infty، ثم — بوضع vn=eunv_n = \eu^{u_n} وتبيان أن vn+1=vn+1+12vn+O(vn2)v_{n+1} = v_n + 1 + \frac{1}{2v_n} + O\bigl(v_n^{-2}\bigr) — أثبت

un=lnn+lnn2n+O(1n).u_n = \ln n + \frac{\ln n}{2n} + O\Bigl(\frac1n\Bigr).
حل

حل التمرين 6.12.

(un)(u_n) متزايدة؛ ولو كانت محدودة لتقاربت إلى قيمة منتهية \ell مع =+e\ell = \ell + \eu^{-\ell}: وهذا مستحيل. ومنه unu_n \to \infty. ولتكن vn=eunv_n = \eu^{u_n} \to \infty: عندئذٍ

vn+1=eun+eun=vne1/vn=vn(1+1vn+12vn2+O(vn3))=vn+1+12vn+O(vn2).v_{n+1} = \eu^{u_n + \eu^{-u_n}} = v_n\,\eu^{1/v_n} = v_n\Bigl(1 + \frac1{v_n} + \frac1{2v_n^2} + O\bigl(v_n^{-3}\bigr)\Bigr) = v_n + 1 + \frac{1}{2v_n} + O\bigl(v_n^{-2}\bigr).

ويعطي جمع vk+1vk=1+O(1)v_{k+1} - v_k = 1 + O(1) أولًا أن vn=n+O(n)v_n = n + O(n)، ومنه vncnv_n \geq cn ابتداءً من رتبة ما؛ وبإعادة الجمع مع 12vk=O(1k)\frac1{2v_k} = O(\frac1k) نجد vn=n+O(lnn)v_n = n + O(\ln n). وبمرور آخر: 12vk=12k(1+O(lnkk))\frac{1}{2v_k} = \frac{1}{2k}\bigl(1 + O\bigl(\tfrac{\ln k}k\bigr)\bigr)، ومنه

vn=n+k<n12k+O(1)=n+lnn2+O(1).v_n = n + \sum_{k<n}\frac1{2k} + O(1) = n + \frac{\ln n}2 + O(1).

وأخيرًا un=lnvn=lnn+ln(1+lnn2n+O(1n))=lnn+lnn2n+O(1n)u_n = \ln v_n = \ln n + \ln\Bigl(1 + \frac{\ln n}{2n} + O\bigl(\tfrac1n\bigr)\Bigr) = \ln n + \frac{\ln n}{2n} + O\bigl(\tfrac1n\bigr).

6.6 مسألة: الإقلاع المتدرّج، من أويلر–ماكلورين إلى الأعداد الأولية

نادرًا ما يسلّم مقدار ضمني أو متراكم مقارباته دفعةً واحدة؛ بل تُستخرج على مرورات، يُغذّي كل مرور فيها العلاقةَ المعرِّفة بالتقدير السابق. وتدرّب مسألة نهاية الأسبوع هذه تلك الحلقة على معادلات جديدة، وتبرهن على صيغة أويلر–ماكلورين من الرتبة الأولى (وهي ترقية شبه المنحرف لمقارنة المتسلسلة بالتكامل، بأشرطة خطأ صارمة)، وتقلب xlnx=nx\ln x = n، وتصرف أشهر شيك للطريقة: فمن مبرهنة الأعداد الأولية المقبولة، نستخرج القانون المقارب pnnlnnp_n \sim n\ln n للعدد الأولي رقم nn.

مسألة 6.1

مسألة نهاية الأسبوع — تصحيح أويلر–ماكلورين ومقاربات العدد الأولي رقم nn

الجزء الأول — حلقة الإقلاع المتدرّج على معادلة جديدة.

  1. برهن على قول الوحدانية في التعريف 6.2: إذا كان f=ikciφi+o(φk)=ikciφi+o(φk)f = \sum_{i\leq k} c_i\varphi_i + o(\varphi_k) = \sum_{i \leq k} c_i'\varphi_i + o(\varphi_k) على السلّم نفسه، فإن ci=cic_i = c_i' من أجل كل ii. ثم ادفع مثال الدرس المختلط رتبة أخرى:

    1xlnx=1x+lnxx2+(lnx)2x3+o((lnx)2x3)(x+),\frac{1}{x - \ln x} = \frac1x + \frac{\ln x}{x^2} + \frac{(\ln x)^2}{x^3} + o\Bigl(\frac{(\ln x)^2}{x^3}\Bigr) \qquad (x \to +\infty),

    واشرح لماذا لا يظهر أي حدّ cx2\frac{c}{x^2}.

  2. بيّن أنه من أجل كل n1n \geq 1 يوجد للمعادلة ex+x=n\eu^x + x = n حلٌّ حقيقي واحد بالضبط xnx_n، وأن xn+x_n \to +\infty مع xnlnnx_n \sim \ln n.
  3. أقلِع تدريجيًا مرتين:

    xn=lnnlnnn(lnn)22n2+o((lnn)2n2).x_n = \ln n - \frac{\ln n}{n} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr).
  4. تحقق عدديًا عند n=1000n = 1000: قارِن x10006.90083x_{1000} \approx 6.90083 بالقيم ذات الحدّ الواحد والحدّين والثلاثة في السؤال 3، إلى خمسة أرقام عشرية.

الجزء الثاني — أويلر–ماكلورين، من الرتبة الأولى.

  1. برهن على متطابقة نواة شبه المنحرف: من أجل gg من الصنف C2C^2 على [0,1]\intcc{0}{1}،

    01g(t) ⁣dt=g(0)+g(1)21201t(1t)g(t) ⁣dt\int_0^1 g(t)\,\dd t = \frac{g(0) + g(1)}{2} - \frac12\int_0^1 t(1 - t)\,g''(t)\,\dd t

    (كامِل 12t(1t)g\frac12 t(1-t)g'' بالتجزئة مرتين).

  2. لتكن ff من الصنف C2C^2 على [1,+)\intco{1}{+\infty} مع 1f<\int_1^\infty \abs{f''} < \infty. بيّن أن

    En=k=1nf(k)1nff(1)+f(n)2E_n = \sum_{k=1}^{n} f(k) - \int_1^n f - \frac{f(1) + f(n)}{2}

    تتقارب إلى ثابت EE، مع حدّ الذيل EEn18nf\abs{E - E_n} \leq \frac18\int_n^\infty\abs{f''}: وهي صيغة أويلر–ماكلورين من الرتبة الأولى.

  3. طبّق هذا على f(t)=1tf(t) = \frac1t: برهن على أن

    Hn=lnn+γ+12n+εn,εn18n2,H_n = \ln n + \gamma + \frac{1}{2n} + \varepsilon_n, \qquad \abs{\varepsilon_n} \leq \frac{1}{8n^2},

    فتقوّي التمرين 6.3 (وعيّن الثابت بالمقدار γ\gamma بالمقارنة مع المثال 6.7).

  4. استخرج المعامل التالي: بيّن أن εn=112n2+o(1n2)\varepsilon_n = -\frac{1}{12n^2} + o\bigl(\frac1{n^2}\bigr) (فتزايدات EnE_n هي 1201t(1t)f(n+t) ⁣dt=112f(n)+o(f(n))\frac12\int_0^1t(1-t)f''(n+t)\dd t = \frac1{12}f''(n) + o(f''(n))؛ واجمع الذيل مع المبرهنة 6.6).
  5. طبّق السؤال 6 على f=lnf = \ln: أعد استنباط تقارب dn=lnn!(n+12)lnn+nd_n = \ln n! - (n + \frac12)\ln n + n في ثلاثة أسطر (وهو الخطوة 1 في المبرهنة 6.13)، مع مكافأة معدل الخطأ dn=d+O(1n)d_n = d + O\bigl(\frac1n\bigr).
  6. طبّق السؤال 6 على f(t)=1tf(t) = \frac{1}{\sqrt t}: بيّن أن

    k=1n1k=2n+c+12n+O(1n3/2)\sum_{k=1}^{n}\frac1{\sqrt k} = 2\sqrt n + c + \frac{1}{2\sqrt n} + O\Bigl(\frac{1}{n^{3/2}}\Bigr)

    من أجل ثابت cc ما، وقيّم كل الحدود عند n=104n = 10^4 (والثابت هو c1.4604c \approx -1.4604).

الجزء الثالث — القلب: المعادلة xlnx=nx\ln x = n.

  1. بيّن أن للمعادلة xlnx=nx\ln x = n حلًا واحدًا بالضبط xn[1,+)x_n \in \intco{1}{+\infty} من أجل n1n \geq 1، وأن xnx_n \to \infty، وأن lnxnlnn\ln x_n \sim \ln n.
  2. استنتج القلب ذا الحدّ الواحد xnnlnnx_n \sim \dfrac{n}{\ln n}، ثم أقلِع تدريجيًا مرة أخرى:

    lnxn=lnnlnlnn+o(1),xn=nlnn(1+lnlnnlnn+o(lnlnnlnn)).\ln x_n = \ln n - \ln\ln n + o(1), \qquad x_n = \frac{n}{\ln n}\Bigl(1 + \frac{\ln\ln n}{\ln n} + o\Bigl(\frac{\ln\ln n}{\ln n}\Bigr)\Bigr).
  3. اختبر عند n=106n = 10^6: الجذر الحقيقي هو x87848x \approx 87\,848؛ قارِن بالقيمة ذات الحدّ الواحد (72382\approx 72\,382) وبالقيمة ذات الحدّين (86140\approx 86\,140)، واشرح بطء المكسب (فوسيط النشر هو lnlnnlnn\frac{\ln\ln n}{\ln n}، ولا يساوي إلا 0.19\approx 0.19 عند n=106n = 10^6).
  4. نقبل الآن مبرهنة الأعداد الأولية: يحقق عدد π(x)\pi(x) الأعداد الأولية التي لا تتجاوز x\leq x العلاقةَ π(x)xlnx\pi(x) \sim \frac{x}{\ln x} حين xx \to \infty (وهي مبرهَنة بأمانة في مجلد السنة الثالثة). وبكتابة pnp_n للعدد الأولي رقم nn، برّر π(pn)=n\pi(p_n) = n، وشغّل قلب السؤالين 11–12 لتبرهن على أن

    pnnlnn.p_n \sim n \ln n .
  5. الأرباح: (أ) بيّن أن knpkn2lnn2\sum_{k \leq n} p_k \sim \frac{n^2\ln n}{2} (قارِن klnk\sum k\ln k بالمقدار tlnt ⁣dt\int t\ln t\,\dd t)؛ (ب) احسب الفرصة التقريبية لأن يكون عدد صحيح عشوائي منتظم ذو 100100 رقمًا أوليًا (ln10100230.26\ln 10^{100} \approx 230.26: نحو واحد من كل 230230).

الجزء الرابع — تصدير الطريقة: xtanx=1x\tan x = 1.

  1. بيّن أنه من أجل كل n1n \geq 1 يوجد للمعادلة tanx=1x\tan x = \frac1x حلٌّ واحد بالضبط xnx_n في (nπ,nπ+π2)\intoo{n\pi}{\,n\pi + \frac\pi2}، وأن zn=xnnπ0+z_n = x_n - n\pi \to 0^+.
  2. حدٌّ واحد: zn1nπz_n \sim \dfrac{1}{n\pi}.
  3. بيّن أن نشر znz_n لا يحتوي أي حدّ cn2\frac{c}{n^2}: أي zn=1nπ+O(1n3)z_n = \frac1{n\pi} + O\bigl(\frac{1}{n^3}\bigr).
  4. ثلاثة حدود: باستعمال arctanu=uu33+O(u5)\arctan u = u - \frac{u^3}3 + O(u^5) و1xn=1nπzn(nπ)2+O(n3zn2)\frac1{x_n} = \frac{1}{n\pi} - \frac{z_n}{(n\pi)^2} + O(n^{-3}\cdot z_n^2)، برهن على أن

    xn=nπ+1nπ43π3n3+o(1n3).x_n = n\pi + \frac{1}{n\pi} - \frac{4}{3\pi^3 n^3} + o\Bigl(\frac{1}{n^3}\Bigr).
  5. تحقق عند n=3n = 3: الجذر الحقيقي x39.5293344x_3 \approx 9.5293344؛ قارِن القيمتين ذات الحدّ الواحد وذات الحدود الثلاثة، وقابِل بجملة واحدة بينها وبين المتتالية tanx=x\tan x = x في الدرس (المثال 6.23): أين تقع كل متتالية في نافذتها، ولماذا.

الجزء الخامس — إقلاع متدرّج حركي، وقواعد اللعبة، والتركيب.

  1. لتكن u0(0,π)u_0 \in \intoo{0}{\pi} وun+1=sinunu_{n+1} = \sin u_n. بيّن أن un0u_n \to 0 تناقصيًا، واحسب نهاية 1un+121un2\dfrac{1}{u_{n+1}^2} - \dfrac{1}{u_n^2} (انشر sin2\sin^{-2} عبر sinu=uu36+o(u3)\sin u = u - \frac{u^3}6 + o(u^3)).
  2. استنتج، عبر متوسطات تشيزارو (مجلد السنة الأولى)، الكلاسيكيةَ

    un3n.u_n \sim \sqrt{\frac{3}{n}} .
  3. (حسابات مصدَّقة) باستعمال الحدّ الصارم في السؤال 7، بيّن أن تقييم lnn+γ+12n\ln n + \gamma + \frac1{2n} عند n=106n = 10^6 يعطي H106H_{10^6} بخطأ لا يتجاوز 1.2510131.25\cdot10^{-13} — أي مجموع بمليون حدّ محسوبًا إلى ثلاثة عشر رقمًا بثلاثة حدود.
  4. (قواعد اللعبة) برهن أو ادحض، بالبراهين أو بالأمثلة المضادة: (أ) إذا كان unvn+u_n \sim v_n \to +\infty فإن lnunlnvn\ln u_n \sim \ln v_n؛ (ب) إذا كان unvnu_n \sim v_n فإن eunevn\eu^{u_n} \sim \eu^{v_n}؛ (ج) إذا كان fgf \sim g عند ++\infty (مع f,gf, g قابلة للاشتقاق) فإن fgf' \sim g'.
  5. (تركيب) بجملة واحدة لكل بند: حلقة الإقلاع المتدرّج في الطريقة 6.22 كما استُعملت في الأجزاء الأول والثالث والرابع؛ وماذا يضيف تصحيح شبه المنحرف إلى المبرهنة 6.6؛ ولماذا يكون قلب xlnxx\ln x هو بالضبط الجسر من π(x)\pi(x) إلى pnp_n؛ وأي قاعدة من قواعد السؤال 24 حرست أي خطوة. وسمِّ القمتين: صيغة أويلر–ماكلورين (من الرتبة الأولى)، والقانون المقارب للعدد الأولي رقم nn.
حل

حل المسألة 6.1.

1. بطرح النشرين: i(cici)φi=o(φk)\sum_i (c_i - c_i')\varphi_i = o(\varphi_k). فإذا اختلف معامل ما، وليكن i0i_0 أولها: فبالقسمة على φi0\varphi_{i_0} واستعمال φj=o(φi0)\varphi_j = o(\varphi_{i_0}) من أجل j>i0j > i_0 نجد ci0ci0=o(1)c_{i_0} - c_{i_0}' = o(1): أي صفرًا، وهو تناقض. وأما النشر: فمع u=lnxx0u = \frac{\ln x}x \to 0،

1xlnx=1x11u=1x(1+u+u2+O(u3))=1x+lnxx2+(lnx)2x3+o((lnx)2x3).\frac{1}{x - \ln x} = \frac1x\cdot\frac{1}{1 - u} = \frac1x\bigl(1 + u + u^2 + O(u^3)\bigr) = \frac1x + \frac{\ln x}{x^2} + \frac{(\ln x)^2}{x^3} + o\Bigl(\frac{(\ln x)^2}{x^3}\Bigr).

ولا يظهر أي حدّ cx2\frac c{x^2} لأن النشر متسلسلة هندسية في u=lnxxu = \frac{\ln x}{x}: فكل حدّ يحمل من قوى lnx\ln x بقدر ما يحمل من قوى 1x\frac1x بعد الأولى؛ ودرجة السلّم 1x2\frac1{x^2} (أي معامل (lnx)0(\ln x)^0) غائبة ببساطة، بمعامل 00.

2. f(x)=ex+xf(x) = \eu^x + x متصلة ومتزايدة تمامًا، ونهايتاها -\infty و++\infty: فهي تقابل RR\R \to \R، ومن ثم فإن xn=f1(n)x_n = f^{-1}(n) موجود ووحيد، و xn+x_n \to +\infty (لأن f1f^{-1} متزايدة إلى ++\infty). ومن exn=nxn\eu^{x_n} = n - x_n: xn=ln(nxn)lnnx_n = \ln(n - x_n) \leq \ln n، ومنه xn/n0x_n/n \to 0 وxn=lnn+ln(1xn/n)=lnn+o(1)lnnx_n = \ln n + \ln(1 - x_n/n) = \ln n + o(1) \sim \ln n.

3. نكتب un=xn/nu_n = x_n/n. المرور الثاني: un=lnn+o(1)nu_n = \frac{\ln n + o(1)}{n}، ومنه

xn=lnn+ln(1un)=lnnun+O(un2)=lnnlnnn+o(lnnn).x_n = \ln n + \ln(1 - u_n) = \ln n - u_n + O(u_n^2) = \ln n - \frac{\ln n}{n} + o\Bigl(\frac{\ln n}n\Bigr).

المرور الثالث: الآن un=lnnnlnnn2+o(lnnn2)u_n = \frac{\ln n}{n} - \frac{\ln n}{n^2} + o\bigl(\frac{\ln n}{n^2}\bigr)، وln(1un)=unun22+O(un3)\ln(1 - u_n) = -u_n - \frac{u_n^2}2 + O(u_n^3):

xn=lnnlnnn+lnnn2(lnn)22n2+o((lnn)2n2)=lnnlnnn(lnn)22n2+o((lnn)2n2),x_n = \ln n - \frac{\ln n}n + \frac{\ln n}{n^2} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr) = \ln n - \frac{\ln n}{n} - \frac{(\ln n)^2}{2n^2} + o\Bigl(\frac{(\ln n)^2}{n^2}\Bigr),

مع ابتلاع الحدّ lnnn2\frac{\ln n}{n^2} في o((lnn)2n2)o\bigl(\frac{(\ln n)^2}{n^2}\bigr).

4. عند n=1000n = 1000: ln10006.90776\ln 1000 \approx 6.90776 (بخطأ 71037\cdot10^{-3})؛ وبحدّين: 6.900856.90085 (بخطأ 21052\cdot10^{-5})؛ وبثلاثة حدود: 6.900826.90082 (بخطأ دون 10510^{-5})، في مقابل x10006.90083x_{1000} \approx 6.90083. ويشتري كل مرور العامل المتوقع lnnn\frac{\ln n}{n} تقريبًا.

5. مكاملتان بالتجزئة، ابتداءً من اليمين: مع  ⁣d ⁣dt[12t(1t)]=12t\frac{\dd}{\dd t}\bigl[\tfrac12t(1-t)\bigr] = \tfrac12 - t وt(1t)t(1-t) منعدمة عند الطرفين،

1201t(1t)g(t) ⁣dt=01(12t)g(t) ⁣dt=[(12t)g]0101g=g(0)+g(1)201g.\frac12\int_0^1 t(1-t)g''(t)\dd t = -\int_0^1\Bigl(\frac12 - t\Bigr)g'(t)\dd t = -\Bigl[\Bigl(\frac12 - t\Bigr)g\Bigr]_0^1 - \int_0^1 g = \frac{g(0) + g(1)}2 - \int_0^1 g .

وبإعادة الترتيب، هذه هي المتطابقة المذكورة.

6. نحسب التزايد، ثم نطبّق السؤال 5 على g(t)=f(n+t)g(t) = f(n + t):

En+1En=f(n+1)nn+1 ⁣ff(n+1)f(n)2=f(n)+f(n+1)2nn+1 ⁣f=1201t(1t)f(n+t) ⁣dt.\begin{align*} E_{n+1} - E_n &= f(n{+}1) - \int_n^{n+1}\!f - \frac{f(n{+}1) - f(n)}2 \\ &= \frac{f(n) + f(n{+}1)}2 - \int_n^{n+1}\!f = \frac12\int_0^1 t(1-t)f''(n+t)\dd t . \end{align*}

وبما أن 0t(1t)140 \leq t(1-t) \leq \frac14: فإن En+1En18nn+1f\abs{E_{n+1} - E_n} \leq \frac18\int_n^{n+1}\abs{f''}، ويتقارب مجموعه على nn بحكم الفرض: ومنه تتقارب (En)(E_n) (فتزايداتها قابلة للجمع بإطلاق) إلى EE ما، مع

EEnknEk+1Ek18nf.\abs{E - E_n} \leq \sum_{k\geq n}\abs{E_{k+1} - E_k} \leq \frac18\int_n^\infty\abs{f''} .

7. f(t)=1tf(t) = \frac1t: f(t)=2t3f''(t) = \frac2{t^3}، 1f=1<\int_1^\infty\abs{f''} = 1 < \infty. والسؤال 6:

Hn=lnn+1+1n2+E+(EnE)=lnn+(E+12)+12n+εn,H_n = \ln n + \frac{1 + \frac1n}{2} + E + (E_n - E) = \ln n + \Bigl(E + \frac12\Bigr) + \frac1{2n} + \varepsilon_n,

مع εn=EnE18n2 ⁣dtt3=18n2\abs{\varepsilon_n} = \abs{E_n - E} \leq \frac18\int_n^\infty\frac{2\dd t}{t^3} = \frac1{8n^2}. وبالمقارنة مع Hn=lnn+γ+o(1)H_n = \ln n + \gamma + o(1) (المثال 6.7) يتعيّن E+12=γE + \frac12 = \gamma.

8. من صيغة التزايد في السؤال 6،

εn=EnE=kn1201t(1t)2 ⁣dt(k+t)3=kn(1k301t(1t) ⁣dt+O(1k4)),\varepsilon_n = E_n - E = -\sum_{k\geq n}\frac12\int_0^1 t(1-t)\,\frac{2\,\dd t}{(k+t)^3} = -\sum_{k \geq n}\Bigl(\frac1{k^3}\int_0^1t(1-t)\dd t + O\Bigl(\frac1{k^4}\Bigr)\Bigr),

باستعمال 1(k+t)3=1k3+O(1k4)\frac{1}{(k+t)^3} = \frac1{k^3} + O\bigl(\frac1{k^4}\bigr) بانتظام من أجل t[0,1]t \in \intcc01. ومع 01t(1t)=16\int_0^1 t(1-t) = \frac16 و kn1k312n2\sum_{k\geq n}\frac1{k^3} \sim \frac{1}{2n^2} (المبرهنة 6.6):

εn=1612n2+o(1n2)=112n2+o(1n2).\varepsilon_n = -\frac16\cdot\frac{1}{2n^2} + o\Bigl(\frac1{n^2}\Bigr) = -\frac{1}{12n^2} + o\Bigl(\frac{1}{n^2}\Bigr).

9. f=lnf = \ln: f(t)=1t2f''(t) = -\frac1{t^2}، وهي قابلة للمكاملة بإطلاق. ويعطي السؤال 6 أن

lnn!=1nlnt ⁣dt+lnn2+E+O(18n ⁣dtt2)=(n+12)lnnn+1+E+O(1n),\ln n! = \int_1^n\ln t\,\dd t + \frac{\ln n}2 + E + O\Bigl( \frac1{8}\int_n^\infty\frac{\dd t}{t^2}\Bigr) = \Bigl(n + \frac12\Bigr)\ln n - n + 1 + E + O\Bigl(\frac1n\Bigr),

ومنه dn=1+E+O(1n)d_n = 1 + E + O\bigl(\frac1n\bigr): أي تقارب (dn)(d_n) — وهو الخطوة 1 في المبرهنة 6.13 — مع المعدل O(1/n)O(1/n). (وتعطي قيمة النهاية عند ستيرلينغ أن E=ln2π1E = \ln\sqrt{2\pi} - 1.)

10. f(t)=t1/2f(t) = t^{-1/2}: f(t)=34t5/2f''(t) = \frac34 t^{-5/2}، وهي قابلة للمكاملة بإطلاق. والسؤال 6:

k=1n1k=2n2+1+1n2+E+O(n3/2)=2n+c+12n+O(n3/2),\sum_{k=1}^n \frac1{\sqrt k} = 2\sqrt n - 2 + \frac{1 + \frac1{\sqrt n}}2 + E + O\bigl(n^{-3/2}\bigr) = 2\sqrt n + c + \frac{1}{2\sqrt n} + O\bigl(n^{-3/2}\bigr),

مع c=E32c = E - \frac32. وعند n=104n = 10^4: 2n=2002\sqrt n = 200، c1.46035c \approx -1.46035، 12n=0.005\frac1{2\sqrt n} = 0.005: فالمتوقع 198.54465198.54465، وفعلًا k104k1/2=198.544645\sum_{k\leq10^4}k^{-1/2} = 198.544645\dots — ثلاثة حدود، وسبعة أرقام.

11. ttlntt \mapsto t\ln t متصلة ومتزايدة تمامًا على [1,)\intco1\infty (بمشتق lnt+11\ln t + 1 \geq 1)، من 00 إلى ++\infty: فيوجد xnx_n وحيد، وxnx_n \to \infty (وإلا لبقيت xnlnxnx_n\ln x_n محدودة). وبأخذ اللوغاريتمات في xnlnxn=nx_n\ln x_n = n: lnxn+lnlnxn=lnn\ln x_n + \ln\ln x_n = \ln n؛ وبما أن lnlnxn=o(lnxn)\ln\ln x_n = o(\ln x_n)، تعطي القسمة على lnxn\ln x_n أن lnnlnxn1\frac{\ln n}{\ln x_n} \to 1: أي lnxnlnn\ln x_n \sim \ln n.

12. من xn=nlnxnx_n = \frac{n}{\ln x_n} ومن lnxnlnn\ln x_n \sim \ln n: xnnlnnx_n \sim \frac{n}{\ln n}. والمرور التالي: lnlnxn=ln(lnn(1+o(1)))=lnlnn+o(1)\ln\ln x_n = \ln\bigl(\ln n\,(1 + o(1))\bigr) = \ln\ln n + o(1)، ومنه lnxn=lnnlnlnn+o(1)\ln x_n = \ln n - \ln\ln n + o(1) و

xn=nlnnlnlnn+o(1)=nlnn11lnlnn+o(1)lnn=nlnn(1+lnlnnlnn+o(lnlnnlnn)).x_n = \frac{n}{\ln n - \ln\ln n + o(1)} = \frac{n}{\ln n}\cdot\frac{1}{1 - \frac{\ln\ln n + o(1)}{\ln n}} = \frac{n}{\ln n}\Bigl(1 + \frac{\ln\ln n}{\ln n} + o\Bigl(\frac{\ln\ln n}{\ln n}\Bigr)\Bigr).

13. عند n=106n = 10^6: nlnn72382\frac{n}{\ln n} \approx 72\,382 (بخطأ 18%18\%)، ويعطي الحدّان 86140\approx 86\,140 (بخطأ 1.9%1.9\%)، في مقابل القيمة الحقيقية x87848x \approx 87\,848. والمكسب في كل مرور ليس إلا العامل lnlnnlnn2.6313.80.19\frac{\ln\ln n}{\ln n} \approx \frac{2.63}{13.8} \approx 0.19: فالسلالم اللوغاريتمية تتقارب ببطء يثير الجنون — وهي واقعة حياة أينما كانت الأعداد الأولية حاضرة.

14. توجد nn عددًا أوليًا بالضبط pn\leq p_n (وهي p1,,pnp_1, \dots, p_n): أي π(pn)=n\pi(p_n) = n. وتعطي مبرهنة الأعداد الأولية (المقبولة؛ مجلد السنة الثالثة) أن n=π(pn)pnlnpnn = \pi(p_n) \sim \frac{p_n}{\ln p_n}، أي pnnlnpnp_n \sim n\ln p_n: وهذه هي المعادلة xlnxnx\ln x \approx n مقروءةً في الاتجاه المعاكس. وبأخذ اللوغاريتمات: lnpn=lnn+lnlnpn+o(1)\ln p_n = \ln n + \ln\ln p_n + o(1)، ويفرض lnlnpn=o(lnpn)\ln\ln p_n = o(\ln p_n) أن lnpnlnn\ln p_n \sim \ln n كما في السؤال 11. وبالتعويض رجوعًا:

pnnlnpn=nlnnlnpnlnnnlnn.p_n \sim n\ln p_n = n\,\ln n\,\frac{\ln p_n}{\ln n} \sim n\ln n .

15. (أ) نثبّت ε>0\varepsilon > 0؛ ومن أجل kk الكبيرة، (1ε)klnkpk(1+ε)klnk(1 - \varepsilon)k\ln k \leq p_k \leq (1 + \varepsilon)k\ln k. وبالمقارنة مع الدالة المتزايدة tlntt\ln t (بحصر من نمط المبرهنة 6.6knklnk=1ntlnt ⁣dt+O(nlnn)=n2lnn2n24+O(nlnn)n2lnn2\sum_{k\leq n}k\ln k = \int_1^n t\ln t\,\dd t + O(n\ln n) = \frac{n^2\ln n}2 - \frac{n^2}4 + O(n\ln n) \sim \frac{n^2\ln n}2. ومنه knpk=n2lnn2(1+O(ε)+o(1))\sum_{k\leq n}p_k = \frac{n^2\ln n}{2}(1 + O(\varepsilon) + o(1)) من أجل كل ε\varepsilon: أي knpkn2lnn2\sum_{k\leq n}p_k \sim \frac{n^2\ln n}2. (ب) بمبرهنة الأعداد الأولية، تكون بين الأعداد الصحيحة حتى 1010010^{100} نسبةٌ 1ln10100=1230.26\sim \frac{1}{\ln 10^{100}} = \frac1{230.26\dots} أولية: فاحتمال أن يكون عدد صحيح عشوائي منتظم ذو 100100 رقمًا أوليًا هو نحو 1230\frac1{230}.

16. على (nπ,nπ+π2)\intoo{n\pi}{n\pi + \frac\pi2}، تكون g(x)=tanx1xg(x) = \tan x - \frac1x متصلة ومتزايدة تمامًا (لأن g=1+tan2x+1x2>0g' = 1 + \tan^2x + \frac1{x^2} > 0)، بالقيمة g1nπ<0g \to -\frac1{n\pi} < 0 عند الطرف الأيسر وg+g \to +\infty عند الأيمن: فجذر واحد بالضبط xnx_n. وبما أن tanzn=tanxn=1xn0\tan z_n = \tan x_n = \frac1{x_n} \to 0 مع zn(0,π2)z_n \in \intoo{0}{\frac\pi2}: فإن zn0+z_n \to 0^+.

17. tanznzn\tan z_n \sim z_n و1xn1nπ\frac1{x_n} \sim \frac1{n\pi}: ومنه zn1nπz_n \sim \frac1{n\pi}.

18. zn=arctan1xnz_n = \arctan\frac1{x_n} وarctanu=u+O(u3)\arctan u = u + O(u^3). ومع zn=O(1n)z_n = O(\frac1n):

1xn=1nπ11+znnπ=1nπznn2π2+O(1n4)=1nπ+O(1n3),\frac1{x_n} = \frac{1}{n\pi}\cdot\frac1{1 + \frac{z_n}{n\pi}} = \frac1{n\pi} - \frac{z_n}{n^2\pi^2} + O\Bigl(\frac1{n^4}\Bigr) = \frac1{n\pi} + O\Bigl(\frac1{n^3}\Bigr),

ومنه zn=1nπ+O(1n3)z_n = \frac1{n\pi} + O\bigl(\frac1{n^3}\bigr): فدرجة السلّم cn2\frac{c}{n^2} تحمل المعامل 00، لأن أول تصحيح للمقدار 1xn\frac1{x_n} هو نفسه من الحجم znn2=O(n3)\frac{z_n}{n^2} = O(n^{-3}).

19. نُدخِل zn=1nπ+O(n3)z_n = \frac1{n\pi} + O(n^{-3}) في الصيغة السابقة:

1xn=1nπ1n3π3+O(1n5),\frac{1}{x_n} = \frac{1}{n\pi} - \frac{1}{n^3\pi^3} + O\Bigl(\frac1{n^5}\Bigr),

ثم zn=arctan1xn=1xn13(1xn)3+O(1n5)=1nπ1n3π313n3π3+O(1n5)z_n = \arctan\frac1{x_n} = \frac1{x_n} - \frac{1}{3}\Bigl(\frac1{x_n}\Bigr)^3 + O\Bigl(\frac1{n^5}\Bigr) = \frac1{n\pi} - \frac{1}{n^3\pi^3} - \frac{1}{3n^3\pi^3} + O\Bigl(\frac1{n^5}\Bigr):

xn=nπ+1nπ43π3n3+O(1n5).x_n = n\pi + \frac{1}{n\pi} - \frac{4}{3\pi^3n^3} + O\Bigl(\frac1{n^5}\Bigr).

20. عند n=3n = 3: حدٌّ واحد يعطي 9.530889.53088، وثلاثة حدود تعطي 9.529299.52929، والجذر الحقيقي 9.529339.52933: فالخطآن 1.51031.5\cdot10^{-3} و 51055\cdot10^{-5}. والمقابلة: من أجل tanx=x\tan x = x يجب أن يجعل الجذر tan\tan ضخمًا، ومن ثم يلتصق بالطرف الأيمن nπ+π2n\pi + \frac\pi2 من النافذة، على مسافة 1nπ\sim\frac1{n\pi} قبل المقارب؛ ومن أجل xtanx=1x\tan x = 1 يجب أن يجعل الجذر tan\tan ضئيلًا، ومن ثم يقع بُعيد الطرف الأيسر nπn\pi، على مسافة 1nπ\sim\frac1{n\pi} بعد الصفر. الطريقة نفسها، وجغرافيا معكوسة.

21. sinu<u\sin u < u على (0,π)\intoo0\pi، وsin\sin يرسل (0,π)\intoo0\pi في (0,1](0,π)\intoc01 \subseteq \intoo0\pi: فبعد خطوة واحدة u1(0,1]u_1 \in \intoc{0}{1}، ثم تتناقص (un)(u_n) وتُحدّ من أسفل بالقيمة 00: فتتقارب إلى نقطة ثابتة للتطبيق sin\sin، أي إلى 00. والنشر: sinu=u(1u26+o(u2))\sin u = u(1 - \frac{u^2}6 + o(u^2))، ومنه

1un+121un2=1un2((1un26+o(un2))21)=1un2(un23+o(un2))13.\frac{1}{u_{n+1}^2} - \frac1{u_n^2} = \frac{1}{u_n^2}\Bigl(\bigl(1 - \tfrac{u_n^2}6 + o(u_n^2)\bigr)^{-2} - 1\Bigr) = \frac{1}{u_n^2}\Bigl(\frac{u_n^2}{3} + o(u_n^2)\Bigr) \longrightarrow \frac13 .

22. بتشيزارو (مجلد السنة الأولى)، يتقارب متوسط التزايدات إلى النهاية نفسها:

1n1un2=1n(1u02+k=0n1(1uk+121uk2))13,\frac{1}{n}\cdot\frac{1}{u_n^2} = \frac1n\Bigl(\frac1{u_0^2} + \sum_{k=0}^{n-1} \Bigl(\frac1{u_{k+1}^2} - \frac1{u_k^2}\Bigr)\Bigr) \longrightarrow \frac13 ,

ومنه un23nu_n^2 \sim \frac3n، وبما أن كل الحدود موجبة، un3/nu_n \sim \sqrt{3/n}.

23. حسب السؤال 7، Hnlnnγ12n18n2\abs{H_n - \ln n - \gamma - \frac1{2n}} \leq \frac1{8n^2}. وعند n=106n = 10^6 يكون هذا الحدّ 181012=1.251013\frac{1}{8\cdot10^{12}} = 1.25\cdot10^{-13}: فثلاثة حدود محسوبة تعطي المجموع التوافقي ذا المليون حدّ إلى ثلاثة عشر رقمًا، بشهادة خطأ صارمة تمامًا — وهو كل مغزى الصيغة المقاربة ذات الباقي الصريح.

24. (أ) صحيح: lnunlnvn=lnunvn0\ln u_n - \ln v_n = \ln\frac{u_n}{v_n} \to 0 في حين أن lnvn+\ln v_n \to +\infty، ومن ثم تؤول نسبة اللوغاريتمين إلى 11. (ب) خاطئ: un=n+1vn=nu_n = n + 1 \sim v_n = n، لكن eun/evn=e1\eu^{u_n}/\eu^{v_n} = \eu \neq 1. فالتكافؤ يحتمل أخطاءً جمعية o(1)o(1) في الأس، لا O(1)O(1). (ج) خاطئ: f(x)=x+sin(x2)g(x)=xf(x) = x + \sin(x^2) \sim g(x) = x عند ++\infty، لكن f(x)=1+2xcos(x2)f'(x) = 1 + 2x\cos(x^2) يتذبذب بلا حدّ في حين أن g=1g' = 1: فلا يلزم أن تكون مشتقات دالتين متكافئتين قابلة للمقارنة إطلاقًا.

25. جرت حلقة الطريقة 6.22 على النحو نفسه ثلاث مرات: وطّن الجذر، واستخرج حدًّا خامًا، وأعِد إدخاله من أجل الرتبة التالية — على ex+x=n\eu^x + x = n (الجزء الأول)، وعلى xlnx=nx\ln x = n (الجزء الثالث)، وعلى xtanx=1x\tan x = 1 (الجزء الرابع). ويرقّي تصحيح شبه المنحرف مقارنةَ المتسلسلة بالتكامل من “الفرق يتقارب” إلى حدٍّ صريح f(1)+f(n)2\frac{f(1) + f(n)}2 بباقٍ مصدَّق O(nf)O(\int_n^\infty \abs{f''}) — أي ثوابت وأشرطة خطأ بدل مجرد التقارب. والجسر إلى الأعداد الأولية قلبٌ محض: فمبرهنة الأعداد الأولية تقول π(x)lnxx\pi(x)\ln x \approx x، ومن ثم فإن pnp_n، المعرَّف بالعلاقة π(pn)=n\pi(p_n) = n، يحلّ معادلة من نمط xlnx=nx\ln x = n — ويرث مقارباتها. وقد شرّعت القاعدة (أ) في السؤال 24 كل انتقال من unvnu_n \sim v_n إلى lnunlnvn\ln u_n \sim \ln v_n (السؤالان 11 و14)؛ وخطأُ (ب) هو السبب في أننا لا نأخذ أسّي التكافؤات أبدًا. والقمتان: صيغة أويلر–ماكلورين من الرتبة الأولى (السؤال 6)، والقانون المقارب pnnlnnp_n \sim n\ln n للعدد الأولي رقم nn (السؤال 14).

المصطلحات المعرَّفة في هذا الفصل

عرض كل المصطلحات (395) في المسرد