Реализовать иерархию классов для работы с геометрическими телами

friend ostream & operator<<(ostream & os, const Ellipsoid & obj)
{
os << "minor axis a, b, c" << obj.a << "," << obj.b << "," << obj.c << " " << "Volume: " << obj.Volume() << " " << "Area: " << obj.Area() <<"\n";
return os;
}

 #include <iostream>
#include <cmath>
#include <cstdlib>
#define pi 3.1416
#define p 1.6075
using namespace std;
 
class Figure
{
public:
    virtual double Area() const {}
    virtual double Volume() const {}
    virtual void Increment(double N, const Figure & obj) {}
    virtual ~Figure() {}
};
 
class Ellipsoid: public Figure
{
    double a, b, c;
public:
    Ellipsoid(double a, double b, double c);
    ~Ellipsoid() {}
    double Area() const;
    double Volume() const;
    void Increment(double N, const Ellipsoid & obj);
    friend ostream& operator<<(ostream & os, const Ellipsoid & obj);
};
 
class Ball: public Figure
{
    double R;
public:
    Ball(double R);
    ~Ball() {}
    double Area() const;
    double Volume() const;
    void Increment(double N, const Ball & obj);
    friend ostream& operator<<(ostream & os, const Ball & obj);
};
 
class Cylinder: public Figure
{
    double R, H;
public:
    Cylinder(double R, double H);
    ~Cylinder() {}
    double Area() const;
    double Volume() const;
    void Increment(double N, const Cylinder & obj);
    friend ostream& operator<<(std::ostream & os, const Cylinder & obj);
};
 
class Cone: public Figure
{
    double R, H;
public:
    Cone(double R, double H);
    ~Cone() {}
    double Area() const;
    double Volume() const;
    void Increment(double N, const Cone & obj);
    friend ostream& operator<<(std::ostream & os, const Cone & obj);
};
 
Ellipsoid::Ellipsoid(double a, double b, double c): a(a), b(b), c(c)
{
    if (a <= 0.0 || b <= 0.0 || c <= 0.0 || a == b == c) {exit(1);}
}
 
double Ellipsoid::Area() const
{
    double S = 4.0*pi*pow((1.0/3.0)*(pow(a * b, p) + pow(a * c, p) + pow(b * c, p)), 1.0/p);
    return S;
}
 
double Ellipsoid::Volume() const 
{
    double V = 4.0*pi*a*b*c/3.0;
    return V;
}
 
void Ellipsoid::Increment(double N, const Ellipsoid & obj)
{
    this -> a = N*obj.a;
    this -> b = N*obj.b;
    this -> c = N*obj.c;
}
 
 
ostream& operator<<(ostream & os, const Ellipsoid & obj)
 {
    os << "minor axis a, b, c" << obj.a << "," << obj.b << "," << obj.c << " " << "Volume: " << obj.Volume() << " " << "Area: " << obj.Area() <<"\n";
    return os;
 }
  
 Ball::Ball(double R): R(R)
{
    if (R <= 0.0) exit(1);
}
 
double Ball::Area() const
{
    double S = 4.0*pi*R*R;
    return S;
}
 
double Ball::Volume() const 
{
    double V = 4.0*pi*R*R*R/3.0;
    return V;
}
 
void Ball::Increment(double N, const Ball & obj)
{
    this -> R = N*obj.R;
}
 
 
ostream& operator<<(ostream & os, const Ball & obj)
 {
    return os << "Radius of ball" << obj.R << " " << "Volume: " << obj.Volume() << " " << "Area: " << obj.Area() <<"\n";
 }
  
 Cylinder::Cylinder(double R, double H): R(R), H(H)
{
    if (H <= 0.0 || R <= 0.0) exit(1);
}
 
double Cylinder::Area() const
{
    double S = 2.0*pi*(R*R + R*H);
    return S;
}
 
double Cylinder::Volume() const 
{
    double V = pi*R*R*H;
    return V;
}
 
void Cylinder::Increment(double N, const Cylinder & obj)
{
    this -> R = N*obj.R;
    this -> H = N*obj.H;
}
 
 
ostream& operator<<(ostream & os, const Cylinder & obj)
 {
    return os << "Radius and Height of cylinder" << obj.R << "," << obj.H << "Volume: " << obj.Volume() << " " << "Area: " << obj.Area() <<"\n";
 }
  
 Cone::Cone(double R, double H): R(R), H(H)
{
    if ( R <= 0.0 || H <= 0.0) exit(1);
}
 
double Cone::Area() const
{
    double S = pi*R*(R + sqrt(R*R + H*H));
    return S;
}
 
double Cone::Volume() const 
{
    double V = pi*R*R*H/3.0;
    return V;
}
 
void Cone::Increment(double N, const Cone & obj)
{
    this -> R = N*obj.R;
    this -> H = N*obj.H;
}
 
 
ostream& operator<<(ostream & os, const Cone & obj)
 {
    os << "Radius and Height of cone" << obj.R << "," << obj.H << "Volume: " << obj.Volume() << " " << "Area: " << obj.Area() <<"\n";
    return os;
 }
  
 int main () 
 {
    double N = 2.0;
    Figure *f[4];
    f[0] = new Cone(1.0, 1.0);
    cout << &f[0] << "," << f[0] << *f[0] << "\n";
    system("pause");
 }