## This script file uses the function planet10

###############################################
## Set the initial conditions 
###############################################
# Duration of simulation, and time-step; and initial time 
 T = 1000 ;  dt = 1 ; t = 0 ; 
################################################
# Initial condition for position and velocity
################################################
 x = [ 93 , 0 ] ;
# x has dimensions [L]
 v = [0,1.1] ;
# v has dimensions [L/T]
################################################
# setting this to a value greater than zero 
#   causes plots to occur, and defines the duration of pause while plotting (in seconds)
plotduring=0.015;
velocityfactor=10; 
  ## velocityfactor has dimensions [T] - it specifies how to render the velocity
  ##                                  (with a vector of size v*velocityfactor


[ x , v , t , logbook ] = planet10( x , v , t , T , dt , plotduring , velocityfactor );