R语言 利用pi/4 = 1 – 1/3 + 1/5 – 1/7 + … 公式求pi的近似值，直到误差小于10-6，要用R语言进行编程

Question

如题所述

Answer 1

#include <stdio.h>

int main(void)
{
double pi = 0.0f, x = 0.0f, y = 0.0f, d = 0.0f; int i = 0, j = 0, k = 0;
do {
j = i + 1; /*求第一项分母*/
k = j + 2; /*求第二项分母*/
x = j;/*转化为浮点数*/
y = k;/*转化为浮点数*/
x = 1.0f / x;/*求第一项*/
y = 1.0f / y;/*求第二项*/
d = x - y;/*求两项差*/
pi += d;/*累加到π*/
i += 4;/*进行下一项*/
}while(x > 1e-6); /*直到某一项差绝对值小于目标误差*/
printf("π = %9.8f\n", pi * 4.0f);
return 0;
}

运行结果为:
π = 3.14159065