Zad2a.m

n = [5 10 50 100 200]';
eps1 = zeros( length(n), 1 );
T = zeros( length(n), 1 );

for k=1:length(n)
    
    % Generowanie macierzy A
    % A(i,j) = 12 dla i=j
    % A(i,j) = 3.8 dla i=j-1 lub i=j+1
    % A(i,j) = 0 dla pozostalych
    A = zeros( n(k) );

    for i=1:n(k)
        A(i,i) = 12;
    end

    for j=2:n(k)
        A(j-1,j) = 3.8;
    end

    for j=1:n(k)-1
        A(j+1,j) = 3.8;
    end
    


    % Generowanie wektora b
    % b(i) = 4.5 - 0.5*i
    b = zeros( n(k), 1 );

    for i=1:n(k)
        b(i) = 4.5 - 0.5*i;
    end
    
    
    
    tic;
    x = GJ( A, b, n(k) );
    T(k) = toc;
    
    
    
    eps1(k) = norm( A*x - b );
end



T = T .* 1000; % aby wyniki byly w ms 



% Wykres bledu rozwiazania od liczby rownan
figure(1);
hold on;

plot(n,eps1,'r.', 'MarkerSize', 25);
grid on;

title('Zależnosc bledu rozwiazania od liczby rownan (czesc a))');
xlabel('Liczba rownan n');
ylabel('Blad rozwiazania eps1');

hold off;



% Wykres czasu dzialania algorytmu od liczby rownan
figure(2);
hold on;

plot(n,T,'r.', 'MarkerSize', 25);
grid on;

title('Zaleznosc czasu dzialania algorytmu od liczby rownan (czesc a))');
xlabel('Liczba rownan n');
ylabel('Czas działania algorytmu [ms]');

hold off;