lorenz.m

来自「本程序是求解神经电另一chay方程的较为简单的例子。还有一个非线性物理中的Lor」· M 代码 · 共 9 行

 function dy = lorenz(t,Y)
dy = zeros(3,1);    % a column vector
sigma  = 10.; b = 8/3.;gama  = 24.74;
 
	 
 dy(1) = sigma*(Y(2) - Y(1));
 dy(2) = gama*Y(1) - Y(2) - Y(1)*Y(3);
 dy(3) = -b*Y(3) + Y(1) * Y(2);