load matrix.dat
a = spconvert( matrix )
load rhs.dat
sol = a \ rhs
save solution.dat sol -ascii -double
exit
