prog2

double det_hilbert(int n) {
	double det_h = 1.0;
	double det_b = 1.0;
	for (int i = 1; i < n - 1; i++){
		det_h *= fact(i);
	}
	det_h = pow(det_h, 4);
	for (int j = 1; j < 2*n - 1; j++){
		det_b *= fact(j);
	}
	double det = det_h / det_b;
	return det;
}