数值分析中的复化梯形公式和复化simpson公式，matlab程序报错

Question

问题相关代码

function I=T_n(f,a,b,n)%复化梯形
format long           %显示15位双精度
h=(b-a)/n;
x=linspace(a,b,n+1);  %生成n+1个点,这些点的间距为(b-a)/n
y=feval(f,x); 
I=h*(0.5*y(1)+sum(y(2:n))+0.5*y(n+1));


function y=fun2(x)
y=(sin(x))/x;

运行结果及报错内容

>> T_n(@fun2,0,1,10)
索引超出矩阵维度。

出错 T_n (line 6)
I=h*(0.5*y(1)+sum(y(2:n))+0.5*y(n+1));