Хотелось бы вообще узнать правильный ли код и в правильном направлении я иду,так сказать

void GetSolutionTridiagSLAE(double* a, double* b, double* c, double* bb, int N, double* x)
{
 
    // 
    for (int n = 1; n < N; n++)
    {
        a[n] = a[n-1]*a[n] - b[n-1]*c[n-1];
        if (n < N-1)
            b[n] = b[n]*a[n-1];
        bb[n] = bb[n]*a[n-1] - bb[n-1]*c[n-1];
    }
    x[N-1] = bb[N-1]/a[N-1];
    // 
    for (int n = N-2; n > -1; n--)
        x[n] = (bb[n] - b[n]*x[n+1])/a[n];
}
 
void MatrixVectorProduct(double **A , double* x, double *b, int N)
{
    for (int i = 0; i < N; i++)
    {
        b[i] = 0;
        for (int j = 0; j < N; j++)
            b[i] = b[i] + A[i][j]*x[j];
    }
}
 
int main()
{
    int N = 50;  
    double *a = new double[N];
    double *b = new double[N-1];
    double *c = new double[N-1];
    double *bb = new double[N];
    double *x = new double[N];
    
    // 
    for (int i = 0; i < N; i++)
    {
        a[i] = (double)i*2;
        bb[i] = (double)i*3; 
            if (i < N-1)
            {
                b[i] = (double)i*4;
                c[i] = (double)i*5;
            }
    }
    
    GetSolutionTridiagSLAE(a, b, c, bb, N, x);
    // 
 
 
    
    // 
    double *testbb = new double[N];
    // 
 
 
    MatrixVectorProduct(A, x, testbb);
    // 
 
    
    delete [] a ;
    delete [] b ;
    delete [] c ;
    delete [] x ;
    delete [] bb;
    delete [] A[i]; 
}