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