PROC SQRTA
SET NODOECHO NOOUTPUT NOECHO
* [AMATRIX]=INPUT MATRIX
* [XMATRIX]=OUTPUT MATRIX
* GET SQRT OF A MATRIX
MATRIX A_=[AMATRIX]
MATRIX V_=EIGVEC(A_)
MATRIX D_=EIGVAL(A_)
MATRIX Q_=SQRT(D_)
MATRIX X_=V_*DIAG(Q_)*V_'
MATRIX [XMATRIX]=X_
DELETE / ALL_
SET DOECHO OUTPUT ECHO
PROCEND