//Błędy maxymalne dla tmax w funkcji kroku przestrzennego
        if (wybor == 1)
        {
            // wyliczenie bledu maksymalnego w funkcji kroku
            double blad_max_KMB = 0.0, blad_max_laasonen = 0.0;
            for (int j = 0; j < liczba_wezl - 1; j++) // Pętla po węzłach
            {
                //KMB
                if (fabs(wynik_KMB[t_bezposrednie - 1][j] - wynik_anali[t_posrednie - 1][j]) > blad_max_KMB) {
                    blad_max_KMB = fabs(wynik_KMB[t_bezposrednie - 1][j] - wynik_anali[t_posrednie - 1][j]);
                }
                //Laasonen
                if (fabs(wynik_laasonen[t_posrednie - 1][j] - wynik_anali[t_posrednie - 1][j]) > blad_max_laasonen) {
                    blad_max_laasonen = fabs(wynik_laasonen[t_posrednie - 1][j] - wynik_anali[t_posrednie - 1][j]);
                }
            }
            fprintf(wykres1, "%.15lf,%.15lf,%.15lf\n", h, blad_max_KMB, blad_max_laasonen);
            cout << h << " " << blad_max_KMB << " " << blad_max_laasonen << endl;
        }