数值稳定算法.txt

舍入误差与数值稳定性 数值稳定算法 数值不稳定算法

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