Add method positional argument

Added compulsory positional argument "method", determining whether to use the Jacobi or Gauss-Seidel iterative methods. This variable is then passed onto the solve function, which determines which method to use.