double trace(M,n)/* return trace of double matrix */register double *M;register int n;{	register int i,con;	register double sum;	sum=0;	con=n+1;	for(i=0;i<n;i++,M+=con)		sum+= *M;	return(sum);}