clc; 
clear all; 
close all; 
format short g
t=[0:0.001:10]';%vector columna de tiempo
a1=8.034; %parámetro del sistema

y=a1*t.*t.*t; %sistema a identificar\n
fi=[t.*t.*t ];
[m, p]=size(fi);%OJO CHECAR ESTE CODIGO
[~,theta] =mincuad(y,fi);
theta%resultado de la identificación