subroutine prnsum to print summary for multiple parameter sets
or multiple epochs
p. macneil august, 1976
maximum of 50 solutions

use of b(i,j) i=1,ncol (=2*ntotal+2), j=1,nparam
b(i,j) =
xnom(j), i = 1
scale(j), i = 2
adjustment(j), i = 3,5,7,....
sigma(j), i = 4,6,8,.....