| This is an update about location of MATLAB files sent to me from Bill
Muth:
> You can find the following MATLAB files in PERCPT::SYS$PUBLIC:[MISC]
>
> MATLAB.SAV - VMS Executables for MATLAB
> MATLAB.ZIP - MS-DOS executables
> MATSRC.ZIP - Source code for compiling under MS-DOS
> MATSRC.SAV - VMS, UNIX sources
>
> The VMS MATLAB.EXE requires the logical $MATLAB to point to the
> directory containing HELP.LIS.
>
> Bill
in MATSRC.sav there is source code for fortran programs used by matlab,
this is outline of some of the routines in lib.for, trying this version
of MATLAB , it does not have a lot of functions that i have on my PC-DOS
version of MATLAB that i bought few days ago from MathWorks in Mass., but
again the above version of MATLAB seem to the 1982 one, very old one.
You Can get a professional version of MATLAB for about $300, if you are
associated with a school, and a student version of it for $50 (the one
i got).
/Nasser
----------------------------------------------------------------------
SUBROUTINE WGECO(AR,AI,LDA,N,IPVT,RCOND,ZR,ZI)
C
C WGECO FACTORS A DOUBLE-COMPLEX MATRIX BY GAUSSIAN ELIMINATION
C AND ESTIMATES THE CONDITION OF THE MATRIX.
C
-------------------------------------------------------------------------
SUBROUTINE WGEFA(AR,AI,LDA,N,IPVT,INFO)
C
C WGEFA FACTORS A DOUBLE-COMPLEX MATRIX BY GAUSSIAN ELIMINATION.
C
-------------------------------------------------------------------------
SUBROUTINE WGESL(AR,AI,LDA,N,IPVT,BR,BI,JOB)
C
C WGESL SOLVES THE DOUBLE-COMPLEX SYSTEM
C A * X = B OR CTRANS(A) * X = B
C USING THE FACTORS COMPUTED BY WGECO OR WGEFA.
C
-------------------------------------------------------------------------
SUBROUTINE WGEDI(AR,AI,LDA,N,IPVT,DETR,DETI,WORKR,WORKI,JOB)
C
C WGEDI COMPUTES THE DETERMINANT AND INVERSE OF A MATRIX
C USING THE FACTORS COMPUTED BY WGECO OR WGEFA.
C
-------------------------------------------------------------------------
SUBROUTINE HTRIDI(NM,N,AR,AI,D,E,E2,TAU)
C THIS SUBROUTINE IS A TRANSLATION OF A COMPLEX ANALOGUE OF
C THE ALGOL PROCEDURE TRED1, NUM. MATH. 11, 181-195(1968)
C BY MARTIN, REINSCH, AND WILKINSON.
C HANDBOOK FOR AUTO. COMP., VOL.II-LINEAR ALGEBRA, 212-226(1971).
C
C THIS SUBROUTINE REDUCES A COMPLEX HERMITIAN MATRIX
C TO A REAL SYMMETRIC TRIDIAGONAL MATRIX USING
C UNITARY SIMILARITY TRANSFORMATIONS.
C
-------------------------------------------------------------------------
SUBROUTINE HTRIBK(NM,N,AR,AI,TAU,M,ZR,ZI)
C
C THIS SUBROUTINE IS A TRANSLATION OF A COMPLEX ANALOGUE OF
C THE ALGOL PROCEDURE TRBAK1, NUM. MATH. 11, 181-195(1968)
C BY MARTIN, REINSCH, AND WILKINSON.
C HANDBOOK FOR AUTO. COMP., VOL.II-LINEAR ALGEBRA, 212-226(1971).
C
C THIS SUBROUTINE FORMS THE EIGENVECTORS OF A COMPLEX HERMITIAN
C MATRIX BY BACK TRANSFORMING THOSE OF THE CORRESPONDING
C REAL SYMMETRIC TRIDIAGONAL MATRIX DETERMINED BY HTRIDI.
C
-------------------------------------------------------------------------
SUBROUTINE IMTQL2(NM,N,D,E,Z,IERR,JOB)
C
C THIS SUBROUTINE IS A TRANSLATION OF THE ALGOL PROCEDURE IMTQL2,
C NUM. MATH. 12, 377-383(1968) BY MARTIN AND WILKINSON,
C AS MODIFIED IN NUM. MATH. 15, 450(1970) BY DUBRULLE.
C HANDBOOK FOR AUTO. COMP., VOL.II-LINEAR ALGEBRA, 241-248(1971).
C
C THIS SUBROUTINE FINDS THE EIGENVALUES AND EIGENVECTORS
C OF A SYMMETRIC TRIDIAGONAL MATRIX BY THE IMPLICIT QL METHOD.
C THE EIGENVECTORS OF A FULL SYMMETRIC MATRIX CAN ALSO
C BE FOUND IF TRED2 HAS BEEN USED TO REDUCE THIS
C FULL MATRIX TO TRIDIAGONAL FORM.
C
-------------------------------------------------------------------------
SUBROUTINE CORTH(NM,N,LOW,IGH,AR,AI,ORTR,ORTI)
C
C THIS SUBROUTINE IS A TRANSLATION OF A COMPLEX ANALOGUE OF
C THE ALGOL PROCEDURE ORTHES, NUM. MATH. 12, 349-368(1968)
C BY MARTIN AND WILKINSON.
C HANDBOOK FOR AUTO. COMP., VOL.II-LINEAR ALGEBRA, 339-358(1971).
C
C GIVEN A COMPLEX GENERAL MATRIX, THIS SUBROUTINE
C REDUCES A SUBMATRIX SITUATED IN ROWS AND COLUMNS
C LOW THROUGH IGH TO UPPER HESSENBERG FORM BY
C UNITARY SIMILARITY TRANSFORMATIONS.
C
-------------------------------------------------------------------------
SUBROUTINE COMQR3(NM,N,LOW,IGH,ORTR,ORTI,HR,HI,WR,WI,ZR,ZI,IERR
* ,JOB)
C*****
C MODIFICATION OF EISPACK COMQR2 TO ADD JOB PARAMETER
C JOB = 0 OUTPUT H = SCHUR TRIANGULAR FORM, Z NOT USED
C = 1 OUTPUT H = SCHUR FORM, Z = UNITARY SIMILARITY
C = 2 SAME AS COMQR2
C = 3 OUTPUT H = HESSENBERG FORM, Z = UNITARY SIMILARITY
C ALSO ELIMINATE MACHEP
C C. MOLER, 11/22/78 AND 09/14/80
C OVERFLOW CONTROL IN EIGENVECTOR BACKSUBSTITUTION, 3/16/82
C*****
C
C THIS SUBROUTINE IS A TRANSLATION OF A UNITARY ANALOGUE OF THE
C ALGOL PROCEDURE COMLR2, NUM. MATH. 16, 181-204(1970) BY PETERS
C AND WILKINSON.
C HANDBOOK FOR AUTO. COMP., VOL.II-LINEAR ALGEBRA, 372-395(1971).
C THE UNITARY ANALOGUE SUBSTITUTES THE QR ALGORITHM OF FRANCIS
C (COMP. JOUR. 4, 332-345(1962)) FOR THE LR ALGORITHM.
C
C THIS SUBROUTINE FINDS THE EIGENVALUES AND EIGENVECTORS
C OF A COMPLEX UPPER HESSENBERG MATRIX BY THE QR
C METHOD. THE EIGENVECTORS OF A COMPLEX GENERAL MATRIX
C CAN ALSO BE FOUND IF CORTH HAS BEEN USED TO REDUCE
C THIS GENERAL MATRIX TO HESSENBERG FORM.
C
-------------------------------------------------------------------------
SUBROUTINE WSVDC(XR,XI,LDX,N,P,SR,SI,ER,EI,UR,UI,LDU,VR,VI,LDV,
* WORKR,WORKI,JOB,INFO)
C
C
C WSVDC IS A SUBROUTINE TO REDUCE A DOUBLE-COMPLEX NXP MATRIX X BY
C UNITARY TRANSFORMATIONS U AND V TO DIAGONAL FORM. THE
C DIAGONAL ELEMENTS S(I) ARE THE SINGULAR VALUES OF X. THE
C COLUMNS OF U ARE THE CORRESPONDING LEFT SINGULAR VECTORS,
C AND THE COLUMNS OF V THE RIGHT SINGULAR VECTORS.
C
-------------------------------------------------------------------------
SUBROUTINE WQRDC(XR,XI,LDX,N,P,QRAUXR,QRAUXI,JPVT,WORKR,WORKI,
* JOB)
C
C WQRDC USES HOUSEHOLDER TRANSFORMATIONS TO COMPUTE THE QR
C FACTORIZATION OF AN N BY P MATRIX X. COLUMN PIVOTING
C BASED ON THE 2-NORMS OF THE REDUCED COLUMNS MAY BE
C PERFORMED AT THE USERS OPTION.
C
-------------------------------------------------------------------------
SUBROUTINE WQRSL(XR,XI,LDX,N,K,QRAUXR,QRAUXI,YR,YI,QYR,QYI,QTYR,
* QTYI,BR,BI,RSDR,RSDI,XBR,XBI,JOB,INFO)
C
C WQRSL APPLIES THE OUTPUT OF WQRDC TO COMPUTE COORDINATE
C TRANSFORMATIONS, PROJECTIONS, AND LEAST SQUARES SOLUTIONS.
C FOR K .LE. MIN(N,P), LET XK BE THE MATRIX
C
C XK = (X(JPVT(1)),X(JPVT(2)), ... ,X(JPVT(K)))
C
C FORMED FROM COLUMNNS JPVT(1), ... ,JPVT(K) OF THE ORIGINAL
C N X P MATRIX X THAT WAS INPUT TO WQRDC (IF NO PIVOTING WAS
C DONE, XK CONSISTS OF THE FIRST K COLUMNS OF X IN THEIR
C ORIGINAL ORDER). WQRDC PRODUCES A FACTORED UNITARY MATRIX Q
C AND AN UPPER TRIANGULAR MATRIX R SUCH THAT
C
C XK = Q * (R)
C (0)
C
C THIS INFORMATION IS CONTAINED IN CODED FORM IN THE ARRAYS
C X AND QRAUX.
C
-------------------------------------------------------------------------
SUBROUTINE MAGIC(A,LDA,N)
C
C ALGORITHMS FOR MAGIC SQUARES TAKEN FROM
C MATHEMATICAL RECREATIONS AND ESSAYS, 12TH ED.,
C BY W. W. ROUSE BALL AND H. S. M. COXETER
C
-------------------------------------------------------------------------
DOUBLE PRECISION FUNCTION URAND(IY)
INTEGER IY
C
C URAND IS A UNIFORM RANDOM NUMBER GENERATOR BASED ON THEORY AND
C SUGGESTIONS GIVEN IN D.E. KNUTH (1969), VOL 2. THE INTEGER IY
C SHOULD BE INITIALIZED TO AN ARBITRARY INTEGER PRIOR TO THE FIRST CALL
C TO URAND. THE CALLING PROGRAM SHOULD NOT ALTER THE VALUE OF IY
C BETWEEN SUBSEQUENT CALLS TO URAND. VALUES OF URAND WILL BE RETURNED
C IN THE INTERVAL (0,1).
C
-------------------------------------------------------------------------
SUBROUTINE RROT(N,DX,INCX,DY,INCY,C,S)
C APPLIES A PLANE ROTATION.
-------------------------------------------------------------------------
|