Jim+Williams 源代码

function P=mult(P) %MULT %Mutiply a row by a matrix by a constant. %Used in row reduction. %Calling format A=mult(A) %Cop…

function d=dot(X,Y) %DOT %The dot product of two vectors X and Y. %Calling format: dot(X,Y) %Copyright Gareth Williams, Ste…

function R=siml(P,Q) %SIM %Performs similarity transf inv(B)*A*B. %Calling format: siml(A,B) %Copyright Gareth Williams, Ste…

function pic(P,w) %PIC %subroutine for drawing graphs in GJpic and Gausspic. %Calling format: pic(A,w) %Copyright Gareth Wil…

function p=ginner(X,Y,A) %GINNER %The inner product of two vectors X and Y in Rn. %Rn has 'inner product' defined by a symmet…

function square(n) %SQUARE %Draws a fractal square. %Calling format: square(n), where n is desired # iterations. %Copyright…

function markov(P,X) %MARKOV %Markov Chain Model. %Transition matrix A and initial state column vector X. %Gives distribution…

function B = adj(A) %ADJ % The adjoint of a matrix A. % Calling format: B = adj(A) %Copyright Gareth Williams, Stetson Univ…

function d = gdist(X,Y,A) %GDIST %The distance between points X and Y in Rn. %Rn has inner product defined by a symmetric mat…

function U = tdot(A,B,C) %TDOT %Triple scalar product of three vectors. %The vectors must all have 3 elements. %Calling forma…