Программа выдает не верный результат

#include <math.h>
#include <fstream>
#include <iostream>
#include <conio.h>
using namespace std;
 
double ** CreateMatrix (int M, int N, ifstream *pFin = NULL) 
    {
    double **A = new double *[M];
    double *Body = new double [M * N], *pRow = Body;
    for (int i = 0; i < M; A[i++] = pRow, pRow += N);
    if (pFin == NULL) return A;
    for (int i = 0; i < M; i++)
        for (int j = 0; j < N; *pFin >> A[i][j++]);
    return A;
    }
 
int main (int n, char * narg[]) 
    {
    setlocale (LC_ALL, "Russian"); 
    if (n != 3) 
        {
        cout << "Неверно указаны параметры:" << endl; 
        cout << "1) Путь до входного файла" << endl << "2) Путь выводного файла";
        getch(); 
        return 1;
        }
    ifstream Fin (narg[1]);
    if (!Fin.is_open())
        {
        cout << narg[1] << "  не открылся";
        getch();
        return 2;
        }
    ofstream Fout (narg[2]);
    if (!Fout.is_open())
        {
        cout << narg[2] << "  не открылся";
        getch();
        return 2;
        }
    int N;
    Fin >> N;
    double **Ain = CreateMatrix (N, N, &Fin); 
    double **A = CreateMatrix (N, N);
    for (int i = 0; i < N; i++)
            for (int j = 0; j < N; j++)
                A[i][j] = Ain[i][j];
    double *V = new double [N];
    double *Vold = new double [N];
    double D;
    Vold[0] = 0.0;
    for (;;)
        { //9
        double S = 0;
 
        for (int i = 0; i < N; i++) //10
            {
            V[i] = 0;
            for (int j = 0; j < N;j++)
                V[i]+= A[i][j];
            S += V[i];
            }
        for (int i = 0; i < N; i++) //11
            V[i]=(V[i] / S);
 
        if (Vold[0] != 0.0) //12
        {
            double t = 0;
            for (int i = 0; i < N; i++)
                t += pow((V[i]-Vold[i]), 2); //12.1
            D = sqrt (t / N);
            if (D < 0.0000001) //12.2
                break;
        }
 
        for (int i = 0; i < N; i++) //13
            Vold[i] = V[i];
 
        for ( int i = 0; i < N; i++) //14
            {
            for ( int j = 0; j < N; j++)
                {
                double t = 0;
                for ( int k = 0; k < N; k++)
                    t += (A[i][k] * A[k][j]);
                A[i][j] = t;
                }
            }
        }//15
 
    double *E = new double [N]; //16
    for (int i = 0; i < N; i++)
        {
        E[i] = 0;
        for (int j = 0; j < N; j++)
            E[i] += Ain[i][j]*V[j];
        }
    for (int i = 0; i < N; i++)
        E[i] = E[i] / V[i];
    
    Fout << "Главный собственный вектор" << endl;
    for (int i = 0; i < N; i++)
        Fout << V[i] << " ";
    Fout << endl << "Оценки главного собственного значения" << endl;
    for (int i = 0; i < N; i++)
        Fout << E[i] << " ";
    }