اویلر و مک لورن یک روش کارا برای محاسبه عددی مجموع های متناهی و نامتناهی به فرم$ S^{0} = \sum_1^N f(x) $ ارائه دادن. برای توضیح روش مجموعهای کلی تر زیر را در نظر میگیریم
$ S^{k} = \sum_1^N f^{(k)} (x) , k=0,1,...$ که $f:(0, \infty ) \rightarrow R$ یک تابع از کلاس $ C^{ \infty } $ است. با استفاده از فرمول تیلور برای یک تابع اولیه$F$ از$f$ تخمین ابتکاری زیر را به دست می آوریم
$ \int_0^N f(x)dx=F(N)-F(0)= \sum_1^N (F(n)-F(n-1)) \approx \sum_1^N \frac{ F' (n)}{1!} - \frac{F''(n)}{2!} + \frac{F'''(n)}{3!} -... $* یعنی
$ \int_0^N f(x)dx \approx \frac{ S^{0} }{1!} - \frac{ S^{1} }{2!} + \frac{ S^{2} }{3!} -... $ به طور مشابه برای مشتقات$f$ میتوان فرمولهای زیر را به دست آورد
$f(N)-f(0) \approx \frac{ S^{1} }{1!} - \frac{ S^{2} }{2!} + \frac{ S^{3} }{3!} -... $ و
$ f'(N)-f'(0) \approx \frac{ S^{2} }{1!} - \frac{ S^{3} }{2!} + \frac{ S^{4} }{3!} -... $ و...
با استفاده از این روابط میتوانیم$ S^{1} $و$ S^{2} $و ... را از * یکی یکی حذف میکنیم و رابطهی زیر را به دست میآوریم
$ \sum_1^N f(n) \simeq \int_0^N f(x)dx + \sum_1^ \infty a_{k} ( f^{(k-1)}(N) - f^{(k-1)}(0) ) $ که ضرایب$ a_{k} $ به سرعت به سمت صفر همگرا میشوند.
$ a_{1} = \frac{1}{2} , a_{2} = \frac{1}{12} , a_{3} =0, a_{4} = \frac{-1}{720} , a_{5} =0,...$
حالا از این فرمول با $f(x)= (x+1)^{-2} $ و $N=19$ میتونید برای مجموع خواسته شده استفاده کنید. تجزیه و تحلیل خطای روش دشواره و نیاز به مقدمات خیلی زیاده داره و اینجا نمیشه ارائه کرد درواقع به چندجملهای های برنولی مربوط میشه.