选列主元的LU分解

function [x] = gaussian_elimination_column_pivot(A, b)
    [n, ~] = size(A);
    augmented_matrix = [A, b];  % Augmented matrix [A|b]

    for k = 1:n-1
        % Partial Pivoting: find the row with the largest absolute value in the current column
        [~, max_row] = max(abs(augmented_matrix(k:n, k)));
        max_row = max_row + k - 1;

        % Swap the current row with the row containing the maximum value
        augmented_matrix([k, max_row], :) = augmented_matrix([max_row, k], :);

        % Gaussian elimination with pivoting
        augmented_matrix(k+1:n,k) = augmented_matrix(k+1:n,k)/augmented_matrix(k,k);
        augmented_matrix(k+1:n,k+1:n+1) = augmented_matrix(k+1:n,k+1:n+1) - ...
        augmented_matrix(k+1:n,k) * augmented_matrix(k,k+1:n+1);
    end

    % Backward substitution
    x = augmented_matrix(: ,n+1);
    for j=n:-1:2
        x(j)=x(j)/augmented_matrix(j,j);
        x(1:j-1)=x(1:j-1)-x(j)*augmented_matrix(1:j-1,j);
    end
    x(1)=x(1)/augmented_matrix(1, 1);

end