%% TP Noté 23%%
close all;
a=0.1;
c=0.5;
T=2;
xmin = -pi;
xmax = pi;
dx =0.1;
dt = 0.01
CFL= c*dt/dx^2
u=difrea(a, c, T, xmin, xmax, dx, dt);
%%
% les 2 schémas donnent des résultats comparables et corrects 
% on augmente légèrement le pas de temps
dt =0.0102
CFL= c*dt/dx^2
u=difrea(a, c, T, xmin, xmax, dx, dt);
% le schéma explicite donne des résultats aberrants
% les schémas implicites donnent des résultats corrects.
% la CFL dépasse le seuil de 1/2,
% cela explique l'instabilité du schéma explicite.