Найти наибольшее по абсолютной величине собственное значение

#include <iostream>
#include <fstream>
#include <math.h>
#include <conio.h> 
using namespace std;
double Method(double **A, int n)
{
    double *Y;  //начальный вектор
    Y = new double [n];
    double *R;  //итерации (нахождение собственных векторов)
    R = new double [n];
    double *Q;  //Y(k - 2)
    Q = new double [n];
    cout<<"Начальный вектор:"<<endl;
    for(int j = 0; j < n; j++)
    {
        Y[j] = 1;   
        cout<<Y[j]<<endl;
    }
    getch();
    
    for(int i = 0; i < n; i++)
    {
        R[i] = 1;
    }
 
    double L1 = 0, L2 = 0, L3 = 0; //старые и новые собственные значения
    do
    {
        L1 = L2;
        for(int i = 0; i < n; i++)
        {
            Q[i] = Y[i];
            Y[i] = R[i];
            R[i] = 0;
        }
        // умножение матрицы на вектор
        for(int i = 0; i < n; i++)
        {
            for(int j = 0; j < n; j++)
                R[i]+=A[i][j]*Y[j];
        }
        //новые собственные значения
 
        double D = 0; //скалярное произведение векторов:нового на старый
        for(int i = 0; i < n; i++)
            D+=R[i]*Y[i];
 
        double D1 = 0; //скалярное произведение старого вектора на себя самого
        for(int j = 0; j < n; j++)
            D1+=Y[j]*Y[j];
        L2 = D/D1;
    }
    while(abs(L2 - L1) > 0.0001);
    cout<<"Первое собственное значение:"<<endl;
    cout<<L2<<endl;
 
    double max, mod;
    for(int i = 0; i < n; i++)
    {
        if(R[i] < 0)
            mod = -R[i];
        else
            mod = R[i];
    }
        max = R[0];
        for(int i = 1; i < n; i++)
        {
            if(max < R[i])
                max = R[i];
        }
    cout<<"Собственный вектор:"<<endl;
    double *vekt;
    vekt = new double [n];
    for(int i = 0; i < n; i++)
    vekt[i] = R[i]/max;
    for(int i = 0; i < n; i++)
    printf("%f\n", vekt[i]);
 
    for(int i = 0; i < n; i++)
    {
        R[i] = 1;
        Y[i] = 1;
    }
    L1 = 0;
    int count = 0; // счетчик
    bool flag = true; //старые знаки
    bool flag1 = true; //новые знаки
 
    do
    {
        for(int i = 0; i < n; i++)
        {
            Q[i] = Y[i];
            Y[i] = R[i];
            R[i] = 0;
        }       
 
        for(int i = 0; i < n; i++)
        {
            for(int j = 0; j < n; j++)
                R[i]+=A[i][j]*Y[j];
        }
 
        double S = 0; //скалярное произведение векторов:нового на старый
        for(int i = 0; i < n; i++)
            S+=R[i]*Y[i];
 
        double S1 = 0; //скалярное произведение старого вектора на себя самого
        for(int j = 0; j < n; j++)
            S1+=Y[j]*Y[j];
 
        double S2 = 0; //скалярное произведение векторов:старого на предыдущий  
        for(int i = 0; i < n; i++)
            S2+=Y[i]*Q[i];
        L3 = L1;
        L1 = (S - L2*S1)/(S1 - L2*S2);  
 
        /*Для нахождения второго собственного значения (так как оно находится неверно) рассматриваем последовательности. 
        Начнем сравнивать со второй итерации: для этого заведем переменную - счетчик итераций.
        Находим значение и сравниваем с первоначальным значением,заносим информацию о знаке в переменную. 
        Добавляем проверку на знак: если у нас не первая итерация и знак предыдущий не равен новому полученному,
        то выходим и берем предыдущее собственного значение,
        иначе сохраняем информацию о новом знаке в переменную,которую завели под информацию о старом знаке, и продолжаем.*/
        if(L1 > L3)
            flag = true;
        else
            flag = false;
        if(count!=0 && flag1!=flag)
            break;
        else
            flag = flag1;
        count++;
    }
        while(abs(L1 - L3) > 0.0001);
        cout<<"Второе собственное значение:"<<endl;
        cout<<L3;
        getch();
    return L3;
}
int main()
{
    setlocale(LC_ALL,"");
    int i, j;
    int n;
    double **A;
    ifstream in;
    in.open("1.txt", ios::in);
    in>>n;
    if(!in)
    {
        cout<<"Error"<<endl;
        return 1;
    }
    A = new double *[n];
    for(int i = 0; i < n; i++)
        A[i] = new double [n];
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            in>>A[i][j];
        }
    }
    //проверка на неправельный ввод
    if (!in.eof())
    {
        cout<<"Неверный ввод данных"<<endl;
        getch();
        return 0;
    }
    cout<<"Матрица имеет вид:"<<endl; // вывод матрицы на экран
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            cout<<A[i][j]<<" ";
        }
        cout<<endl;
    }
    getch();
    in.close();
    Method(A, n);
    cout<<endl;
    return 0;
}