файлы матлаб / check

clear,clc;
a=0.95;
b=0.7;
c=0.6;
d=3.5;
e=0.25;
f=0.1;
syms x1 x2 x3;
dx=@(t,x)[(x(3)-b)*x(1)-d*x(2)...
;d*x(1)+(x(3)-b)*x(2)...
;c+a*x(3)-(x(3)^3)/3-(x(1)^2+x(2)^2)*(1+e*x(3))+f*x(3)*x(1)^3];
x0=0;
ddx=0.001;
% Второй метод
[t,y]=ode45(dx, [0 500], [-0.1387; -0.2527; 1.3259]);
[t1,y1]=ode45(dx, [0 500], [0.1917; 0.6143+ddx; 1.7180]);

t2=[0:0.02:500];
z=interp1(t,y,t2);
z2=interp1(t1,y1,t2);
tou=2485;

Mas=[];
% for i=1:1:tou
zch=z2(1867,:);
A11=zch(3)-b;
A12=-d;
A13=zch(1);
A21=d;
A22=zch(3)-b;
A23=zch(2);
% A31=-2*zch(1)-2*zch(1)*zch(3)*e+3*zch(1)^2*zch(3)*f;
A31=-(2*zch(1)+0)*(1+e*zch(3))+3*zch(1)^2*zch(3)*f;
%A32=-2*zch(2)-2*zch(2)*zch(3)*e;
A32=-(0+ 2*zch(2))*(1+e*zch(3));
% A33=a-zch(3)^2- (zch(1)^2)*e-e*zch(2)^2+f*zch(1)^3;
A33=a-zch(3)^2-(zch(1)^2+ zch(2)^2)*(1+e)+f*zch(1)^3;
A=[A11, A12, A13;A21, A22, A23;A31, A32, A33];
Mas=[eig(A)];
% end
% figure(3)
% plot(t2(1:tou),Mas)