중앙대/수치해석 과제2 matlab코드포함

목차

1. 풀이
2. 프로그램
3. 프로그램
4. 프로그램

본문내용

2. 프로그램
function [L,U] = lum(A)
[m,n] = size(A);
if m~=n, error('matrix A must be square'); end

save lum_output.txt
fin01 = fopen('lum_output.txt','wt');
iter = 0;
fprintf(fin01, 'U and L (forward elimination)\n')

L=eye(n);

% forward elimination
for k = 1:n-1
for i = k+1:n
factor = A(i,k)/A(k,k);
A(i,k:n)= A(i,k:n)-factor*A(k,k:n);
L(i,k) = factor;

iter = iter+1;
fprintf(fin01,'iter%d\n', iter);
for y = 1:n
for x =1:n
fprintf(fin01, '%f\t', A(y,x));
end

참고 자료

없음