X=[1 2 2 2 1 0 3 3 ]'
N=8;
n=4;
t=0
r=0
l
al=0
for k=n+1:N              
       e=t-n+k
    r0=r+X(k)*X(k)+X(e)*X(e)
end

  R=r0
  L=sqrt(r0)
w=0

for t=1:n
r=0
    for k=n+1:N
       w=t-n+k
    r=r+X(k)*X(k)+X(w)*X(w)
    end
 a=R
  c=r0
 r1=0
 R=sym([a c;c' r])
for k=n+1:N
  r1=r1+X(k-i)*X(k)+X(k-n+i)*X(k-n)
end
if t>1
    c=[c r1]
else c=r1
end
al=L
 l=c'*inv(L')
 Q=sqrt(r-l*l')
 o=zear(k,1)
L=sym([a1 o;l Q])
 
end