Model pro WinPP:

# obycejne diferencialni rovnice - rovnice samocisteni jezer

# model peti kanadskych jezer

# Horejsi jezero a jeho pritok

x1' = -f13/V1*x1 + p1*f1

# Michignaske jezero

x2' = -f23/V2*x2 + p2*f2

# Huronske jezero - sem tece voda z predchozich dvou jezer

x3' = f13/V1*x1 + f23/V2*x2 - f34/V3*x3 + p3*f3

# Erijske jezero - sem tece voda z Huronskeho jezera

x4' = f34/V3*x3 - f45/V4*x4 + p4*f4

# Ontario

x5' = f45/V4*x4 - f56/V5*x5 +p5*f5

#pocatecni podminky (v procentech pocatecniho znecisteni)

x1(0)=100

x2(0)=0

x3(0)=0

x4(0)=0

x5(0)=0

# nastaveni parametru - objemy jezer

par V1=2900, V2=1180, V3=850, V4=116, V5=393

# nastaveni parametru - koncentrace necistot v pritocich jezer

# do vsech jezer priteka cista voda, jenom do tretiho jezera priteka

# p3*f3=0.075 necistot za rok (mereno v procentech necistot, ktere jsou

# na pocatku v Horejsim jezere)

par p1=0, p2=0, p3=0.005, p4=0, p5=0

# nastaveni parametru - pritoky jezer

par f1=15, f2=38, f3=15, f4=17, f5=14

# nastaveni parametru - toky mezi jezerama

par f13=15, f23=38, f34=68, f45=85, f56=99

# parametry zobrazovani

# bude nas zajimat Horejsi jezero - x1

@ bound=1000, total=1000, yp=x1, xlo=0,ylo=0,xhi=1000,yhi=100

# vyplavovani necistot probiha velice dlouhou dobu,

# proto nastavime hrubsi krok integrace

@ dt=0.5

done

Průběh znečištění v Hořejším jezeře.

Průběh znečištění v ostatních jezerech. Všímejte si toho, kde dosáhne znečištění nejmenší hodnoty (Erijske jezero) a jak to souvisí s objemem jezera a s rychlostí, jakou z jezera vytéká voda.

.