Реализация класса полином

#ifndef CLASS_POLINOM_CPP
#define CLASS_POLINOM_CPP
class Polynom
{
private:
    unsigned int deg; //bazmandami astichany
    double *coef; //cucich gortsakitsneri zangvatsi vra
    void CorrectDeg(); 
public:
    Polynom(); //aranc parametreri constructor
    Polynom(unsigned int, double*); //parametrerov
    Polynom(const Polynom &); //copy constructor
    ~Polynom(); //destructor
 
    unsigned int GetDeg(); //astichan stanalu function
    double GetCoef(unsigned int); //i-rd astichani bazmandami gortsakocneri stacum
    unsigned int SetCoef(double, unsigned int); //i-rd astichani bazmandami gortsakci taly
    Polynom operator + (const Polynom &);   
    Polynom operator = (const Polynom &);   
    friend Polynom MultConst(double, Polynom &); //bazmandamneri bazmapatkum
    void InputPolynom();    //neratselu function
    void OutputPolynom();   //artatselu function
};
#endif
 
#include <iostream>
using namespace std;
 
// 0 astichani ev 0 gortsakicnerov
Polynom::Polynom()
{
    deg=0;
    coef=new double[deg+1];
    coef[0]=0.0;
}
 
Polynom::Polynom(unsigned int new_deg, double *new_coef)
{
    deg=0;
    for(unsigned int i=0;i <new_deg+1;i++)
        if(new_coef[i]!=0) 
deg=i;   
    coef=new double[deg+1];
    for(unsigned int i=0;i<=deg;i++)
        coef[i]=new_coef[i]; 
}
 
Polynom::Polynom (const Polynom &f)
{
    deg=f.deg;
    coef=new double[deg+1];
    for(unsigned int i=0;i<=deg;i++)
        coef[i]=f.coef[i];
}   
 
Polynom::~Polynom()
{
    delete[] coef;
}
unsigned int Polynom::GetDeg()
{
    return deg;
}
 
double Polynom::GetCoef(unsigned int i)
{
    if(i<=deg) 
        return coef[i];
    else 
        return 0.0;
}
 
unsigned int Polynom::SetCoef(double new_coef, unsigned int i)
{
    if(i<=deg) 
coef[i]=new_coef;
    else 
return deg;
    CorrectDeg();   
    return deg;
}
 
Polynom Polynom::operator + (const Polynom &t)
{
    Polynom *result;
 
    if(deg>=t.deg)
{ 
        result=new Polynom(deg,coef);
        for(unsigned i=0;i<=t.deg;i++)
            result->coef[i]=result->coef[i]+t.coef[i];
    }
    else
{
        result=new Polynom(t.deg,t.coef);
        for(unsigned i=0;i<=deg;i++)
            result->coef[i]=result->coef[i]+coef[i];
    }
    result->CorrectDeg();
    return *result;
}
 
Polynom Polynom::operator = (const Polynom &t)
{
    deg=t.deg;
    delete[] coef;
    coef=new double[deg+1];
    for(unsigned int i=0;i<=deg;i++)
        coef[i]=t.coef[i];
    return *this;
}
 
//bazmandamy tvov bazmapatkeleu function
Polynom MultConst(double k, Polynom &t)
{
    if(k==0)
   {
        Polynom result;
        return result;
    }
    else
    {
        int deg=t.deg;
        double *tmp_coef=new double[deg+1];
        for(int i=0;i<=deg;i++)
            tmp_coef[i]=k*t.coef[i];
        Polynom result(deg,tmp_coef);
        delete[] tmp_coef;
        return result;
    }
}
/*astichani chshgrtum 
;orinak ete deg=4,bayc dra gortsakicy 0 e,ays funkcian deg-y kdardzni 3 ev aydpes minchev vor gortsakicy chlini 0 */
void Polynom::CorrectDeg()
{
    if(coef[deg]==0)
       while(deg && coef[deg]==0)
        deg--;
 }
 
 
void Polynom::InputPolynom()
{
    cout << "Input degree: ";
    cin >> deg;
    delete[] coef;
    coef=new double[deg+1];
    for(unsigned int i=0;i<deg;i++)
        cin >> coef[i];        
}
 
void Polynom::OutputPolynom()
{
    if(coef[deg]==1)
        cout << "x^" << deg;
    else 
 
    if(coef[deg]==-1)
        cout << "-x^" << deg;
    else
        cout << coef[deg] << "x^" << deg; 
    for(int i=deg-1;i>0;i--)
{
        if(coef[i]>0)
       {
            if(coef[i]==1)
                cout << " + " << "x^" << i;
            else
                cout << " + " << coef[i] << "x^" << i;
        }
else 
if(coef[i]<0)
            if(coef[i]==-1)
                cout << " - " << "x^" << i;
            else
                cout << " - " << (-1)*coef[i] << "X^" << i;
    }
 
    if(coef[0]>0)
        cout << " + " << coef[0] << endl;
    else 
        if(coef[0]<0)
             cout << " - " << (-1)*coef[0] << endl;
}
 
 
 
//main
#include <iostream>
using namespace std;
#include "Class_Polinom.h"
int main()
{
    Polynom a;
    a.InputPolynom();
    a.OutputPolynom();
    Polynom b;
    b.InputPolynom();
    b.OutputPolynom();
    Polynom sum;
    sum=a+b;
    sum.OutputPolynom();
    return 0;
}