Matlab以矩阵形式显示Cotes系数（可以改n的值）

clc,clear 
syms t;
for n=1:8

for k=1:n+1    %计算Cotes系数
    l=1;
    for j=1:n+1
        if(j~=k)
            l=l*(t-j+1)/(k-j);
        end
   
    end
    s(k)=int(l,t,0,n)/n;
    
    B(n,k)=s(k);%以矩阵形式显示C-cotes系数
end

end
B