各路大神能不能帮我看看这个微分方程组的龙格库塔法数值解要怎么在MATLAB里写出来,十分感谢

方程组如下,MATLAB我还没学到这一块但是要急用了
①dy/dt=0.518x
②dx/dt=1.107x(1-x-y)-0.518x
如果可以的话在MATLAB里打出来截图,谢谢大神!!
(另外我的MATLAB是2019b的)

f = @(t,y)[1.107*y(1)*(1-y(1)-y(2))-0.518*y(1);0.518*y(1)];

x0 = [1;-1];%由于没有给初值,我随便设的,不同初值不同结果

tspan = [0,20];

[t,y] = ode45(f,tspan,x0);

plot(t,y(:,1),'r--',t,y(:,2),'b-')

legend('x','y')


温馨提示:答案为网友推荐,仅供参考