Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

barg

Autor NIE ponosi JAKIEJKOLWIEK odpowiedzialności za skutki użytkowania tego programu oraz za wykorzystanie zawartych tu informacji.

A=[4 8 -4
8 17 -1
-4 -1 57]

B=[1,1,1]

X=A^-1*B'

L=zeros(size(A));

n=3;

for k=1:n
s=0;
for p=1:k-1
s=s+L(k,p)^2;
end

L(k,k)=sqrt(A(k,k)-s);

for i=k+1:n
s=0;
for p=1:k-1
s=s+L(i,p)*L(k,p);
end
L(i,k)=(A(i,k)-s)/L(k,k);
end

end

y=zeros(1,n);
x=y;

for j=1:n
s=0;
for i=1:j-1
s=s+L(j,i)*y(i);
end

y(j)=(B(j)-s)/L(j,j);
end

for i=n:-1:1
s=0;
for j=i+1:n
s=s+L(j,i)*x(j);
end
x(i)=(y(i)-s)/L(i,i);
end

L
L*L'
x
X