Program FlexPDE w zastosowaniu do symulacji pól temperatury

1 = 0.4 { polowa szerokosci bloku }
c1 = 0.6 { polowa wysookosci bloku }
a2 = 0.25 { polowa grubosci wypelnienia }
b2 = 0.35 { polowa szerokosci wypelnienia }
T0 = 15 { temperatura poczatkowa }
Tpow = 25 { temperatura powietrza }
Tgleba = 2 { temperatura gleby }
hT = 10 { wspolczynnik wnikania ciepla }
V = VOL_INTEGRAL(1)
{ calkowita objetosc zbiornika wraz z
wypelnieniem }
Vz = VOL_INTEGRAL(1,'zbiornik') { objetosc zbiornika bez wypelnienia }
Vw = VOL_INTEGRAL(1,'wypelnienie') { objetosc wypelnienia }
INITIAL VALUES
Tp = T0
EQUATIONS
div(q) = 0
EXTRUSION
SURFACE 'denko zbiornika' z = -c1
SURFACE 'wieczko zbiornika' z = c1
BOUNDARIES
SURFACE 'denko zbiornika' VALUE(Tp) = Tgleba
SURFACE 'wieczko zbiornika' NATURAL(Tp) = hT*(Tp - Tpow)
REGION 1 'zbiornik' { plastik }
k = 0.22 { przewodność cieplna scianek }
START(0,0)
NATURAL(Tp) = 0 LINE TO (a1,0)
NATURAL(Tp) = hT*(Tp - Tpow) LINE TO (a1,b1)
LINE TO (0,b1)
NATURAL(Tp) = 0 LINE TO CLOSE
REGION 2 'wypelnienie' { styropian }
k = 0.033 { przewodność cieplna wypelnienia }