Очистка памяти при работе с функциями

double* Function(T_GSM_element* a , double *F, int N)
{
double *X; //вектор решения
double *X_first=new double[N];         //начальное приближение
for (int i = 0; i < N; i++) {X_first[i]=0;} // (в дальнейшем будет являться предыдущим решением)
double *s,*r_old,*r_new,a1,b1; // вспомогательные переменные при решении
r_old=r_new=s=vector_summ(Mult_Matrix_and_Vecror (a,X_first,N), F, N, 0);
 
for (int i=0;i<10000; i++)
{
a1=scalar_product(..)/scalar_product( ..);//не буду писать довольно громоздкий код
//функция scalar_product - return double              
X=vector_summ(...);
//функция  vector_summ - double *d=new double[count]...;return d;
X_first=X;
r_new=vector_summ(...);
b1=scalar_product(..)/scalar_product(..);
r_old=r_new;
s=vector_summ(...);
}
return X;
}