我编的LUD算法 结果出来没办法控制用 format short没有效果
  以下是代码
X=[0.2 0.34 0.23 1.2 -0.45 -1.23 0 0.3 ]'
N=8;
n=4;
e=0.01
t=0
r=0
l=0
al=[]
r01=0
p=[]
r1=0
b=0

for t=n+1:N              
   
    r0=r+X(t)*X(t)+X(t)*X(t)
end

for t=n+1:N
       w=t-1
    r01=r01+X(t)*X(t-1)+X(w)*X(w+1)
 end
  c=r01
  R=r0
  L=sqrt(r0)
w=0

   
for k=1:n
r=0
    for t=n+1:N
       w=t-k
    r=r+X(k)*X(k)+X(w)*X(w)
    end
 a=R
 
 
 
b=c'
 R=sym([a b;c r])
for t=n+1:N
       w=t-k
       p=k-1
    r1=r1+X(k)*X(k-p)+X(w)*X(w+p)
 end

al=L
 l=c*inv(L')

 Q=sqrt(r-l*l')
 o=zeros(k,1)
L=sym([al o;l Q])
 c=sym([c r1])
end

Fai=0

Fai=e\L\L'