Code covered by the BSD License  

Highlights from
Method of elimination of Gauss with pivoting partial

Method of elimination of Gauss with pivoting partial

by

 

23 Oct 2006 (Updated )

Solve Linear Equation in format Ax=b

test.m
%  Progetto Metodo di Gauss con Pivoting Parziale
%
%             Programma elaborato da
%
%      Giovanni DI CECCA & Virginia BELLINO
%           50 / 887           408 / 466
%
%             http://www.dicecca.net


% Funzione di Test

% Lo scopo di questa funzione  quello di andare ad eseguire 
% tutti i test previsti nell'esercizio, in modo automatico 
% ed in successione.
%
% Per evitare possibili errori nell'uso delle matrici, 
% ad ogni test viene eseguito un reset della memoria, 
% sia video (il clc) sia fisica (il clear)

clc % pulisci schermo

clear % Pulisci memoria

disp(' Progetto Metodo di Gauss con Pivoting Parziale')
disp(' ')
disp('             Programma elaborato da')
disp(' ')
disp('      Giovanni DI CECCA & Virginia BELLINO')
disp('            50 / 887         408 / 466')
disp(' ')
disp('             http://www.dicecca.net')
disp(' ')
disp(' ')

disp('test - 1')

a=[2 0 -2 0
    9 10 -4 -1
    3 4 -2 1
    1 -2 -2 3]

b=[0
    -5
    -3
    1]

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 2')

a=[3 10 9
    8 4 -10
    -18 8 48]

b=[1
    0
    -2]

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 3')

a=[10 1 2 3
    1 2 0 0
    2 0 6 1
    3 0 1 5]
    
b=[17
    -6
    12
    31]

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 4')

a=[7 -2 3
    1 11 3
    3 12 15]
    
b=[8
    35
    42]

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 5.1')

q=5;

% Ciclo di creazione di A
for t=1: q
    for z=1: q
        a(t,z)=(t)^(z-1);
    end
end

% Ciclo di creazione di B
for t=1: q
    for z=1: q
        b(t,1)=a(t,q)/2;
    end
end


[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 5.2')

q=10;

% Ciclo di creazione di A
for t=1: q
    for z=1: q
        a(t,z)=(t)^(z-1);
    end
end

% Ciclo di creazione di B
for t=1: q
    for z=1: q
        b(t,1)=a(t,q)/2;
    end
end


[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 5.3')

q=20;

% Ciclo di creazione di A
for t=1: q
    for z=1: q
        a(t,z)=(t)^(z-1);
    end
end

% Ciclo di creazione di B
for t=1: q
    for z=1: q
        b(t,1)=a(t,q)/2;
    end
end

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 5.4')

q=50;

% Ciclo di creazione di A
for t=1: q
    for z=1: q
        a(t,z)=(t)^(z-1);
    end
end

% Ciclo di creazione di B
for t=1: q
    for z=1: q
        b(t,1)=a(t,q)/2;
    end
end

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 6')

a=[10E-15 3
    2 3]

b=[3
    5]

    
[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc % pulisci schermo

clear % Pulisci memoria

disp('test - 7')

a=[1 1
    1 1.001]

b=[1
    0]

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)

disp('Premere un tasto per continuare')
pause

%---------------------------------------------------------------------

clc

clear

disp('test - 7.1')

a=[1 (1+4/1000)
    1 1.001]

b=[1
    0]

[A,B,pivot,det,iflag,x]=myGaussSolve2(a,b)


Contact us