romberg算法计算积分 matlab程序

fprintf（‘下面将用Romberg算法求某个函数从a到b积分值的近似解。\n‘）
%输入参数
syms x
f=input（‘输入被积函数的表达式:‘）
a=input（‘输入积分下限:‘）
b=input（‘输入积分上限:‘）
epsilon=input（‘输入精度:‘）
%变量初始化。
h=b-a
k=1
T（11）=（h/2）*（subs（fa）+subs（fb））
%迭代过程.
while（1）
    T（k+11）=0.5*（T（k1）+h*sum（subs