Question: Solving a Linear System with the SOR method

I have a coefficient matrix A and the constant matrix b.

I want to find the unknowns x_1,x_2,...,x_8. When I put my code and run it in Maple I get an error. Any ideas as to why this is happening? The code I tried is as follows:

 

 

with(Student[NumericalAnalysis]):

A := Matrix([[-1, 0, 0, (1/2)*sqrt(2), 1, 0, 0, 0], [0, -1, 0, (1/2)*sqrt(2), 0, 0, 0, 0], [0, 0, -1, 0, 0, 0, 1/2, 0], [0, 0, 0, -(1/2)*sqrt(2), 0, -1, -1/2, 0], [0, 0, 0, 0, -1, 0, 0, 1], [0, 0, 0, 0, 0, 1, 0, 0], [0, 0, 0, -(1/2)*sqrt(2), 0, 0, (1/2)*sqrt(3), 0], [0, 0, 0, 0, 0, 0, -(1/2)*sqrt(3), -1]]);

b := Vector([0, 0, 0, 0, 0, 10000, 0, 0]);

LinearSolve(A, b, maxiterations = 300, method = SOR(1.25), tolerance = 10^(-2));

 

I also tried it with a SOR variable of 1.5 instead of 1.25, but that gives me the same error.

Please Wait...