//Potovanje delca po pravokotniku [-a/2,a/2]x[-b/2,b/2] po zakonih Brownovega gibanja.
//Skripta vrne verjetnost, da delec, ki zacne potovanje v (0,0), zapusti pravokotnik skozi levo ali desno stranico.

//stranici pravokotnika
real a=10,b=10; 

//opis robov pravokotnika
border leva(t=-b/2,b/2){x=-a/2; y=-t;label=1;};
border spodnja(t=-a/2,a/2){x=t; y=-b/2;label=2;};
border desna(t=-b/2,b/2){x=a/2; y=t;label=3;};
border zgornja(t=-a/2,a/2){x=-t; y=b/2;label=4;};

real n=100;//stevilo delilnih tock na vsaki stranici
mesh Th = buildmesh (leva(n)+spodnja(n)+desna(n)+zgornja(n));

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

solve Poisson(u,v,solver=LU) = 
int2d(Th)(dx(u)*dx(v) + dy(u)*dy(v))
+on(1,3,u=1)+on(2,4,u=0);

plot(u,fill=1,dim=3,wait=1,nbiso=100);
cout << "Verjetnost, da gre delec skozi levi ali desni del roba, je " << u(0,0) << endl;