📄 neweq1.m
字号:
function OUT=neweq1(t,X)
a=10; b=8/3; c=28; r=-6;
x(1)=X(1);x(2)=X(2);x(3)=X(3);x(4)=X(4);x(5)=X(5);
Q=[ X(6), X(11), X(16), X(21), X(26);
X(7), X(12), X(17), X(22), X(27);
X(8), X(13), X(18), X(23), X(28);
X(9), X(14), X(19), X(24), X(29);
X(10), X(15), X(20), X(25), X(30);];
dx(1)=a*(x(2)-x(1))+x(4)+x(5);
dx(2)=c*x(1)-x(2)-x(1)*x(3)+x(5);
dx(3)=x(1)*x(2)-b*x(3)+x(5);
dx(4)=-x(2)*x(3)+r*x(4);
dx(5)=-3*x(1);
DX1=[dx(1);dx(2);dx(3);dx(4);dx(5)];
J=[ -a, a, 0, 1, 1;
c-x(3), -1, -x(1), 0, 1;
x(2), x(1), -b, 0, 1;
0, -x(3), -x(2), r, 0;
-3, 0, 0, 0, 0;];
F=J*Q;
OUT=[DX1; F(:)];
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -