我不能给出具体的解决方法，你程序提供的东西太少。说说我的想法，S=S+....检验一下是不是满足矩阵的加法和乘法。
改成S(i)=S(i)+(y(i)/(h(i)^3))*p1+(y(i+1)/(h(i)^3))*p2+(m(i)/(h(i)^2))*p3+(m(i+1)/(h(i)^2))*p4;试试
估计是维数不匹配
