function  theta =mincuadm(y,fi,Nob,p,n)
theta=[1:p]'; %vector columna de parámetros
theta(1)=0;  %condición inicial del vector de parámetros
 psi=zeros(p,n); %vector columna de observaciones
P=eye(p,p)*10e20; %matriz de covarianza P
I=eye(n,n); %matriz identidad
ys=zeros(n,1);

  for k=1:Nob
    for j=1:n
        for i=1:p %se forma el regresor
            psi(i,j)=fi(k+Nob*(j-1),i);
        end
         for i=1:n 
              ys(i,1)=y(k+Nob*(i-1));    
         end
    end
    e=ys-psi'*theta; %error de regresión
    theta= theta+P*psi*(I+psi'*P*psi)^(-1)*e;%vector estimado
    P=P-(P*psi*(I+psi'*P*psi)^(-1)*(psi')*P);%Matriz de covarianza
  end
end