código en matlab para la simulación de oscilador de chua

% la primera parte son las ecuaciones de modelo
function dx=chua(t,x)
global alfa beta gamma a b
fx1=b.*x(1)+0.5*(a-b).*(abs(x(1)+1)-abs(x(1)-1));
dx(1)=alfa.*(x(2)-x(1)-fx1);
dx(2)=x(1)-x(2)+x(3);
dx(3)=-beta.*x(2)-gamma.*x(3);
dx=dx';
end
% guarde en una carpeta la primera parte del codigo te va a salir como chua.m
%abra otro nuevo script y copiar y pegar esto y estas listo

close all
clear
clc
global alfa beta gamma a b
%para cambiar el punto de operacion cambia las siguientes variables
alfa=10.61;
beta=17.719;
gamma=0.13699;
a=-1.40726;
b=-0.816308;
X0=[10 0 0];
T=[0 500];
[t x]=ode45(@chua,T,X0);
figure
plot(x(:,1),x(:,2))
figure
plot(x(:,2),x(:,3))
figure
plot(x(:,2),x(:,3))
figure
plot3(x(:,1),x(:,2),x(:,3))