matlab非单一rhs维度问题，求前辈们知点一下

[m0_38143617]

主程序如此，报错的就是（x：，i+1）那一句，不知道怎么改，
function [x,t]=LONGE1(x0,t0,a,b,h)
n=(b-a)/h;
x(1)=x0;
t(1)=t0;
for i=1:n
K1=fun(x(i),t(i));
K2=fun(x(i)+h/2*K1,t(i)+h/2);
K3=fun(x(i)+h/2*K2,t(i)+h/2);
K4=fun(x(i)+h*K3,t(i)+h);
x(:,i+1)=x(:,i)+h*(K1+2*K2+2*K3+K4)/6;
t(i+1)=t(i)+h;
end

end
函数子程序如下
function z=fun(x,t)
ss=[-10 -100 0 0;1 0 0 0;0 0 -5 -25;0 0 1 0];
F=10;
z=ss*x+F*sin(10*t);
就是函数里多了个矩阵就无法执行了