% test_konsistenz.m
%----------------------------------------
% Testfunktion für die Konsistenzordnung
%    verschiedener Diskretisierungen
%       der 1D Poisson-Gleichung
%----------------------------------------

clear all
format long

disp('******************************')
disp(' Vergleich: Exakte Lösung vs. ')
disp('  2. Ordnung Diskretisierung  ')
disp('******************************')

fehler = 1;

for i = 2:10
  
    N = 2^i;
    fehleralt = fehler;

    [uex] = uexakt(N);
    [u2o] = pois1d_2o(N);

    fehler = norm(uex-u2o,Inf);
    
    if ( i > 2 )
       N
       Faktor = fehleralt/fehler
    end   
    
end   

disp('  ')
disp('******************************')
disp('  ')
disp('  ')


clear all
format long

disp('******************************')
disp(' Vergleich: Exakte Lösung vs. ')
disp('  4. Ordnung Diskretisierung  ')
disp('******************************')

fehler = 1;

for i = 2:10
  
    N = 2^i;
    fehleralt = fehler;

    [uex] = uexakt(N);
    [u4o] = pois1d_4o(N);

    fehler = norm(uex-u4o,Inf);
    
    if ( i > 2 )
       N
       Faktor = fehleralt/fehler
    end   
    
end   

disp('  ')
disp('******************************')
disp('  ')
disp('  ')


clear all
format long

disp('******************************')
disp(' Vergleich: Exakte Lösung vs. ')
disp('     MSV Diskretisierung      ')
disp('******************************')

fehler = 1;

for i = 2:10
  
    N = 2^i;
    fehleralt = fehler;

    [uex] = uexakt(N);
    [uMSV] = pois1d_MSV(N);

    fehler = norm(uex-uMSV,Inf);
    
    if ( i > 2 )
       N
       Faktor = fehleralt/fehler
    end   
    
end   

disp('  ')
disp('******************************')
disp('  ')
disp('  ')