clc; 
clear all; 
close all;

global u %variable global almacena la entrada del sistema
format short g
ti=0; 
h=0.001;  
tf = 5;
    ts=ti:h:tf; %vector tiempo de simulación
    opciones=odeset('RelTol',1e-3, 'InitialStep',2.5e-3,'MaxStep',2.5e-3);
    %solución numérica del sistema dinámico lineal
    [t,x]=ode45('cap7_se1',ts,0,opciones);
    xp=cap7_se1(t,x); %se obtiene la velocidad
    fi=[-x, u];%regresor
    %algoritmo de mínimos cuadrados
    [r,theta] =mincuad(xp,fi);
    theta %resultado de la identificación