[C++] Pętla for wzorst znacznika

#include <iostream>
#include <fstream>
#include <cmath>
using namespace std;

int main()
{
    //**** Stale ****
    double lambda = 0.75;
    double c = 3000000;
    double cb = 3996200;
    double Gb = 0.00125;
    double Tb = 37;
    double k = 4995.25;
    double Qmet = 245;
    double L = 0.035;
    double delT = 5;
    int n = 10;
    double Q;
    //**** Obliczenie warunku stabilności ****
    double WS;
    double a, h;
    Q = k * Tb + Qmet;
    a = lambda / c;
    h = L / n;
    WS = 1 -(( 2 * a * delT ) /( h * h ) ) -(( k * delT ) / c );
    if( WS >= 0 )
         cout << "Warunek stabilnosci zostal spelniony" << endl;
   
    else
         cout << "Nie spelniono warunku" << endl;
    //**** Fourier ****
    long double const1 = 1 -(( 2 * a * delT ) /( h * h ) ) -(( k * delT ) / c );
    long double const2 =( a * delT ) /( h * h );
    long double const3 =( Q * delT ) / c;
    long double T[ 9 ][ 21 ];
    int I, F; // z capslockiem zeby sie odroznialy od licznikow dla ciepla
   
    for( I = 1; I <= 7; I++ ) // warunek poczatkowy - w kazdym wezle (I) dla zerowej chwili czasu (F=0) -> T=36,6
   
    {
        T[ 0 ][ 0 ] = 37.5;
        T[ 8 ][ 0 ] = 32;
        T[ I ][ 0 ] = 36.6;
    }
   
   
    for( F = 1; F <= 20; F++ )
    {
        for( I = 1; I <= 8; I++ )
        {
           
            T[ 0 ][ F ] = 37.5;
            T[ I ][ F ] =( T[ I - 1 ][ F - 1 ] + T[ I + 1 ][ F - 1 ] ) / 2;
            T[ 8 ][ F ] = 32;
           
           
        }
       
    }
   
    for( F = 0; F <= 5; F++ )
    {
        cout << "Rozklad temperatur F[" << F << "] " << endl;
        for( I = 0; I <= 8; I++ )
        {
            cout << T[ I ][ F ] << endl;
        }
       
    }
    system( "pause" );
    return 0;
}