File: invhilb-1.m

package info (click to toggle)
octave2.1 1%3A2.1.73-13
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 37,028 kB
  • ctags: 20,874
  • sloc: cpp: 106,508; fortran: 46,978; ansic: 5,720; sh: 4,800; makefile: 3,186; yacc: 3,132; lex: 2,892; lisp: 1,715; perl: 778; awk: 174; exp: 134
file content (8 lines) | stat: -rw-r--r-- 257 bytes parent folder | download | duplicates (9) 
1
2
3
4
5
6
7
8
 
result4 = [16, -120, 240, -140;
	   -120, 1200, -2700, 1680;
	   240, -2700, 6480, -4200;
	   -140, 1680, -4200, 2800];

(invhilb (1) == 1 && invhilb (2) == [4, -6; -6, 12]
 && invhilb (4) == result4
 && abs (invhilb (7) * hilb (7) - eye (7)) < sqrt (eps))