如何为符号矩阵获得正确的双重形式？

Question

如何从MatLab中的这个符号矩阵中获得正确的双重形式？我已尝试过所有内容，但我不喜欢使用feval或内联函数，因为它们不推荐使用

这是获取矩阵的代码

\'{a}

并有这样的输入

function T = Romberg (a, b, m, f)

T = zeros(m, m);
T = sym(T);

syms f(x) c h;

f(x) = f;

c = (subs(f,a)+subs(f,b)) / 2;

h = b - a;

T(1,1) = h * c;

som = 0 ;

n = 2;

for i = 2 : m

    h = h / 2;

    for r = 1 : n/2

        som = som + subs(f,(a + 2*(r-1)*h));

        T(i,1) = h * (c + som);

        n = 2*n;

    end

end

r = 1;

for j = 2 : m

    r = 4*r;

    for i = j : m

        T(i,j) = (r * T(i, j-1) - T(i-1,j-1)/(r-1)); 
    end

end

end

我得到一个带有所有

的符号矩阵输出

Romberg(0, 1, 4, '2*x')

我想要双输出。

你能帮帮我吗？ 非常感谢你提前！