I am trying to solve a matrix system to find the relative arrival rates of a queueing network using Gauss-Seidel.The maple commands are below:
restart;
with(Student[NumericalAnalysis]); with(LinearAlgebra);
A := Matrix([[1, -.333, -.333, -.333], [0, 1, -.333, -.333], [0, -.333, 1, -.333], [0, -.333, -.333, 1]]);
IsDefinite(A, 'query' = 'positive_semidefinite');
true
b := Vector([1, 1, 1, 1]);
IterativeApproximate(A, initialapprox = Vector([1, 1, 1, 1]), tolerance = 10^(-3), maxiterations = 20, stoppingcriterion = relative(infinity), method = gaussseidel);
Error, (in Student:-NumericalAnalysis:-IterativeApproximate) check that the augmented matrix has the correct dimensions
I do not understand this error as the matrix is 4x4 as shown. Can anyone see where I went wrong?