przewodzenie ciepła w płyctce 2d

#include <iostream>
#include<string>
#include<vector>
#include<ctime>
#include<cstdlib>
#include<fstream>
using namespace std;

const int wymiar = 560; // rozmiar tablicy
const int ileRazy = 50; // ilość iteracji

struct punkt2D {
    int x;
    int y;
    punkt2D( int _x, int _y )
        : x( _x )
         , y( _y )
    { }
};

int main() {
    clock_t start = clock();
    srand( time( NULL ) );
    int kier; // kierunek poruszania się
    double suma; //suma do średniej
    ////////////////////////////
    //inicjalizacja  i wypełnienie
    double ** tab2 = new double *[ wymiar ]; //alokacja pamieci
    for( int i = 0; i < wymiar; ++i ) {
        tab2[ i ] = new double[ wymiar ]; //alokacja pamieci
        for( int j = 0; j < wymiar; ++j ) { //wpisanie wartosci do tablicy
            tab2[ i ][ j ] = 0;
        }
    }
    for( int i = 0; i < wymiar; ++i, cout << endl ) {
        tab2[ i ][ wymiar - 1 ] = 100;
    }
   
   
    ////////////////////////////
    // tablica do zapisu wyników
    double ** ntab = new double *[ wymiar ]; //alokacja pamieci
    for( int i = 0; i < wymiar; ++i ) {
        ntab[ i ] = new double[ wymiar ]; //alokacja pamieci
        for( int j = 0; j < wymiar; ++j ) //wpisanie wartosci do tablicy
             ntab[ i ][ j ] = 0;
       
    }
    ///////////////////////////
    //liczenie metodą montecarlo
    double aaa = ileRazy;
    double srednia;
   
    for( int q = 1; q < wymiar - 1; q++ ) { //pętla po wierszach
        for( int c = 1; c < wymiar - 1; c++ ) { //pętla po kolumnach
            for( int i = ileRazy; i > 0; i-- ) { //ilość proóbek
                punkt2D punkt( c, q ); // punkt startowy
                //pętla losująca drogę
                while(( punkt.x > 0 ) &&( punkt.x < wymiar - 1 ) &&( punkt.y > 0 ) &&( punkt.y < wymiar - 1 ) )
                {
                   
                    kier =( rand() % 4 ) + 1;
                    if( kier == 1 ) {
                        punkt.y++;
                    }
                    else if( kier == 2 ) {
                        punkt.y--;
                    }
                    else if( kier == 3 ) {
                        punkt.x++;
                    }
                    else if( kier == 4 ) {
                       
                        punkt.x--;
                    }
                }
                suma += tab2[ punkt.y ][ punkt.x ]; //suma próbek
            }
            ntab[ c ][ q ] =( suma / aaa ); //zapis średniej do talicy
            suma = 0;
            cout << "czas wykonania: " << clock() - start << " msek" << endl; //czyszczenie zmiennej
        }
    }
    /*for(int i=0;i<wymiar;i++){
    for(int j=0;j<wymiar;j++){
    cout <<ntab[i][j]<<" ";
    }
    cout<< endl;
    }*/
    cout << "czas wykonania: " << clock() - start << " msek";
    ///////////////////////////////
    //zapis do pliku
    fstream zapis;
    zapis.open( "wynik_montecarlo1.txt", ios::out );
    /*for(int i=0;i<wymiar;i++){
    for(int j=0;j<wymiar;j++){
    zapis <<ntab[i][j]<<", ";
    }
    zapis << endl;
    }*/
    zapis << " czas obliczen: " <<( clock() - start ) / 60000 << endl;
    zapis.close();
   
    //czyszczenie pamięci
    for( int i( 0 ); i < wymiar; ++i )
         delete[] tab2[ i ]; //uwolnienie pamieci
   
    delete[] tab2; //uwolnienie pamieci
    tab2 = NULL;
   
    for( int i( 0 ); i < wymiar; ++i )
         delete[] ntab[ i ]; //uwolnienie pamieci
   
    delete[] ntab; //uwolnienie pamieci
    tab2 = NULL;
   
   
    return 0;
}
/ cpp ]