Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement