📄 gepivot.mht

📁 it is a very essential matlab code.

💻 MHT

字号:

From: <Saved by Windows Internet Explorer 7>
Subject: 
Date: Tue, 12 May 2009 09:49:51 -0700
MIME-Version: 1.0
Content-Type: text/html;
	charset="Windows-1252"
Content-Transfer-Encoding: 7bit
Content-Location: http://www.mece.ualberta.ca/Courses/mec390/390code/gepivot.m
X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2900.3350

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=windows-1252">
<META content="MSHTML 6.00.6000.16825" name=GENERATOR></HEAD>
<BODY><PRE>% partial pivoting routine called by gauselim
%
% USAGE:   [A, err] = gepivot(A, k)
%
% where        A = augmented matrix
%              k = current pivot row
%            err = error code
%                = 0 = no error
%                = 1 = matrix is singular

function [A, err] = gepivot(A, k)

err = 0;
[n, m] = size(A);

piv = A(k,k);
newpivot = k;

for i = k+1:n            % check all rows below pivot
   if abs(A(i,k)) &gt; abs(piv) 
       newpivot = i;
       piv = A(i,k);
   end
end

if piv == 0              % trouble - can't avoid 0 pivot
   err = 1;
elseif newpivot ~= k     % best pivot elsewhere, swap rows
   for j = k:m
     temp1 = A(k,j);
     temp2 = A(newpivot,j);
     A(k,j) = temp2;
     A(newpivot,j) = temp1;
   end
end


</PRE></BODY></HTML>

⌨️ 快捷键说明

复制代码 Ctrl + C

搜索代码 Ctrl + F

全屏模式 F11

切换主题 Ctrl + Shift + D

显示快捷键 ?

增大字号 Ctrl + =

减小字号 Ctrl + -