我將Lotka Voterra模型（捕食捕食者）寫成Scilab函數，並用ODE解決它。我的問題是我想看到流量的演變。我通過在決議中包含流程（在下面的腳本中僅包含一個）而找到了解決方案，但它確實很「重」。任何人有更好的解決方案？ //parameters
x=[0.04,0.005,0.3,0.2]
n = x(1);//birth rate of prey
c = x(2);//captu
這是我的LU分解CROUT方法代碼： function [L, U] = croutreduction(A)
[row,column]=size(A);
L=eye(row,column);
//A = 3x3
if row==3 then
U(1,1)=A(1,1); U(1,2)=A(1,2); U(1,3)=A(1,3);
L