ash/amuchastegui   nov 1968    subroutine snorml
inversion of coefficient matrix and solution of double precision
normal equations with extended precision intermediate steps
using fact that coefficient matrix is symmetric

solve normal equations and determine inverse of coefficient
matrix with double precision arithmetic using fact that
matrix is symmetric
gauss-jordan method with pivots on diagonal

*  start=4000
calculate standard deviations and correlations,
printout correlation matrix, read normal equations
back into storage