매트랩 프로그래밍-수치해석

소개글

매트랩 프로그래밍-수치해석

목차

Problem 1. Answer)
Problem 2. Answer)
Problem 3. Answer)
Compare your results with the previous results.

본문내용

Problem 1. Answer)

>> A=[1 1/2 1/3 1/4 1/5; 1/2 1/3 1/4 1/5 1/6; 1/3 1/4 1/5 1/6 1/7; 1/4 1/5 1/6 1/7 1/8; 1/5 1/6 1/7 1/8 1/9]
B=[1; 2; 3; 4; 5]
U=A;
for k=1:5
for i=k+1:5
pi=U(i,k)/U(k,k);
for j=1:5
U(i,j)=U(i,j)-(U(k,j)*pi);
end
end
end
M=U*inv(A);
L=inv(M);
L
U
Z=[0;0;0;0;0];
Z(1)=B(1);
Z(2)=B(2)-Z(1)*L(2,1);
Z(3)=B(3)-Z(2)*L(3,2)-Z(1)*L(3,1);
Z(4)=B(4)-Z(3)*L(4,3)-Z(2)*L(4,2)-Z(1)*L(4,1);
Z(5)=B(5)-Z(4)*L(5,4)-Z(3)*L(5,3)-Z(2)*L(5,2)-Z(1)*L(5,1);
Z
X=[0;0;0;0;0];
X(5)=Z(5)/U(5,5);
X(4)=(Z(4)-U(4,5)*X(5))/U(4,4);
X(3)=(Z(3)-U(3,5)*X(5)-U(3,4)*X(4))/U(3,3);
X(2)=(Z(2)-U(2,5)*X(5)-U(2,4)*X(4)-U(2,3)*X(3))/U(2,2);
X(1)=(Z(1)-U(1,5)*X(5)-U(1,4)*X(4)-U(1,3)*X(3)-U(1,2)*X(2))/U(1,1);
X
A*X

참고 자료

없음