MATLAB: Please help with the code for Gauss Seidel Method

gauss-seidel

Please help with the code , I do not know what I am doing wrong. If you can please give me an advice. the first approximation should be x1=-0.2; x2=0.156 and x3=-.508. I have attached the code.

Please follow the format that I have already , because I will use it to do a bigger matrix 33 by 9

Best Answer

Your assigment to T should be within for loops where both sibscripts are changing

It is uncommon to call norm() passing in a single value. Also those subscripts imply the assigmnent to norm2 should be inside for loops.

Related Solutions

MATLAB: Update values in while loop and store values from a while loop in arrays

% Exercise 2.1 Order of Convergence
clear
close
clc
% Define Variables
F = 12;
x_old = 0;
itr = 0;
epsilon_new = 1e30;
epsilon_target = 1e-6;
while epsilon_new > epsilon_target
    itr = itr+1;
   
    x_new = springcalc(F, x_old);
    epsilon_new = abs(x_new - x_old);
    epsilon(itr)=epsilon_new;
    x_old = x_new;  
end
semilogy(1:length(epsilon), epsilon);
function [x_new] = springcalc(F, x_old)
k_0 = 30;
x_new = F/(k_0*(1+x_old^(0.1)));
end

MATLAB: Place array in a matrix

Perhaps I'm misunderstanding but if

 x = [1.0000         0
      0.5000    1.0000
      0.5000    1.5000
           0   -1.0000
     -1.0000    1.0000
      0.5000   -1.0000];

which is a 6x2 array, then wouldn't

 x*x'

be the desired 6x6 array?

Best Answer

Related Solutions

MATLAB: Update values in while loop and store values from a while loop in arrays

MATLAB: Place array in a matrix

Related Question