بدأ التحليل المقارب — وهو فنّ تعويض مقدار معقّد بمقدار بسيط مع خطأ محكوم — في مجلد السنة الأولى مع نشور تايلور. ويجعله هذا الفصل تخصّصًا قائمًا بذاته: النشور على سلالم عامة، ومقارنة المتسلسلة بالتكامل بكل قوتها المقاربة، وصيغة ستيرلينغ (مبرهَنةً كاملةً)، والدراسة المنهجية للمتتاليات المعرَّفة ضمنيًا. وهذه التقنيات هي الخبز اليومي للتحليل المقارب، وكل فصل لاحق يقدّر أي شيء — متسلسلات أو تكاملات أو احتمالات — يأكل من هذه المائدة.
6.1 علاقات المقارنة والسلالم
تعريف 6.1
قرب نقطة a (مع a∈R أو ±∞)، ومن أجل الدوال (أو المتتاليات، مع n→∞): تُعرَّف f=o(g) وf=O(g) وf∼g كما في مجلد السنة الأولى. وسلّم المقارنة عند a عائلة من الدوال الموجبة، قابلة للمقارنة مثنى مثنى، ومرتَّبة ترتيبًا كليًا بالعلاقة o(⋅) — والسلّم القياسي عند +∞ هو
xα(lnx)β(α,β∈R),
مرتَّبًا معجميًا في (α,β)، ويُنقَّح عند الحاجة بالأسّيات eγx.
تعريف 6.2(النشر المقارب)
تقبل fالنشر المقارب
f=c1φ1+c2φ2+⋯+ckφk+o(φk)(φi+1=o(φi) في السلّم)
إذا حققت البواقي المتتالية التقديرات المعروضة. وتكون المعاملات عندئذٍ وحيدة: c1=limf/φ1، ثم بالتراجع ci+1=lim(f−∑j≤icjφj)/φi+1.
مثال 6.3
نشور تايلور نشورٌ مقاربة على السلّم (x−a)k عند a. لكن المفهوم أوسع تمامًا: فعند +∞،
x−lnx1=x1⋅1−xlnx1=x1+x2lnx+o(x2lnx),
وهو نشر على سلّم مختلط — ولا تنطبق أي مبرهنة تايلور، بل المتسلسلة الهندسية وحساب المقادير o فقط.
مثال 6.4(السلّم القياسي مرتَّب فعلًا)
يحتاج الادعاء المعجمي في التعريف 6.1 سطرَ برهان لكل حالة. لنقارن xα(lnx)β وxα′(lnx)β′ عند +∞. إذا كان α<α′: فالنسبة xα−α′(lnx)β−β′→0، لأن قوة سالبة للمقدار x تسحق أي قوة للمقدار lnx (ضع x=et: e(α−α′)ttβ−β′→0 بنهاية “الأسّي يغلب كثير الحدود” في مجلد السنة الأولى). وإذا كان α=α′ وβ<β′: فالنسبة (lnx)β−β′→0 مباشرةً. ومن ثم فإن الأزواج (α,β)، مرتَّبةً معجميًا، ترتّب السلّم بالعلاقة o(⋅) — والتعويض x=et هو الحيلة الصالحة لكل المقارنات المختلطة بين القوى واللوغاريتمات.
مثال 6.5(ترتيب حديقة حيوان)
يجب أن تكون السلالم مرتَّبة؛ وهذا هو التمرين القياسي. عند +∞، نقارن n10 وelnn⋅n و2n وnlnn بأخذ اللوغاريتمات:
10lnn≪(lnn)2≪nlnn≪nln2,
حيث تعني an≪bn أن an=o(bn)؛ والمدخل الثاني هو ln(nlnn). وتحافظ الأسّيات على هذه الفجوات التامة (فإذا كان lnun−lnvn→−∞ فإن un/vn→0)، ومنه
n10=o(nlnn),nlnn=o(enlnn),enlnn=o(2n).
والعبرة مزدوجة: قارِن دائمًا عبر اللوغاريتمات (بفروق اللوغاريتمات، لا بنسبها)، ولا تستنتج أبدًا un∼vn من lnun∼lnvn — فالزوج n10 وnlnn نسبة ln فيه تؤول إلى ∞، لكن 2n و4n نسبة ln فيهما تساوي 2 بالضبط وهما غير متكافئين إطلاقًا.
6.2 مقارنة المتسلسلة بالتكامل، مقاربيًا
مبرهنة 6.6
لتكن f متصلة وموجبة ومتناقصة على [1,+∞).
إذا تقاربت ∫1∞f، فإن البواقي تحقق
∫n+1∞f≤k>n∑f(k)≤∫n∞f.
وإذا تباعدت ∫1∞f، فإن المجاميع الجزئية تحقق ∑k=1nf(k)=∫1nf+C+o(1) من أجل ثابت C ما: أي إن الفرق ∑k≤nf(k)−∫1nfيتقارب.
برهان. كان الحصر f(k+1)≤∫kk+1f≤f(k) (بالتناقص) هو حيلة السنة الأولى؛ ويعطي جمعه على k≥n+1 أو على k≥n البندَ (1). ومن أجل (2)، نضع uk=f(k)−∫kk+1f: فبالحصر، 0≤uk≤f(k)−f(k+1)، ومن ثم فإن المجاميع الجزئية للمتسلسلة ∑uk محدودة بالمجموع التلسكوبي f(1)−f(n+1)≤f(1): فتتقارب المتسلسلة. زيادةً على ذلك، تكون المتتالية (∫nn+1f)n غير متزايدة (لأن f متناقصة) وغير سالبة، ومن ثم فهي متقاربة. وبكتابة
k=1∑nf(k)−∫1nf=k=1∑nuk+∫nn+1f,
يتقارب الطرف الأيمن حين n→∞: أي إن الفرق يتقارب إلى ثابت C، وهو البند (2). ∎
مثال 6.7(نشر المتسلسلة التوافقية)
من أجل f(t)=t1: Hn=lnn+γ+o(1)، فنستعيد ثابت أويلر (مجلد السنة الأولى) ببرهان أنظف. وبدفعها رتبة أخرى (التمرين 6.3):
Hn=lnn+γ+2n1+o(n1).
وتُظهر الأعداد المكسب عند n=10: فلدينا H10=2.928968… وln10=2.302585…، ومن ثم فالتقدير الخام للمقدار γ هو H10−ln10=0.626383، بخطأ 0.049؛ وطرح التصحيح 201 يعطي 0.576383، وهو يبعد عن γ=0.577216 بمقدار 8.3⋅10−4 فقط — وهو نفسه الحدّ التالي 12⋅1001 في النشر، كما تبرهن مسألة نهاية الأسبوع (السؤال 8).
مثال 6.8(تقدير خام للمقدار ln(n!) بلا ستيرلينغ)
تحدّد حيلة الحصر وحدها موقع ln(n!) بالفعل. فبما أن ln متزايدة،
∫k−1klntdt≤lnk≤∫kk+1lntdt,
وبالجمع على k=2,…,n (مع ∫1nln=nlnn−n+1):
nlnn−n+1≤ln(n!)≤(n+1)ln(n+1)−n.
والسياجان كلاهما nlnn−n+O(lnn): ومنه ln(n!)=nlnn−n+O(lnn)، وعلى الخصوص ln(n!)∼nlnn. وما تضيفه صيغة ستيرلينغ هو الدرجتان التاليتان — المقدار 21lnn والثابت ln2π — وهما يكلّفان الحصر التلسكوبي الأدقّ في المبرهنة 6.13. ومعرفةُ أي دقة يشتريها كل أداة نصفُ صنعة التحليل المقارب.
مثال 6.9(التلاشي يقتضي النشور)
لنحسب نهاية n2+n−n. فالحدّان كلاهما ∼n، والكتابة “∼n−n” لا معنى لها: إذ لا يمكن طرح المكافئات. فننشر بدل ذلك:
فالنهاية 21، مع سرعة الاقتراب 8n1 مكافأةً. وتستحق الآلية اسمًا: فالفرق بين مقدارين كبيرين متكافئين يعيش كله في حدّيهما التاليين، ومن ثم يجب النشر إلى أول رتبة يختلف عندها الطرفان — مع حمل الباقي للتصديق على أن لا شيء آخر ينجو عند تلك الرتبة.
مثال 6.10(مقارنة متباعدة، منفَّذة)
من أجل f(t)=tlnt1 على [2,+∞) (وهي متصلة وموجبة ومتناقصة): ∫2xf=lnlnx−lnln2→∞، ومن ثم حسب المبرهنة 6.6 (2)،
k=2∑nklnk1=lnlnn+C+o(1)
من أجل ثابت C ما. ودرسان. أولًا، التباعد حقيقي لكنه جليدي: فالمجموع الجزئي يتجاوز 4 أول مرة قرب n≈ee4−C، وهو عدد فلكي. ثانيًا، أعطت الدالة الأصلية الشكلlnlnn، ولم يُخمَّن: فمن أجل الحدود الرتيبة، يكون التكامل أداة الجمع القانونية، والثابت C — مثل ثابت أويلر γ — هو ذاكرة الحدود الأولى.
6.3 صيغة ستيرلينغ
مبرهنة مساعدة 6.11(تكاملات واليس، من جديد)
لتكن Wn=∫0π/2sinntdt. عندئذٍ nWnWn−1=2π من أجل n≥1، و(Wn) متناقصة، وWn∼2nπ.
برهان. تعطي المكاملة بالتجزئة nWn=(n−1)Wn−2 (مع n≥2)، ومن ثم فإن nWnWn−1 ثابت في n، ويساوي 1⋅W1W0=2π. التناقص: sinn+1≤sinn على [0,2π]. والحصر، بالتفصيل: تعطي الرتابة Wn+1≤Wn≤Wn−1، وبالقسمة على Wn−1>0،
n+1n=Wn−1Wn+1≤Wn−1Wn≤1,
والمتطابقة اليسرى من التراجع عند الدليل n+1. ويؤول الحدّان إلى 1: ومنه Wn∼Wn−1، ومن ثم
nWn2∼nWnWn−1=2π⟹Wn∼2nπ.
∎
مثال 6.12(تكاملات واليس الأولى)
من W0=2π وW1=1 والتراجع nWn=(n−1)Wn−2:
W2=4π,W3=32,W4=163π,W5=158,W6=325π.
تحمل الأدلة الزوجية عاملًا π، والفردية ناطقة — وهما الجداءان المتداخلان للصيغ المغلقة. وعدديًا W6≈0.4909 في مقابل المقارب π/12≈0.5116: فعند n=6 يكون المكافئ ضمن 5% بالفعل، والمتطابقة الجدائية دقيقة عند كل n: 6W6W5=6⋅325π⋅158=2π. وجداول صغيرة كهذه أرخص طريقة لاصطياد زلة جبرية قبل أن تعدي حجة مقاربة.
مبرهنة 6.13(ستيرلينغ)
n!∼2πn(en)n.
برهان.الخطوة 1: n!∼Cn(n/e)n من أجل ثابت C>0 ما. نضع
بنشر تايلور للدالة ln(1+n1). ومن ثم تتقارب المتسلسلة ∑(dn−dn+1) بإطلاق (بالمقارنة مع ∑n−2)، فتتقارب (dn)، ولتكن نهايتها d؛ وبأخذ الأسّي، n!∼Cn(n/e)n مع C=ed.
الخطوة 2: C=2π عبر واليس. تتضافر الصيغة المغلقة W2p=4p(p!)2(2p)!⋅2π (الآتية من التراجع، وهي حساب السنة الأولى معادًا في إطار المبرهنة المساعدة 6.11) مع الخطوة 1:
واحتمال أن يعود مسير عشوائي متناظر إلى 0 عند الزمن 2n هو ∼πn1 — وهو إعلان عن الفصل 22.
ملاحظة 6.15(آفاق داخل هذا المجلد)
يتكلم كل فصل كمّي آتٍ لغةَ هذا الفصل. فالفصل الفصل 7 يصنّف المتسلسلات بمقارنة حدودها بالسلّم n−α(lnn)−β — وترسم مسألة نهاية الأسبوع فيه تلك الحدود كاملةً. ويفعل الفصل 9 الشيء نفسه من أجل التكاملات المعتلة، بالسلّم نفسه في المتغيّر المتصل. ويحسب الفصل 11 أنصاف أقطار التقارب من limsup∣an∣1/n، وهو تمرين في مكافئ الجذور من الرتبة n يكون فيه ستيرلينغ المفتاحَ القياسي (nn!∼en، التمرين 6.4). وتصرف فصول الاحتمال ستيرلينغ مباشرةً: فالتقديرات المحلية في الفصل 22 للمعاملات الثنائية هي المثال 6.14 و المثال 6.21 حرفيًا. فالتحليل المقارب ليس فصلًا هنا؛ بل هو لكنة المجلد.
طريقة 6.16(قائمة تدقيق الإقلاع المتدرّج)
قبل الوثوق بنشر متدرّج، دقّق أربع نقاط. (1) الوجود أولًا: يجب تثبيت الجذر أو المتتالية (بالرتابة والقيم الوسطى) قبل أي نشر — فالرموز بلا مرجع تنتشر بجمال ولا تعني شيئًا. (2) رتبة واحدة في كل مرور: لا يُوثق بكل تعويض إلا إلى رتبة التقدير المُدخَل؛ واستخراج حدّين جديدين من مرور واحد هو المصدر الكلاسيكي للمعاملات الخاطئة. (3) البواقي ترافقك: احمل المقدار o(⋅) عبر كل خطوة جبرية ودَع الابتلاع (ابتلاع البواقي الأكبر للحدود الأصغر) يقع في النهاية، صراحةً. (4) التدقيق العددي: قيّم عند قيمة أمينة واحدة للمقدار n؛ فخطأ المعامل ينجو من إعادة الاستنباط الجبري كثيرًا على نحو مفاجئ، ولا ينجو من الحساب أبدًا تقريبًا.
ملاحظة 6.17(مزالق شائعة)
(1) تُجمَع المكافئات جمعًا سيئًا: فمن un∼n+lnn و vn∼−nلا يجوز استنتاج un+vn∼lnn؛ فالتلاشيات تقتضي نشورًا ببواقٍ صريحة، لا مكافئات مجرّدة. (2) لا تأخذ أسّي تكافؤ أبدًا: فلدينا n+1∼n لكن en+1∼en؛ والاتجاه الآمن هو أخذ لوغاريتمات المكافئات المؤولة إلى +∞ (مسألة نهاية الأسبوع في هذا الفصل، السؤال 24). (3) النشر المقارب مرتبط بسلّم: فكتابة f=x1+o(x21) تدّعي أكثر مما تدّعيه f=x1+o(x1)، وخلط الاثنين يُبطل الجبر اللاحق. (4) في الإقلاع المتدرّج، عوّض النشر الجاري كله، بما فيه الباقي — فإسقاط مقدار o(⋅) في منتصف المرور يُنتج معاملات معقولة لكنها خاطئة. (5) تحتاج مقارنة المتسلسلة بالتكامل إلى الرتابة: فهي تفشل تمامًا من أجل الحدود المتذبذبة (قارِن ∑ksink، الفصل 7).
مثال 6.18(ستيرلينغ بالأعداد)
عند n=10: تعطي الصيغة 20π(10/e)10≈3598696 في مقابل 10!=3628800: أي خطأ نسبي 8.3⋅10−3، وهو مذهل من أجل قول “مقارب” عند n=10. وللخطأ بنية — وهي التنقيح الدقيق n!=2πn(n/e)n(1+12n1+O(n−2)) — ويشرح تصحيحه الأول 1201≈8.3⋅10−3 الفجوةَ المرصودة بدقة تقريبًا. وآلة أويلر–ماكلورين في مسألة نهاية الأسبوع هي بالضبط المصدر المنهجي لحدود التصحيح هذه.
ملاحظة 6.19(أين يُستعمَل هذا الفصل)
المقارنة المقاربة هي نحو كل ما هو كمّي فيما بعد: محكّات التقارب وبانوراما برتران في الفصل 7، ومحكّات قابلية المكاملة في الفصل 9، وحسابات نصف قطر التقارب في الفصل 11، ومبرهنات النهاية في الفصل 22 (حيث يُشغّل ستيرلينغ تقديرات دي موافر–لابلاس). ويُصنّع مجلد السنة الثالثة الفكرةَ الواحدة التي نبرهن عليها هنا يدويًا — استخرج الحدّ الرئيسي وحُدّ الباقي — فتصير طريقة لابلاس والتقارب المهيمَن عليه.
مثال 6.20(تكامل مقارَن بنفسه: ∫2xlntdt)
تعمل صندوق أدوات المقارنة على التكاملات أيضًا. لتكن F(x)=∫2xlntdt (والدالة المكامَلة متصلة على [2,∞)). نكامل بالتجزئة:
ومنه F(x)∼lnxx. وسيتعرّف القراء الذين لقوا مبرهنة الأعداد الأولية في مسألة نهاية الأسبوع في هذا الفصل على F: فهو التكامل اللوغاريتمي، وهو المقدِّر الأفضل للمقدار π(x)، ويبيّن الحساب أنه يتفق مع lnxx إلى الرتبة الأولى.
مثال 6.21(ستيرلينغ على معامل ثنائي غير متوازن)
يعطي روتين العواملات الثلاثة نفسه المستعمل من أجل المثال 6.14، ومن أجل (n3n)=n!(2n)!(3n)!:
والمعدل الأسّي 427=2233 هو e3nH(1/3) بترميز الإنتروبيا في نظرية المعلومات: فالمعاملات الثنائية غير المتوازنة تنمو أبطأ تمامًا من المعامل المركزي 4n لكل خطوتين — وهنا (27/4)1/3≈1.89<2 لكل خطوة. وكل مقارب ثنائي في التوافقيات والاحتمال (الفصل 22) هو هذا الحساب الواحد بأوزان مختلفة.
6.4 المتتاليات المعرَّفة ضمنيًا
طريقة 6.22
لإيجاد مقاربات حلول xn لمعادلة F(x,n)=0:
وطّن: برهن على وجود xn ووحدانيته في فترة محدَّدة (بالرتابة ومبرهنة القيم الوسطى)، وجد سلوكه الخام (النهاية، ورتبة النموّ).
أقلِع تدريجيًا: عوّض الصورة الخامة xn=(الحدّ الرئيسي)(1+εn) في المعادلة وحُلّ من أجل الرتبة التالية للمقدار εn؛ وكرّر، فينقّح كل مرور رتبةً واحدة.
مثال 6.23
من أجل n≥1، للمعادلة tanx=x حلٌّ واحد بالضبط xn في (nπ−2π,nπ+2π) (فالدالة tanx−x متزايدة من −∞ إلى +∞ هناك، ومشتقها tan2x≥0). الخام:xn=nπ+2π−yn مع yn∈(0,π)؛ وبما أن xn→∞ وtanxn=xn→+∞، فإن xn تقترب من المقارب من اليسار: yn→0. الإقلاع المتدرّج:tanxn=cotyn=tanyn1∼yn1، وتعطي المعادلة cotyn=xn∼nπ أن yn∼nπ1. ومنه
حُلّ x+lnx=n مقاربيًا. التوطين:x↦x+lnx متزايدة من −∞ إلى +∞ على (0,+∞): فجذر وحيد xn، وxn→∞. الخام: يعطي lnxn=o(xn) أن xn∼n. الإقلاع المتدرّج: من xn=n−lnxn ومن lnxn=lnn+o(1) (لوغاريتمات مكافئات، وكلا الطرفين →∞):
xn=n−lnn+o(1);
وبمرور آخر، مع lnxn=ln(n−lnn+o(1))=lnn−nlnn+o(nlnn):
xn=n−lnn+nlnn+o(nlnn).
(وللتحقق عند n=100: الجذر هو x≈95.4415؛ وتعطي صيغة الحدود الثلاثة 100−4.6052+0.0461=95.4409، وصيغة الحدّين 95.3948 — فيكسب كل مرور الرتبة المتوقعة.) الحلقة نفسها، ومنظر ثالث: فطريقة الطريقة 6.22 لا تعبأ بشكل المعادلة، بل بأن يعزل كل مرور المجهولَ المهيمن فحسب.
6.5 تمارين
تمرين 6.1★
انشر عند +∞، حدّين بعد الحدّ الرئيسي:
x2+x+1,ln(x2+x)−2lnx,x−lnxx+sinx.
حل
حل التمرين 6.1.
x2+x+1=x1+x1+x21=x+21+83⋅x1+o(x1) (بالنشر الثنائي: يعطي 21u−81u2 مع u=x1+x21 أن 2x1+2x21−8x21=2x1+8x23، ثم نضرب في x).
ونرتّب الإسهامات على السلّم عند +∞: xlnx≫x1≥xsinx≫x2(lnx)2. ومن ثم فالحدّان اللذان يليان الحدّ الرئيسي 1 هما xlnx، ثم حدّ التذبذب المحدود xsinx:
x−lnxx+sinx=1+xlnx+xsinx+O(x2(lnx)2).
تمرين 6.2★
أعطِ الطبيعة (تقارب أو تباعد) ومقاربات الحدّ الرئيسي عند التباعد للمقدار ∑k≤nkα من أجل α>−1 وα=−1 وα<−1، عبر المبرهنة 6.6.
حل
حل التمرين 6.2.
f(t)=tα (t≥1).
α>−1: تباعد، وحسب المبرهنة 6.6 (2)، ∑k≤nkα=α+1nα+1+C+o(1) إذا كان α<0 (حيث تتناقص f)؛ ومن أجل α≥0 (حيث تتزايد f) يعطي الحصر نفسه بمتراجحات معكوسة أن ∑k≤nkα∼α+1nα+1.
من أجل n≥2، برهن على أن للمعادلة xn+x=1 حلًا وحيدًا xn∈(0,1)، وأن xn→1، وأثبت
xn=1−nlnn+o(nlnn).
(من xnn=1−xn: خذ اللوغاريتمات وأقلِع تدريجيًا مع xn=1−εn.)
حل
حل التمرين 6.5.
g(x)=xn+x−1 متزايدة تمامًا على [0,1] من −1 إلى 1: فجذر وحيد xn. وبما أن xnn=1−xn∈(0,1): فلو كان xn≤c<1 على متتالية جزئية، لكان xnn≤cn→0، ومنه 1−xn→0: وهو تناقض مع xn≤c. ومنه xn→1.
نكتب xn=1−εn، εn→0+. وتُقرأ المعادلة (1−εn)n=εn، أي
nln(1−εn)=lnεn⟹−nεn(1+o(1))=lnεn.
ومنه nεn=−lnεn(1+o(1))→+∞، و بأخذ اللوغاريتمات من جديد: lnn+lnεn=ln(−lnεn)+o(1). وبما أن ln(−lnεn)=o(ln(1/εn))، يعطي هذا أن lnεn∼−lnn، ومنه εn=n−lnεn(1+o(1))∼nlnn:
(أخرِج العامل n: Sn=n1∑k(1+nklnn)−1؛ وتعرّف على مجموع من نمط ريمان بوسيط بطيء التغيّر t=lnn، واحسب ∫011+tudu=tln(1+t)، واختم بأن Sn∼lnnlnlnn.)
حل
حل التمرين 6.9.
نُخرِج العامل n ونضع t=lnn:
Sn=n1k=1∑n1+tnk1.
ومن أجل t ثابت، يكون المجموع مجموعَ ريمان للدالة u↦1+tu1 على [0,1]؛ والدالة رتيبة في u، ومن ثم يُحصر مجموع ريمان بالتكامل مزاحًا بخطوة واحدة:
∫011+tudu−n1≤Sn≤∫011+tudu+n1
(وهي مقارنة مجاميع ريمان لدالة رتيبة بتكاملها، وهي صالحة من أجل كل n بالعدد t=lnn الخاص به). والآن ∫011+tudu=tln(1+t)، و n1=o(tlnt): ومنه
Sn=lnnln(1+lnn)+O(n1)∼lnnlnlnn.
تمرين 6.10★
برهن على المتطابقة (lnn)lnn=nlnlnn، ثم رتّب ما يلي ترتيبًا تصاعديًا بالعلاقة o(⋅) عند اللانهاية، مع البراهين: n2، (lnn)lnn، 2n، n!، nn.
حل
حل التمرين 6.10.
المتطابقة: (lnn)lnn=elnnlnlnn=(elnn)lnlnn=nlnlnn. والترتيب: نقارن اللوغاريتمات. ln(n2)=2lnn؛ ln((lnn)lnn)=lnnlnlnn؛ ln(2n)=nln2؛ ln(n!)=nlnn−n+O(lnn) (بستيرلينغ، أو بالحصر الأخشن lnn!∼nlnn)؛ ln(nn)=nlnn. وبما أن 2lnn=o(lnnlnlnn)، lnnlnlnn=o(n)، nln2=o(nlnn−n)، وnlnn−n∼nlnn لكن n!/nn→0 (ففرق اللوغاريتمين −n+O(lnn)→−∞):
n2=o((lnn)lnn),(lnn)lnn=o(2n),2n=o(n!),n!=o(nn).
(ففي كل خطوة: يؤول فرق اللوغاريتمين إلى +∞، ومن ثم تؤول النسبة إلى 0.)
تمرين 6.11★★
(ذيل المقدار ∑1/k2، بحدّين) باستعمال التلسكوب الدقيق ∑k>nk(k+1)1=n+11 والتفكيك k21=k(k+1)1+k2(k+1)1، برهن على أن
k>n∑k21=n1−2n21+O(n31).
حل
حل التمرين 6.11.
نفكك k21=k(k+1)1+k2(k+1)1 ونجمع من أجل k>n:
k>n∑k21=n+11+k>n∑k2(k+1)1,
فالمجموع الأول تلسكوبي بالضبط (k(k+1)1=k1−k+11). وأما الثاني: فلدينا k2(k+1)1=k31+O(k41) (لأن k2(k+1)1−k31=k3(k+1)−1)، وبمقارنة التكامل ∑k>nk31=2n21+O(n31)، ∑k>nk41=O(n31). ومنه
ويعطي جمع vk+1−vk=1+O(1) أولًا أن vn=n+O(n)، ومنه vn≥cn ابتداءً من رتبة ما؛ وبإعادة الجمع مع 2vk1=O(k1) نجد vn=n+O(lnn). وبمرور آخر: 2vk1=2k1(1+O(klnk))، ومنه
6.6 مسألة: الإقلاع المتدرّج، من أويلر–ماكلورين إلى الأعداد الأولية
نادرًا ما يسلّم مقدار ضمني أو متراكم مقارباته دفعةً واحدة؛ بل تُستخرج على مرورات، يُغذّي كل مرور فيها العلاقةَ المعرِّفة بالتقدير السابق. وتدرّب مسألة نهاية الأسبوع هذه تلك الحلقة على معادلات جديدة، وتبرهن على صيغة أويلر–ماكلورين من الرتبة الأولى (وهي ترقية شبه المنحرف لمقارنة المتسلسلة بالتكامل، بأشرطة خطأ صارمة)، وتقلب xlnx=n، وتصرف أشهر شيك للطريقة: فمن مبرهنة الأعداد الأولية المقبولة، نستخرج القانون المقارب pn∼nlnn للعدد الأولي رقم n.
مسألة 6.1
مسألة نهاية الأسبوع — تصحيح أويلر–ماكلورين ومقاربات العدد الأولي رقم n
الجزء الأول — حلقة الإقلاع المتدرّج على معادلة جديدة.
برهن على قول الوحدانية في التعريف 6.2: إذا كان f=∑i≤kciφi+o(φk)=∑i≤kci′φi+o(φk) على السلّم نفسه، فإن ci=ci′ من أجل كل i. ثم ادفع مثال الدرس المختلط رتبة أخرى:
x−lnx1=x1+x2lnx+x3(lnx)2+o(x3(lnx)2)(x→+∞),
واشرح لماذا لا يظهر أي حدّ x2c.
بيّن أنه من أجل كل n≥1 يوجد للمعادلة ex+x=n حلٌّ حقيقي واحد بالضبط xn، وأن xn→+∞ مع xn∼lnn.
أقلِع تدريجيًا مرتين:
xn=lnn−nlnn−2n2(lnn)2+o(n2(lnn)2).
تحقق عدديًا عند n=1000: قارِن x1000≈6.90083 بالقيم ذات الحدّ الواحد والحدّين والثلاثة في السؤال 3، إلى خمسة أرقام عشرية.
الجزء الثاني — أويلر–ماكلورين، من الرتبة الأولى.
برهن على متطابقة نواة شبه المنحرف: من أجل g من الصنف C2 على [0,1]،
∫01g(t)dt=2g(0)+g(1)−21∫01t(1−t)g′′(t)dt
(كامِل 21t(1−t)g′′ بالتجزئة مرتين).
لتكن f من الصنف C2 على [1,+∞) مع ∫1∞∣f′′∣<∞. بيّن أن
En=k=1∑nf(k)−∫1nf−2f(1)+f(n)
تتقارب إلى ثابت E، مع حدّ الذيل ∣E−En∣≤81∫n∞∣f′′∣: وهي صيغة أويلر–ماكلورين من الرتبة الأولى.
اختبر عند n=106: الجذر الحقيقي هو x≈87848؛ قارِن بالقيمة ذات الحدّ الواحد (≈72382) وبالقيمة ذات الحدّين (≈86140)، واشرح بطء المكسب (فوسيط النشر هو lnnlnlnn، ولا يساوي إلا ≈0.19 عند n=106).
نقبل الآن مبرهنة الأعداد الأولية: يحقق عدد π(x) الأعداد الأولية التي لا تتجاوز ≤x العلاقةَ π(x)∼lnxx حين x→∞ (وهي مبرهَنة بأمانة في مجلد السنة الثالثة). وبكتابة pn للعدد الأولي رقم n، برّر π(pn)=n، وشغّل قلب السؤالين 11–12 لتبرهن على أن
pn∼nlnn.
الأرباح: (أ) بيّن أن ∑k≤npk∼2n2lnn(قارِن ∑klnk بالمقدار ∫tlntdt)؛ (ب) احسب الفرصة التقريبية لأن يكون عدد صحيح عشوائي منتظم ذو 100 رقمًا أوليًا (ln10100≈230.26: نحو واحد من كل 230).
الجزء الرابع — تصدير الطريقة: xtanx=1.
بيّن أنه من أجل كل n≥1 يوجد للمعادلة tanx=x1 حلٌّ واحد بالضبط xn في (nπ,nπ+2π)، وأن zn=xn−nπ→0+.
حدٌّ واحد: zn∼nπ1.
بيّن أن نشر znلا يحتوي أي حدّ n2c: أي zn=nπ1+O(n31).
ثلاثة حدود: باستعمال arctanu=u−3u3+O(u5) وxn1=nπ1−(nπ)2zn+O(n−3⋅zn2)، برهن على أن
xn=nπ+nπ1−3π3n34+o(n31).
تحقق عند n=3: الجذر الحقيقي x3≈9.5293344؛ قارِن القيمتين ذات الحدّ الواحد وذات الحدود الثلاثة، وقابِل بجملة واحدة بينها وبين المتتالية tanx=x في الدرس (المثال 6.23): أين تقع كل متتالية في نافذتها، ولماذا.
الجزء الخامس — إقلاع متدرّج حركي، وقواعد اللعبة، والتركيب.
لتكن u0∈(0,π) وun+1=sinun. بيّن أن un→0 تناقصيًا، واحسب نهاية un+121−un21(انشر sin−2 عبر sinu=u−6u3+o(u3)).
استنتج، عبر متوسطات تشيزارو (مجلد السنة الأولى)، الكلاسيكيةَ
un∼n3.
(حسابات مصدَّقة) باستعمال الحدّ الصارم في السؤال 7، بيّن أن تقييم lnn+γ+2n1 عند n=106 يعطي H106 بخطأ لا يتجاوز 1.25⋅10−13 — أي مجموع بمليون حدّ محسوبًا إلى ثلاثة عشر رقمًا بثلاثة حدود.
(قواعد اللعبة) برهن أو ادحض، بالبراهين أو بالأمثلة المضادة: (أ) إذا كان un∼vn→+∞ فإن lnun∼lnvn؛ (ب) إذا كان un∼vn فإن eun∼evn؛ (ج) إذا كان f∼g عند +∞ (مع f,g قابلة للاشتقاق) فإن f′∼g′.
(تركيب) بجملة واحدة لكل بند: حلقة الإقلاع المتدرّج في الطريقة 6.22 كما استُعملت في الأجزاء الأول والثالث والرابع؛ وماذا يضيف تصحيح شبه المنحرف إلى المبرهنة 6.6؛ ولماذا يكون قلب xlnx هو بالضبط الجسر من π(x) إلى pn؛ وأي قاعدة من قواعد السؤال 24 حرست أي خطوة. وسمِّ القمتين: صيغة أويلر–ماكلورين (من الرتبة الأولى)، والقانون المقارب للعدد الأولي رقم n.
حل
حل المسألة 6.1.
1. بطرح النشرين: ∑i(ci−ci′)φi=o(φk). فإذا اختلف معامل ما، وليكن i0 أولها: فبالقسمة على φi0 واستعمال φj=o(φi0) من أجل j>i0 نجد ci0−ci0′=o(1): أي صفرًا، وهو تناقض. وأما النشر: فمع u=xlnx→0،
ولا يظهر أي حدّ x2c لأن النشر متسلسلة هندسية في u=xlnx: فكل حدّ يحمل من قوى lnx بقدر ما يحمل من قوى x1 بعد الأولى؛ ودرجة السلّم x21 (أي معامل (lnx)0) غائبة ببساطة، بمعامل 0.
2.f(x)=ex+x متصلة ومتزايدة تمامًا، ونهايتاها −∞ و+∞: فهي تقابل R→R، ومن ثم فإن xn=f−1(n) موجود ووحيد، و xn→+∞ (لأن f−1 متزايدة إلى +∞). ومن exn=n−xn: xn=ln(n−xn)≤lnn، ومنه xn/n→0 وxn=lnn+ln(1−xn/n)=lnn+o(1)∼lnn.
4. عند n=1000: ln1000≈6.90776 (بخطأ 7⋅10−3)؛ وبحدّين: 6.90085 (بخطأ 2⋅10−5)؛ وبثلاثة حدود: 6.90082 (بخطأ دون 10−5)، في مقابل x1000≈6.90083. ويشتري كل مرور العامل المتوقع nlnn تقريبًا.
5. مكاملتان بالتجزئة، ابتداءً من اليمين: مع dtd[21t(1−t)]=21−t وt(1−t) منعدمة عند الطرفين،
مع c=E−23. وعند n=104: 2n=200، c≈−1.46035، 2n1=0.005: فالمتوقع 198.54465، وفعلًا ∑k≤104k−1/2=198.544645… — ثلاثة حدود، وسبعة أرقام.
11.t↦tlnt متصلة ومتزايدة تمامًا على [1,∞) (بمشتق lnt+1≥1)، من 0 إلى +∞: فيوجد xn وحيد، وxn→∞ (وإلا لبقيت xnlnxn محدودة). وبأخذ اللوغاريتمات في xnlnxn=n: lnxn+lnlnxn=lnn؛ وبما أن lnlnxn=o(lnxn)، تعطي القسمة على lnxn أن lnxnlnn→1: أي lnxn∼lnn.
12. من xn=lnxnn ومن lnxn∼lnn: xn∼lnnn. والمرور التالي: lnlnxn=ln(lnn(1+o(1)))=lnlnn+o(1)، ومنه lnxn=lnn−lnlnn+o(1) و
13. عند n=106: lnnn≈72382 (بخطأ 18%)، ويعطي الحدّان ≈86140 (بخطأ 1.9%)، في مقابل القيمة الحقيقية x≈87848. والمكسب في كل مرور ليس إلا العامل lnnlnlnn≈13.82.63≈0.19: فالسلالم اللوغاريتمية تتقارب ببطء يثير الجنون — وهي واقعة حياة أينما كانت الأعداد الأولية حاضرة.
14. توجد n عددًا أوليًا بالضبط ≤pn (وهي p1,…,pn): أي π(pn)=n. وتعطي مبرهنة الأعداد الأولية (المقبولة؛ مجلد السنة الثالثة) أن n=π(pn)∼lnpnpn، أي pn∼nlnpn: وهذه هي المعادلة xlnx≈n مقروءةً في الاتجاه المعاكس. وبأخذ اللوغاريتمات: lnpn=lnn+lnlnpn+o(1)، ويفرض lnlnpn=o(lnpn) أن lnpn∼lnn كما في السؤال 11. وبالتعويض رجوعًا:
pn∼nlnpn=nlnnlnnlnpn∼nlnn.
15. (أ) نثبّت ε>0؛ ومن أجل k الكبيرة، (1−ε)klnk≤pk≤(1+ε)klnk. وبالمقارنة مع الدالة المتزايدة tlnt (بحصر من نمط المبرهنة 6.6)، ∑k≤nklnk=∫1ntlntdt+O(nlnn)=2n2lnn−4n2+O(nlnn)∼2n2lnn. ومنه ∑k≤npk=2n2lnn(1+O(ε)+o(1)) من أجل كل ε: أي ∑k≤npk∼2n2lnn. (ب) بمبرهنة الأعداد الأولية، تكون بين الأعداد الصحيحة حتى 10100 نسبةٌ ∼ln101001=230.26…1 أولية: فاحتمال أن يكون عدد صحيح عشوائي منتظم ذو 100 رقمًا أوليًا هو نحو 2301.
16. على (nπ,nπ+2π)، تكون g(x)=tanx−x1 متصلة ومتزايدة تمامًا (لأن g′=1+tan2x+x21>0)، بالقيمة g→−nπ1<0 عند الطرف الأيسر وg→+∞ عند الأيمن: فجذر واحد بالضبط xn. وبما أن tanzn=tanxn=xn1→0 مع zn∈(0,2π): فإن zn→0+.
17.tanzn∼zn وxn1∼nπ1: ومنه zn∼nπ1.
18.zn=arctanxn1 وarctanu=u+O(u3). ومع zn=O(n1):
ومنه zn=nπ1+O(n31): فدرجة السلّم n2c تحمل المعامل 0، لأن أول تصحيح للمقدار xn1 هو نفسه من الحجم n2zn=O(n−3).
19. نُدخِل zn=nπ1+O(n−3) في الصيغة السابقة:
xn1=nπ1−n3π31+O(n51),
ثم zn=arctanxn1=xn1−31(xn1)3+O(n51)=nπ1−n3π31−3n3π31+O(n51):
xn=nπ+nπ1−3π3n34+O(n51).
20. عند n=3: حدٌّ واحد يعطي 9.53088، وثلاثة حدود تعطي 9.52929، والجذر الحقيقي 9.52933: فالخطآن 1.5⋅10−3 و 5⋅10−5. والمقابلة: من أجل tanx=x يجب أن يجعل الجذر tan ضخمًا، ومن ثم يلتصق بالطرف الأيمنnπ+2π من النافذة، على مسافة ∼nπ1 قبل المقارب؛ ومن أجل xtanx=1 يجب أن يجعل الجذر tan ضئيلًا، ومن ثم يقع بُعيد الطرف الأيسرnπ، على مسافة ∼nπ1 بعد الصفر. الطريقة نفسها، وجغرافيا معكوسة.
21.sinu<u على (0,π)، وsin يرسل (0,π) في (0,1]⊆(0,π): فبعد خطوة واحدة u1∈(0,1]، ثم تتناقص (un) وتُحدّ من أسفل بالقيمة 0: فتتقارب إلى نقطة ثابتة للتطبيق sin، أي إلى 0. والنشر: sinu=u(1−6u2+o(u2))، ومنه
23. حسب السؤال 7، Hn−lnn−γ−2n1≤8n21. وعند n=106 يكون هذا الحدّ 8⋅10121=1.25⋅10−13: فثلاثة حدود محسوبة تعطي المجموع التوافقي ذا المليون حدّ إلى ثلاثة عشر رقمًا، بشهادة خطأ صارمة تمامًا — وهو كل مغزى الصيغة المقاربة ذات الباقي الصريح.
24. (أ) صحيح: lnun−lnvn=lnvnun→0 في حين أن lnvn→+∞، ومن ثم تؤول نسبة اللوغاريتمين إلى 1. (ب) خاطئ: un=n+1∼vn=n، لكن eun/evn=e=1. فالتكافؤ يحتمل أخطاءً جمعية o(1) في الأس، لا O(1). (ج) خاطئ: f(x)=x+sin(x2)∼g(x)=x عند +∞، لكن f′(x)=1+2xcos(x2) يتذبذب بلا حدّ في حين أن g′=1: فلا يلزم أن تكون مشتقات دالتين متكافئتين قابلة للمقارنة إطلاقًا.
25. جرت حلقة الطريقة 6.22 على النحو نفسه ثلاث مرات: وطّن الجذر، واستخرج حدًّا خامًا، وأعِد إدخاله من أجل الرتبة التالية — على ex+x=n (الجزء الأول)، وعلى xlnx=n (الجزء الثالث)، وعلى xtanx=1 (الجزء الرابع). ويرقّي تصحيح شبه المنحرف مقارنةَ المتسلسلة بالتكامل من “الفرق يتقارب” إلى حدٍّ صريح 2f(1)+f(n) بباقٍ مصدَّق O(∫n∞∣f′′∣) — أي ثوابت وأشرطة خطأ بدل مجرد التقارب. والجسر إلى الأعداد الأولية قلبٌ محض: فمبرهنة الأعداد الأولية تقول π(x)lnx≈x، ومن ثم فإن pn، المعرَّف بالعلاقة π(pn)=n، يحلّ معادلة من نمط xlnx=n — ويرث مقارباتها. وقد شرّعت القاعدة (أ) في السؤال 24 كل انتقال من un∼vn إلى lnun∼lnvn (السؤالان 11 و14)؛ وخطأُ (ب) هو السبب في أننا لا نأخذ أسّي التكافؤات أبدًا. والقمتان: صيغة أويلر–ماكلورين من الرتبة الأولى (السؤال 6)، والقانون المقارب pn∼nlnn للعدد الأولي رقم n (السؤال 14).