I am trying to apply the steps attached to the pdf file into my code

Question

daniel choudhry 2020-10-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/605272-i-am-trying-to-apply-the-steps-attached-to-the-pdf-file-into-my-code

回答： Alan Stevens 2020-10-5

采纳的回答： Alan Stevens

matlabccc.pdf

在 MATLAB Online 中打开

1.I am having trouble trying to put the steps attached to the pdf file into my matlab code.

2. After translating the steps from the pdf file I would have to put this into command window in order to generate my answer:

>> rand('seed',1)

>>A=rand(10,10);

>>d=my_det(A)

>>dd=det(A)

>>abs(d-dd)/abs(dd)

3. lu_fac_pp(A) is another function i am calling witin this code.

function d = my_det(A)
n = size(A,1);
p=(1:n)';
tr = zeros(1,n);
s=0;
[L, U, p]= lu_fac_pp(A)
for k = 1:n-1
    
    if A(n,n) == 0
        det(A)=0;
    end
    
    
    [r m]=max(abs(A(k:n,k)));
    m=m+k-1;
    
    if (A(m,k)~=0)
        
        if (m~=k)
            A([k m],:)=A([m k],:);
            p([k m])=p([m k]);
        end
        
        i=k+1:n;
        A(i,k)=A(i,k)/A(k,k);
        
        j=k+1:n;
        A(i,j)=A(i,j)-A(i,k)*A(k,j);
    end
end
if A(n,n) ==0
    disp('A is not invertible');
    return
end
L=tril(A,-1)+eye(n,n);
U=triu(A);
    
    
end