Переделать код с console wizard на application

#include <vcl.h>
 
#include <iostream.h>
#include <iomanip.h>
 
 
using namespace std;
 
 
int J, N;
double L0, L9, E, E1, X0, X1;
double** A;
double** B;
double* X;
double* F;
double* L;
 
void MATR(int N, double **&A);
void DET(int N, double **&A, double &S);
void TAB(int N, double *&X, double *&F);
void CPN(int N, double *&X, double *&F);
void URAV(double X0, double &X, double E, double E1);
void RandomMatrix(double **&M, int n);
void OutputMatrix(double **&M, int n);
 
double PN(double X1);
 
int main()
{
 
    cout << "Enter size of matrix( >3 or =3) :" << endl;
 
    N = 2;
 
    while (N < 3)
    {
        cin >> N;
    }
    N++;
 
    A = new double *[N];
    for (int i = 0; i < N; i++)
        A[i] = new double[N];
 
    B = new double *[N];
    for (int i = 0; i < N; i++)
        B[i] = new double[N];
    X = new double[N];
    F = new double[N];
    L = new double[N];
 
 
    cout << "Enter boundary eigenvalue (L0):" << endl;
    cin >> L0;
    cout << "Enter boundary eigenvalue (L9):" << endl;
    cin >> L9;
    cout << "Enter first error (E) :" << endl;
    cin >> E;
    cout << "Enter second error (E1) :" << endl;
    cin >> E1;
 
    RandomMatrix(B, N);
 
    cout << "Output Matrix :" << endl;
    OutputMatrix(B, N);
 
    TAB(N, X, F);
    CPN(N, X, F);
    X0 = L0;
 
    for (J = 0; J < N - 1; J++)
    {
        X1 = L9;
        URAV(X0, X1, E, E1);
        L[J] = X1;
        cout << " L [" << J << "] = " << X1 << endl;
    }
 
    system("pause");
}
 
void TAB(int N, double *&X, double *&F)
{
    int I, J, M;
    double H, R, S;
 
    H = (L9 - L0) / N;
    R = L0;
 
    for (M = 0; M < N ; M++)
    {
        for (I = 1; I < N; I++)
        {
            for (J = 1; J < N; J++)
            {
                if (I == J)
                {
                    A[I][I] = B[I][I] - R;
                }
                else
                {
                    A[I][J] = B[I][J];
                }
            }
        }
        DET(N, A, S);
        X[M] = R;
        F[M] = S;
        R += H;
    }
}
 
double PN(double X1)
{
    int I;
    double P;
    P = F[N-1];
 
    for (I = N - 2; I > -1; I--)
    {
        P = F[I] + (X1 - X[I]) * P;
    }
 
    for (I = 0; I < J; I++)
    {
        P = P / (L[I] - X1);
    }
    return P;
}
 
void DET(int N, double **&A, double &S)
{
    int I, J, K, K1;
    double P, R;
 
    P = 1.0;
    for (K = 0; K < N - 1; K++)
    {
        K1 = K + 1;
        S = A[K][K];
        J = K;
        for (I = K1; I < N; I++)
        {
            R = A[I][K];
            if (abs(R) > abs(S))
            {
                S = R;
                J = I;
            }
        }
 
        if (S == 0) return;
 
        if (J != K)
        {
            P = -P;
            for (I = K; I < N; I++)
            {
                R = A[K][I];
                A[K][I] = A[J][I];
                A[J][I] = R;
            }
        }
 
        for (J = K1; J < N; J++) A[K][J] = A[K][J] / S;
 
        for (I = K1; I < N; I++)
        {
            R = A[I][K];
            for (J = K1; J < N; J++) A[I][J] = A[I][J] - A[K][J] * R;
        }
        P = P * S;
    }
    S = P * A[N - 1][N - 1];
}
 
void CPN(int N, double *&X, double *&F)
{
    int I, J;
    double A, B;
    for (J = 1; J < N; J++)
    {
        A = F[J - 1];
        B = X[J - 1];
        for (I = J; I < N; I++) F[I] = (A - F[I]) / (B - X[I]);
    }
}
 
void URAV(double X0, double &X, double E, double E1)
{
    double D, Y, R;
    R = X - X0;
    D = PN(X0);
 
    do
    {
        Y = PN(X);
        if (abs(Y) < E1) return;
        R = R / (D - Y) * Y;
        D = Y;
        X = X + R;
 
    } while (abs(R)< E);
}
 
void RandomMatrix(double **&M, int n)
{
 
    for (int i = 1; i < n; i++)
    {
        for (int j = 1; j < n; j++)
        {
            M[i][j] = (double)(rand() % 100 + 1);
        }
    }
 
}
void OutputMatrix(double **&M, int n)
{
    for (int i = 1; i < n; i++)
    {
        for (int j = 1; j < n; j++)
        {
            cout << setw(5)<< left << M[i][j];
        }
        cout << endl;
    }
}