Function solveGS

Solve Ax = b iteratively using Gauss-Seidel method

void solveGS(T, ulong N) (
  Matrix!(T,N) a,
  ref Vector!(T,N) x,
  Vector!(T,N) b,
  uint iterations = 10
);