TITLE	'Transient heat flow around an Insulating blob'
SELECT
	Bold = On
VARIABLES
	U				{ the temperature }
DEFINITIONS
	k = 1		{ default conductivity }
	c = 1		{ default heat capacity }
	r = 0.5		{ blob radius }
EQUATIONS
	c*dt(U) - dxx(k*U) - dyy(k*U) = 0
BOUNDARIES
	REGION 1
		START (-1,-1)
			VALUE(U) = 0			LINE TO (1,-1)
			NATURAL(U) = 0		LINE TO (1,1)
			VALUE(U) = 100		LINE TO (-1,1)
			NATURAL(U) = 0		LINE TO CLOSE
	REGION 2	{ the embedded blob }
		k = 0.001
		c = 0.1
		START (r,0)
			ARC(CENTER=0,0) ANGLE=360 TO CLOSE
TIME 0 TO 15
PLOTS
	FOR t = 1 BY 1 TO 15
	CONTOUR(U)
HISTORIES
	HISTORY(U) AT (0,r/2) (0,r) (0,3*r/2) FIXED RANGE (0, 100)
END