Jim+Williams 源代码

function A=add(A) %ADD %Add multiple of one row of a matrix to another row. %Used in row reduction. %Calling format: A=add(A)…

function P=div(P) %DIV %Divide a row by a matrix by nonzero constant. %Used in row reduction. %Calling format: A=div(A) %…

function u = gangle(X,Y,A) %GANGLE %The angle between two vectors X and Y, in degrees, in Rn. %Rn has inner product defined by …

function X=cramer(A,B) %CRAMER %Solution to system AX=B using Cramer's Rule. %Calling format: cramer(A,B) %Copyright Gare…

function pow(P) %POW %Computes powers of a matrix. %Option of rational number format. %Calling format: pow(A) %Copyright Ga…

function n = mag(X) %MAG %The magnitude of a vector X. %Calling format: mag(X) %Copyright Gareth Williams, Stetson Universit…

function cn(P) %CN %Condition number of matrix based on 1-norm. % cond(A)= ||A||*||inv(A)|| %Also gives related concep…

function n = gnorm(X,A) %GNORM %The magnitude of a vector X in Rn. %Rn has 'inner product' defined by a symmetric matrix A. …

function C=cross(A,B) %CROSS %Cross product of two vectors. %The vectors must have 3 elements each. %Calling format: cross(A,…

function Q=perminv(P) %PERMINV %Inverse of a permutation defined by a matrix A. %Calling format: B=perminv(A) %Copyright Gar…