Solve the following using the Gauss-Seidel Method. [A] = [4 -2 1; 1 5 -3; 2 2 5]

Question

Explanation / Answer

coeff = [4 -2 1; 1 5 -3; 2 2 5]; constants = [11;-6;7]; tol = 0.001; n = length(constants); checktol = ones(n,1);%Used for Error term evaluation X = zeros(n,1); %Gauss-Siedel iterations iteration = 0; while max(checktol) > tol iteration = iteration + 1; Z = X; for i = 1:n j = 1:n; j(i) = []; Xtemp = X; Xtemp(i) = []; X(i) = (constants(i) - sum(coeff(i,j) * Xtemp)) / coeff(i,i); end Xsolution(:,iteration) = X; checktol = sqrt((X - Z).^2);%Recalculating the error term end solution = [X] ++++++++++++++++++++++++++++++++++++++++ A nice explanation on Gauss-siedel method using Matlab available from here. You can give a try. It's good article. https://ece.uwaterloo.ca/~dwharder/NumericalAnalysis/04LinearAlgebra/gauss/

Related Questions

Navigate

• Browse (All)
• Browse S
• Subjects

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.