55（数值积分）龙贝格积分公式

【问题描述】利用龙贝格积分公式计算函数f(x)=(x^2+x+1)cos(x)，在区间[0, pi/2]范围内的定积分近似值。生成逼近表R，起始时表中只有R(0,0)；迭代一次时，顺序生成R(1,0)和R(1,1)；迭代第二次时，顺序生成R(2,0), R(2,1)和R(2,2)；迭代第J次时，顺序生成R(J,0), R(J,1), …, R(J,J)。当R(J,J)与R(J-1,J-1)的差的绝对值小于给定阈值时，迭代结束。输出迭代次数及R(J,J)的值。

【输入形式】在屏幕上输入所求积分的精度值d。

【输出形式】输出迭代次数，以及积分近似值（保留d位小数）。

【样例1输入】

3

【样例1输出】

3

2.038

【样例1说明】

输入：在屏幕上输入所求积分的精度值d=3，从而计算得到阈值为10^(-3)=0.001。

输出：经过3次迭代，得到满足精度要求的积分近似值为2.038。

【评分标准】根据输入得到的输出准确

#include<iostream>
#include<cmath>
#include<iomanip>
#define pi 3.1415926
using namespace std;
int d;//精度
double D;//阈值D=0.001
double a = 0, b = (pi / 2);
int n = 50;//积分表最大行数
int maxx = 0, p = 0;//maxx是最后迭代表的最大行值 p是迭代次数
double R[50]