Освобождение динамически выделенной памяти - Visual C++ - Ответ 755244

double time_delay(VPath* route, double** Pst, int L, double* CLoad) // CLoad (Channel Loading) загрузка канала, L - длина пакета, Ts - время задержки, V - матрица пропускных способностей ЛЦТ-ов, Route - набор виртуальных путей;
{
    double pro = 0, pro2 = 0;
    int Vb = 64000;                                // скорость вокодера;
    double Tq = Td - (double)(L - Hip)/(double)Vb; // заданное время пребывания пакета в пути (исключая время пакетизации);
    double intensity;                              // intensity - интенсивность потока;
    double result = 0;
    factor* A;
    do{
        A = new factor[route->lenght + 1];   // вектор коэффициентов;
        for(int num_A = 0; num_A< route->lenght + 1; num_A++){
            A[num_A].arg = new double[route->lenght];
            A[num_A].max_degree = 0;
            A[num_A].denominator = 0;
            A[num_A].denominator_degree = 1;
            for(int j=0;j< route->lenght + 1; j++){
                if(j != 0) A[num_A].arg[j] = 0;
                else A[num_A].arg[j] = 1;
            }
        }
        for(int i=0; i < route->lenght; i++){            // i - номер вершины в пути
            intensity = ((double)V[route->knot[i]-1][route->knot[i+1]-1] / (double)(L + Hna))*(1 - CLoad[ElementNum (route->knot[i]-1, route->knot[i+1]-1, num_knots)]);
            for(int num_A = 0; num_A < route->lenght; num_A++){         // num_A - номер коэффициента
                if(num_A != i){                 // проверка умножения №1
                    if(A[num_A].denominator != intensity){  // проверка умножения №2
                        for(int k = A[num_A].max_degree + 1; k>=0; k--){          // k - степень аргумента S
                            if(k == 0)
                                A[num_A].arg[k] = A[num_A].arg[k] * intensity;
                            else{
                                if(k == A[num_A].max_degree)
                                    A[num_A].arg[k] = A[num_A].arg[k-1] + intensity;
                                else{
                                    if(k == A[num_A].max_degree + 1)
                                        A[num_A].arg[k] = 1;
                                    else
                                        A[num_A].arg[k] = A[num_A].arg[k] * intensity + A[num_A].arg[k-1];
                                }
                            }
                        }
                        A[num_A].max_degree++;
                    }
                    else
                        A[num_A].denominator_degree++;
                }
                else
                    A[num_A].denominator = intensity;
            }
            A[route->lenght].arg[0] = A[route->lenght].arg[0] * intensity;
        }
        for(int num_A = 0; num_A < route->lenght; num_A++){
            for(int s = 0; s < route->lenght; s++)
                if(s != num_A){
                    if (A[num_A].arg[num_A] == 0){
                        for (int i = 0; i < route->lenght + 1; i++)
                            A[i].arg[num_A] = A[i].arg[num_A] + A[i].arg[num_A+1];
                    }
                    intensity = A[num_A].arg[s]/A[num_A].arg[num_A];
                    for(int i = 0; i < route->lenght + 1; i++)
                        A[i].arg[s] = A[i].arg[s] - intensity * A[i].arg[num_A];
                }
        }
        for(int num_A = 0; num_A < route->lenght; num_A++){
            A[num_A].arg[num_A] = A[route->lenght].arg[num_A]/A[num_A].arg[num_A];
            if(A[num_A].arg[num_A] != 0){
                intensity = ((double)V[route->knot[num_A]-1][route->knot[num_A+1]-1] / (double)(L + Hna))*(1 - CLoad[ElementNum (route->knot[num_A]-1, route->knot[num_A+1]-1, num_knots)]);
                pro = factorial(A[num_A].denominator_degree - 1);
                pro2 = integral(A[num_A].denominator_degree - 1,Tq,intensity);
                result = result + A[num_A].arg[num_A] * Pst[route->knot[0]-1][route->knot[route->lenght]-1] / factorial(A[num_A].denominator_degree - 1) * integral(A[num_A].denominator_degree - 1,Tq,intensity);
            }
        }
        //for(int num_A = 0; num_A< route->lenght + 1; num_A++)
        //  A[num_A].~factor();
        delete[] A;
        A = (factor*)NULL;
        route = route->next_vpath;
    }while(route != (VPath*)NULL);
    return result = 1 - result;
}

struct factor
{
    double* arg;   // коэффициент при комплексном аргументе s ;
    int max_degree;  // старшая степень комплексного аргумента s;
    double denominator; // модуль корня знаменателя простейшей дроби;
    int denominator_degree; // степень знаменателя простейшей дроби;
    ~factor(){
        delete[] arg;
    }
};