TITLE 'Elektro-termo-mechanika'
VARIABLES
Temp Volt
 DEFINITIONS
lambda ro
k=10
Vg=100 Va=0
Ta=273
j=-ro*grad(Volt)
q=-lambda*grad(Temp)
heat=dot(j,j)/ro
EQUATIONS
Volt: Div(j)=0
Temp: Div(q)=heat
BOUNDARIES
region 1 {Prostokat}
 lambda=1 ro=1
 start (0,0)
 natural(Temp)=k*(Temp-Ta)natural(Volt)=0
 line to (3,0)
 natural(Temp)=0 value(Volt)=Va
 line to (3,1)
 natural(Temp)=k*(Temp-Ta)natural(Volt)=0
 line to (0,1)
 natural(Temp)=0 value(Volt)=Vg
 line to close
region 2 {Kolo}
 lambda=10 ro=10
 start(1.5,0.25)
 arc(center=1.5,0.5) angle=-360
PLOTS
contour(Volt) vector(j)
contour(Temp) vector(q)
END