在Matlab中解决大型ODE系统

Question

我是Matlab的新手，正在尝试以数值方式解决以下系统：

我可以解决单个S的情况，但是我需要在N = 1000（我从1到N）的情况下进行处理，而且我不确定如何扩展代码。

clearvars;
tSpan = [0 10];

S0 = [1/2*sind(45)*cosd(45) 1/2*sind(45)*sind(45)  1/2*cosd(45)];

[t,S] = ode45(@eom,tSpan,S0);

plot(t,S(:,1),t,S(:,2),t,S(:,3))
legend('Sx','Sy','Sz')

%plot(S(:,1),S(:,2));

在文件eom.m中

function dSdt = eom(t,S)
    G=5;
    k=1;
    deltax = -G*k*S(1);
    deltay =-G*k*S(2);
    B(1) = 2*k*deltax;
    B(2) = 2*k*deltay;
    B(3) = k*k;

    dSdt = [B(2)*S(3)-B(3)*S(2);B(3)*S(1)-B(1)*S(3);B(1)*S(2)-B(2)*S(1)];
end