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