mat-1.m

global L m P k c

k=100;
c=0.4;
m=0.3;
L=1;
P=99;

t=0:0.01:20;
y=[0.01 0 0 0];
[t,q]=ode45(@derivadanolin,t,y);
figure(1)
subplot(2,1,1)
plot(t,q(:,1),'LineWidth',1)
grid on
set(xlabel('$Tiempo(seg)$'),'interpreter','latex','fontsize',14)
set(ylabel('$\varphi _1(rad)$'),'interpreter','latex','fontsize',14)
title('Respuesta sistema  linealizado','fontsize',16)
subplot(2,1,2)
plot(t,q(:,2),'LineWidth',1)
grid on
set(xlabel('$Tiempo(seg)$'),'interpreter','latex','fontsize',14)
set(ylabel('$\varphi _2(rad)$'),'interpreter','latex','fontsize',14)
title('Respuesta sistema linealizado','fontsize',16)
figure(2)
subplot(1,2,1)
plot(q(:,1),q(:,3),'LineWidth',1)
grid on
set(xlabel('$x_1$'),'interpreter','latex','fontsize',14)
set(ylabel('$x_3$'),'interpreter','latex','fontsize',14)
title('Diagrama de fase','fontsize',16)
subplot(1,2,2)
plot(q(:,2),q(:,4),'LineWidth',1)
grid on
set(xlabel('$x_2$'),'interpreter','latex','fontsize',14)
set(ylabel('$x_4$'),'interpreter','latex','fontsize',14)
title('Diagrama de fase','fontsize',16)