# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.

# MW_MKIDMD -- Make the identity matrix.

procedure mw_mkidmd (ltm, ndim)

double	ltm[ndim,ndim]		#O set to the identity matrix
int	ndim			#I dimension of (square) matrix

int	i, j

begin
	do j = 1, ndim {
	    do i = 1, ndim
		ltm[i,j] = 0.0
	    ltm[j,j] = 1.0
	}
end