Потеря памяти в программе

double* direct_multiple(double angle1, double angle2)
{
    double* new_matrixS = new double[4];// пока что с двойкой, потом двойка будет входить в переменную для весов
    
    new_matrixS[0]=2*cos(angle1)*cos(angle2);
    new_matrixS[1]=2*cos(angle1)*sin(angle2);
    new_matrixS[2]=2*sin(angle1)*cos(angle2);
    new_matrixS[3]=2*sin(angle1)*sin(angle2);
    return new_matrixS;
}
 
 
 
double* funkg(double koord[],double mu, double lambda)//,double Z, double fi)
{
    
    double* g = new double[2];
    g[0] = /*1 + koord[0] - koord[0]*koord[1];//*/1 + pow(koord[0],2) - ((4*mu + 2*lambda)/(mu + lambda))*koord[0]*koord[1];
    g[1] = /*1 + koord[0] - koord[0]*koord[1];//*/1 + pow(koord[1],2) - ((4*mu + 2*lambda)/(mu + lambda))*koord[0]*koord[1];
    return g;
}
 
double* get_B(double alpha,double fi, double R, double Rc[], double Z, double y[])
{
    double thetta;
    thetta = atan2(y[1]-Rc[1],y[0]- Rc[0]);
    double* matr = new double[4];   
 
    double* Ef = new double[4];
    double* matrfi = new double[4];
    double* matrthfi = new double[4];
    
    Ef[0]=1;
    Ef[1]=0;
    Ef[2]=0;
    Ef[3]=1;
    
    double p0 = 0;
 
    p0 = (alpha*Z - R)/R;
 
    double a1; // случайная величина от 0 до 1 используется для вычисления матрицы B
    double a2; // случайная величина от 0 до 1 используется для вычисления матрицы B если a1 > alpha/(2+alpha)
 
    a1 = rand();
    a2 = rand();
 
    if (thetta < 0) {thetta = thetta + 2*pi;}
 
    matrfi = direct_multiple(fi,fi);
    matrthfi = direct_multiple((thetta + fi)/2,(thetta + fi)/2);
    
    if (a1 < alpha/(2+alpha))
    {
        matr = matrfi;
    }
    else if (2*R - alpha*Z > 0) 
    {
        if (a2 < (2*R - alpha*Z)/2*R)
        {
            matr = Ef;
        }
        else 
        {
            matr = matrthfi;
        }
    }
    else if (a2 < (alpha*Z - 2*R)/(2*R*p0))
    {
        matr[0] = -1*p0;
        matr[1] = 0;
        matr[2] = 0;
        matr[3] = -1*p0;
    }
    else
    {
        matr[0] = matrthfi[0]*p0;
        matr[1] = matrthfi[1]*p0;
        matr[2] = matrthfi[2]*p0;
        matr[3] = matrthfi[3]*p0;
    }
    //delete[] matrfi;
    //delete[] matrthfi;
    return matr;
}

double* matrix = new double[4];

matrix = get_B(...);

double* matr = new double[4];
//...
matr = matrfi;