class derivative{
  Vector3D dp;      // change in pos
  Vector3D dv;      // change in vel
  
  derivative(){
    dp = new Vector3D(0,0,0);
    dv = new Vector3D(0,0,0);
  }
}