#ifndef _CORPOCELESTE
#define _CORPOCELESTE

#include "vettore.h"
#include <vector>

const double G = 6.672e-11;

using namespace std;

class CorpoCeleste {
 public:
  CorpoCeleste();
  CorpoCeleste(double mass);
  CorpoCeleste(double mass, vettore x0);
  CorpoCeleste(double mass, vettore x0, vettore v0);
  void setMass(double mass) { M = mass; };
  void setX0(vettore x0) { x = x0; };
  void setV0(vettore v0) { v = v0; };
  vettore X() { return x; };
  vettore V() { return v; };
  double m() { return M; };
  vettore F(vector<CorpoCeleste*> list);  // restituisce la forza agente sul 
                                          // corpo prodotta dall'insieme list
  virtual void evolve(double dT, vettore F);
  virtual void evolve(double dT, vector<CorpoCeleste*> list);
 protected:
  vettore x;
  vettore v;
  double M;
};

#endif