(*
VL: Rechnernutzung in der Physik, WS11/12
gram.m

Berechne Gram-Determinante einer Matrix
*)

(* ***

<<kap_5_gram.m

mat//MatrixForm

Timing[Det[mat]]

*** *)

dim = 10;
mat = {};
For[i=1,i<=dim,i++,
    row = {};
    For[j=1,j<=dim,j++,
	row = Join[row,{v[i,j]}];
    ];
    mat = Join[mat,{row}];
];

(v[i_,j_]/;i>j) := v[j,i];