📄 i0_c.mht
字号:
From: <由 Microsoft Internet Explorer 5 保存>
Subject:
Date: Wed, 27 Sep 2006 22:50:01 +0800
MIME-Version: 1.0
Content-Type: text/html;
charset="gb2312"
Content-Transfer-Encoding: 7bit
Content-Location: http://www.oxbad.com/DSP/maindoc/arithmetic_source/c/I0.C
X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2900.2962
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=gb2312">
<META content="MSHTML 6.00.2900.2963" name=GENERATOR></HEAD>
<BODY><PRE>/* I0.c - Modified Bessel Function I0(x)
*
* I0(x) = \sum_{k=0}^{\infty}[(x/2)^k / k!]^2
*
*/
#include <math.h>
#define eps (1.E-9) /* \(\ep=10\sp{-9}\) */
double I0(x) /* usage: y = I0(x) */
double x;
{
int n = 1;
double S = 1, D = 1, T;
while (D > eps * S) {
T = x / (2 * n++);
D *= T * T;
S += D;
}
return S;
}
</PRE></BODY></HTML>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -