@everywhere using MarkovProcesses
using Plots

@everywhere load_model("ER")
@everywhere load_automaton("automaton_F")

A_F_R2 = create_automaton_F(ER, 50.0, 75.0, 0.05, 0.075, :P)
sync_ER = A_F_R2*ER 
pm_sync_ER = ParametricModel(sync_ER, (:k3, Uniform(0.0, 100.0)))

@timev r = automaton_abc(pm_sync_ER)
@show r.nbr_sim

histogram(r.mat_p_end', bins = :sqrt, weights = r.weights, normalize = :density, xlims = (0.0, 100.0))
savefig("R2_hist.svg")