#include <vector>
#include <CLHEP/Matrix/Vector.h>

using namespace std;

#include "swim.h"

#define PI 3.1415926

int main() 
{
  double charge = 1.0, p = 1.0 , theta = PI/6.0, phi = PI/4.0;
  vector<HepVector*> traj;
  swim(charge, p, theta, phi, traj);
}