x = [ 0   14 15 16  17 18 19  20 21     35 ] ; N = 10 ; 
## extra displacement to add to the particles so they appear to bounce
## off each other
xtweak = [ 0 : (N-1) ] * 5 ; 
ytweak = 6;

## in cradle2 the masses are 4,9,1,1,1,1,1,1,1
## and the initial situation has all the energy in the left two masses
v = zeros(size(x) ) ;
im = zeros(size(x) ) ;
im(2:9) = 1.0 ;
im(2) = 1.0/4.123 ;
im(3) = 1.0/9.0123 ;

v(2:3) = - 2.0  ; 

## this specifies how often in simulation time we make a plot
Tc = 0.2030123091;  Number_samples_required = 30000 ;  doplot = 1 ; 
## this specifies how often in real time the plots will appear
time_between_plots = 0.025;

bounce2;

