//Prevajanje toplote na pravokotni plosci, z dvema izoliranima robovoma
//in dvema robovoma, kjer temperaturo vzdrzujemo na predpisani vrednosti.

//konstante in funkcije
real k=0.1;//toplotna prevodnost
func q0=1;//izvor toplote znotraj plosce
func T0=0;//temperatura na zgornji in desni stranici kvadrata
func Tn=0;//normalni odvod na levi in spodnji stranici kvadrata

//konstante geometrije obmocja
real x0=0,x1=1;
real y0=0,y1=1;
int m=10,n=10;
mesh Th=square(m,n,[x0+(x1-x0)*x,y0+(y1-y0)*y]);//mreza mxn na kvadratu [x0,x1]x[y0,y1]


fespace Vh(Th,P1);
Vh T,v;

solve Poisson(T,v,solver=LU) = 
int2d(Th)(dx(T)*dx(v) + dy(T)*dy(v))
-int2d(Th)(q0*v/k)
//-int1d(Th,4)(Tn*v)
//-int1d(Th,1)(Tn*v)
+on(1,4,T=T0);


plot(T,fill=1,dim=3,wait=0,nbiso=100,value=1);