exercise_1.cpp

来自「求解误差」· C++ 代码 · 共 41 行

#include<iostream>
#include <iomanip>
#include<cmath>
#define N 100000000
using namespace std;
double item(int n);
int main()
{	
	cout.precision(12);	
	int i=1,a[2]={0};
	
	double X=1.0;
	for(i=2;i<N;i++)
	{	
		X=X+item(i);
		if(fabs(X-693147190546e-12)<5e-6)
		{
			a[0]=i;
			break;
		}
	}
	X=1.0;
	for(i=2;i<N;i++)
	{	
		X=X+item(i);
		if(fabs(X-693147190546e-12)<5e-7)
		{
			a[1]=i;
			break;
		}
	}
	cout<<"误差小于5e-6,n取值为："<<a[0]<<endl<<"误差小于5e-6,n取值为："<<a[1]<<endl;
	return 0;
}
double item(int k)
{
	double ans=1.0;
	if(k%2==0)ans=-ans/k;
	else      ans=ans/k;
	return ans;
}