x = [ 1:10 ] ; N = 10 ; 
v = zeros(size(x) ) ;
im = zeros(size(x) ) ;
im(2:5) = 1.0/9.0 * [1.1 1.0 0.9 1.0 ] ;
im(6:9) = 1.0     * [1.1 1.0 0.9 1.0 ] ;

v(2:5) = 2.0  ; 
v(3) = 2.21 ; 
v(4) = 2.621 ; 
v(6:9) = -0.2 ; 
v(7) = -0.27 ;
v(8) = -0.37 ;
## to make a nice brief movie use this
Tc = 0.030123091;  Number_samples_required = 1000 ;  doplot = 1 ; 

bounce1;

## see RUNbounce1.m for a run that gives nice histograms

