MaxIter =10,000; 
Initialize the flux v(0) with random value generator; 
for k =1:MaxIter 
v(k) = v(k 1); 
for l=1:n 
Calculate μ_{i} and Σ_{i} according to equation (10); 
Generate the sample θ from the Gaussian distribution
N(μ_{i},Σ_{i})
within the constraint Φ(∙); 
Update v_{i}(k) with the sample θ, v_{i}(k)=θ ; 
end 
end 
