PRB . problem   = 8    ; % select Godunov first-order upwind scheme
PRB . init      = 1    ; % select square wave for initial data
PRB . Ncells    = 100  ; % select 100 cell inteval subdivision
PRB . CFL       = 0.49 ; % advance with DT do not excede 0.9 * CFL
PRB . a         = 0.25 ; % speed of the signal
PRB . StepTime  = 1    ; % save the frame every 0.5 second
PRB . FinalTime = 200  ; % compute the solution op to 3 second
% call the solver
[X,UFRAMES,TFRAMES] = LagenSolver(PRB) ;
% generate movie of the solution
LagenMovie(X,UFRAMES,TFRAMES) ;

% plot the solution al frame n. 4
%[,N] = size(TFRAMES) ; 
%plot( X, UFRAMES(N,:) );
