% nettoyage memoire 
clear all
%% compactifie l'affichage 
format short e
format compact 
%% s'arreter dans la fonction qui produit une erreur 
%dbstop if error
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% parametres 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% dimension 1 
d=1;
%% fixons le pas en espace 
dx=0.05;
%% le pas en temps 
dt=dx; 
%% la taille du domaine en espace 
xmin=-1 ; 
xmax=1 ;
%% La grille de difference finies 
xx=[xmin:dx:xmax] ; 
%% sx est la taile de la grille 
[a,sx]=size(xx) ;
%% l'horizon en temps la grille en temps 
%% et sa taille ...
tmax= 2; 
tspan = [0:dt:tmax]; 
[a,st]= size(tspan) ; 
%% allocation de la variable de sortie (phase) 
%% pas reellement necessaire 
psi=zeros(st,sx); 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% initialisation  (gradpsi0 pas utile ici)
[psi0,grpsi0]=fi0(xx) ;
%% resolution du systeme d'edo issus de la semi 
%% discretisation en espace.
%% on passe en argument le pas et la grille
%% l'argument flag=[] n'est pas utilise   
[t psi]=ode23('semidi',tspan,psi0,[],dx,xx) ;
%[t psi]=monode('semidi',tspan,psi0,[],dx,xx) ;
%% la solution est dans psi. 






