%梯形法的计算函数
function f=fun(x)
i = (60-x).^2 + (60-x).*sin(x.^0.5);
f=  i.*(10.*i + 2.*i.^(2/3));