四阶变步长积分公式的MATLAB代码，急求！

求大神用MATLAB编写一下四阶Cotes的变步长公式推导程序，可以的话，能加一下注释吗？

谢谢，急求！