数值不稳定算法.txt

来自「舍入误差与数值稳定性 数值稳定算法 数值不稳定算法」· 文本 代码 · 共 17 行

#include<stdio.h>
#include<math.h>
void main()
{
	double y_0=log(6.0/5.0),y_1;
	int n=1;
	printf("y[0]=%-20f",y_0);
	while(1)
	{
		y_1=1.0/n-5*y_0;
		printf("y[%d]=%-20f",n,y_1);
		if(n>=40) break;
		y_0=y_1;
		n++;
		if(n%3==0) printf("\n");
	}
}