close all
clear all
clc

Time=0:0.2:100;  % Asse di tempi

global beta k1 k2 M1 M2

beta=.4;
k1=1.2;
k2=1.4;
M1=1;
M2=2;

u=[0 0 0 0]; % Per la risposta libera l'ingresso è nullo

x0=[0 1 0 1]'; % Vettore delle condizioni iniziali

OPTIONS=odeset('RelTol',1e-3,'AbsTol',1e-6);
[t,xmolla]=ode23('molla',Time,x0,OPTIONS,u); %Integrazione numerica

%Grafici delle traiettorie degli stati e delle uscite

subplot(2,2,1)
plot(t,xmolla(:,1))
grid
xlabel('t')
ylabel('x1')
subplot(2,2,2)
plot(t,xmolla(:,2))
grid
xlabel('t')
ylabel('x2')
subplot(2,2,3)
plot(t,xmolla(:,3))
grid
xlabel('t')
ylabel('x3')
subplot(2,2,4)
plot(t,xmolla(:,4))
grid
xlabel('t')
ylabel('x4')