Romberg求定积分 输入：积分区间[a,b]，被积函数f(x,y,z) 输出：积分结果-Romberg integral to seek input : integral interval [a, b], the plot function f (x, y, z) Output : integral results