Table 1

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

Table 1: Basic procedure for Gibbs sampling algorithm.