//Potovanje delca po krogu s srediscem S(0,0) in radijem r  po zakonih Brownovega gibanja.
//Skripta vrne verjetnost, da delec, ki zacne potovanje v (0,0), zapusti krog skozi del
// loka robne kroznice, ki je dolocen s srediscnim kotom fi.

real r=1.; //radij kroznice
real fi=pi/2.;//krozni lok, ki pridada kotu fi

//rob kroga je kroznica
border lok1(t=-fi/2,fi/2){x=r*cos(t); y=r*sin(t);label=1;};
border lok2(t=fi/2,2*pi-fi/2){x=r*cos(t); y=r*sin(t);label=2;};

real n1=100;//stevilo delilnih tock na loku
real n2=200;//stevilo delilnih tock na komplementu loka

mesh Th = buildmesh (lok1(n1)+lok2(n2));

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,u=1)+on(2,u=0);

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