1.0        L      ! dlka intervalu, tj. x in [0,L]
0.2        T      ! doba een, tj. t in [0,T]
500        N      ! poet bunk
20         M      ! M+1 je poet asovch vrstev pro vstup vsledk
0e-2       eps    ! koeficient difze
3          met    ! metoda 
0.9        CFL    ! CFL konstanta 
1          IE     ! metoda asov diskretizace
1          PP     ! typ poaten podmnky 
0.2        CL     ! pro PP = 1: poloha lev pepaky
0.4        CP     ! pro PP = 1: poloha prav pepky
2          uL     ! pro PP = 1: hodnota vlevo od CL, tj. pro x < CL
5          uS     ! pro PP = 1: hodnota uprosted, tj. pro  CL < x < CP
1          uP     ! pro PP = 1: hodnota vpravo od CP, tj. pro x : CP

Pozn8mky.

met volme takto: 1 ... Godunov
                  2 ... Van Leer
                  3 ... Engquist-Osher
                  4 ... Rusakov

Zadme-li M, pak dostaneme een v M+1 asovch hladinch
pslunch asm i*T/M, i=0,1,...,M.

Je-li eps > 0, pak je mon vliv difznho lenu zapotat buto
uitm explicitn Eulerovy metody, to volme IE=0, nebo pomoc
implicitn Eulerovy metody, v tom ppad volme IE=1.

CFL konstantu je teba volit <= 1, doporuen volba je CFL = 0.9

PP uruje, jak se v m-funkci u0.m nadefinuje esen v ase t = 0.

Tak nap. pro PP=1, CL=CP, eps=0 je v u0.m  Riemannova uloha, kdy 
vlevo od zvolen polohy C prepky je u = uL a vpravo je u = uR. 
Zvolme-li uL > uP, dostaneme rzovou vlnu, kter se pohybuje 
rychlosti s = 0.5*(uL + uP). Zvolme-li uL < uP, dostaneme 
spojitou vlnu zedn, kdy se elo vlny pohybuje rychlost
uP a konec vlny rychlost uL.

Pro PP=1, CL<CP, eps=0 mme v ase t=0 ti konstantn stavy:
u=uL pro x<CL,
u=uS pro CL<x<CP,
u=uP pro CP<x.
Postupn vzniknou dv vlny,
lev vlna:  pro uL<uS vlna zedn, elo (x-CL)/t=uL, tl (x-CL)/t=uS
            pro uL>uS rzov vlna, (x-CL)/t=0.5(uL+uS) 
prav vlna: pro uS<uP vlna zedn, elo (x-CP)/t=uS, tl (x-CP)/t=uP
            pro uS>uP rzov vlna, (x-CP)/t=0.5(uS+uP) 
Pro CL=0.2,CP=0.4, uL=2,uS=5, uP=1 dostihne tl vlny zedn 
rzovou vlnu v case, kter spoteme takto:

t=(x-0.4)/5, kde x je een rovnice (x-0.2)/5=(x-0.4)/3.

Dostaneme x=0.7 a odtud t=0.1. Pak se pika vlny sniuje, 
pro dostaten velk L,T (teba L=5,T=2) dostaneme rzovou vlnu
spojujc lev stav uL=2 s pravm stavem uP=1, pohybujc se 
rychlost 1.5.