load matrix.dat;
a = spconvert( matrix );
a = a + tril( transpose( a ), -1 );
b=size(a)
n = nnz(a)
load rhs.dat;
sol = a \ rhs;
save solution.dat sol -ascii -double;
