в чем ошибка ? не считает y[n]

#include <iostream>
 
using namespace std;
 
int main()
{
    const int n = 3;
    double sum = 0;
    double A[n][n], L[n][n], U[n][n], R[n][n];
    double f[n];
 
    //задаем матрицу A[][] ...
 
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cout << "\nA[" << i << "][" << j << "]: ";
            cin >> A[i][j];
 
            L[i][j] = 0;
            U[i][j] = 0;
 
            if (i == j)
                U[i][j] = 1;
        }
    }
    cout << endl;
 
    //==============================================
 
    //находим первый столбец L[][] и первую строку U[][]
 
    for (int i = 0; i < n; i++)
    {
        L[i][0] = A[i][0];
        U[0][i] = A[0][i] / L[0][0];
    }
 
    //дальше вычисляем L[][], U[][] по формуле
 
    for (int i = 1; i < n; i++)
    {
        for (int j = 1; j < n; j++)
        {
            if (i >= j) //нижний треугольник
            {
                sum = 0;
                for (int k = 0; k < j; k++)
                    sum += L[i][k] * U[k][j];
 
                L[i][j] = A[i][j] - sum;
            }
            else // верхний
            {
                sum = 0;
                for (int k = 0; k < i; k++)
                    sum += L[i][k] * U[k][j];
 
                U[i][j] = (A[i][j] - sum) / L[i][i];
            }
        }
    }
 
    for (int i = 0; i < n; i++)
    {
        cout << "\nf[" << i << "]: ";
        cin >> f[i];
    }
    cout << endl;
 
    //====================================================
    cout << "\n\nFirst matrix:" << endl;
    for (int i = 0; i < n; i++) // вывод матрицы А
    {
        for (int j = 0; j < n; j++)
            cout << "  " << A[i][j] << "   ";
        cout << "\n\n";
    }
 
    cout << "\n\nL matrix:" << endl; 
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
            cout << "  " << L[i][j] << "   ";
        cout << "\n\n";
    }
 
    cout << "\n\nU matrix: " << endl; 
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
            cout << "  " << U[i][j] << "   ";
        cout << "\n\n";
    }
 
    cout << "\n\nf vector: " << endl;
    for (int j = 0; j < n; j++)
    {
        cout << "  " << f[j] << endl;
    }
    cout << endl;
 
    //====================================================
    
    double y[n];    //Находим y[n] по формуле (2.11)
    for (int i = 0; i < n; i++) {
        sum = 0;
        for (int k = 0; k < i; k++) {
            sum += L[i][k] * y[k];
            y[i] = f[i] - sum;
        }
        
        cout << y[n];
    }
    cout << endl;
 
 
    system("pause");
    return 0;
}

cout << y[i] << endl;

double y[n];    //Находим y[n] по формуле (2.11)
 
 
    for (int i = 0; i < n; i++) {
        sum = 0;
        for (int k = 0; k < i; k++) {
            sum += L[i][k] * y[k];
    //      y[i] = f[i] - sum;
        }
        y[i] = (f[i] - sum)/L[i][i];
 
        cout << y[i] << "   ";
    }