A =[2-14-31;-11213;4233-1;-31324;13-144];
b =[11;14;4;16;18];
c =[11111];
L =zeros(5)+diag(c);
U =zeros(5);
y =zeros(5,1);
x =zeros(5,1);%Doolittle分解
for i =1:5U(1, i)=A(1, i);L(i,1)=A(i,1)/U(1,1);
end
d =0;
e =0;for i =2:5for j = i:5for k =1:(i-1)
d = d +L(i, k)*U(k, j);
e = e +L(j, k)*U(k, i);
end
U(i, j)=A(i, j)- d;L(j, i)=(A(j, i)- e)/U(i, i);
d =0;
e =0;
end
end
%Gauss消去法求解Ly=b
y(1)=b(1);for i =2:5for j =1:(i-1)
d = d +y(j)*L(i, j);
end
y(i)=(b(i)-d)/L(i,i);
d =0;
end
%Gauss消去法求解Ux=y,x即Ax=b的解
x(5)=y(5)/U(5,5);for i =4:-1:1for j =5:-1:(i+1)
d = d +x(j)*U(i, j);
end
x(i)=(y(i)-d)/U(i,i);
d =0;
end