hilbert.mata

来自「是一个经济学管理应用软件 很难找的 但是经济学学生又必须用到」· MATA 代码 · 共 20 行

*! version 1.0.0  10sep2004
version 9.0

mata:

real matrix Hilbert(real scalar n) 
{
	real matrix H
	real scalar un, i, j

	un = trunc(n)
	H = J(un, un, .)
	for (i=1; i<=un; i++) {
		for (j=1; j<=un; j++) H[i,j] = 1/(i+j-1)
	}
	return(H)
}

end