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