chap5.cpp

romberger 积分

#include<iostream.h>
#include"math.h"
#include<fstream.h>
#include<iomanip.h>
#include"chap5.h"


void main()
{
	//程序说明：
	cout<<"*****************************************************************************"<<endl;
	cout<<"                         程序说明："<<endl<<endl;
	cout<<"本程序用Romberg求积法计算椭圆x^2/4+y^2=1的周长"<<endl
		<<"使结果具有n位有效数字"<<endl
		<<"首先要求用户输入所要求有效数位数"<<endl
		<<"然后将积分结果输出到屏幕上"<<endl;
    cout<<"*****************************************************************************"<<endl;
    cout<<endl;

	int nn;//用户输入所要求有效数位数

	do
	{
		
		cout<<"请输入积分要求的有效数位：(需为自然数）; 退出请按0"<<endl;
		cin>>nn;

		if(nn>0)
		{
			romberg r( pow(10, (-1)*nn )/2 );
	        r.output();
		}
		
		else if(nn<0)
		{
			cout<<"有效数位数字输入错误！"<<endl;
			continue;
		}			
	}while(nn!=0);

}