Hello,

We are investigating Gap Junctions channels in cells (during cardiac arrhythmia also propagation of excitation in neural networks).

I am trying to implement our model in Julia, but first version ran about 10x slower than C++. Is this normal? because https://aaronang.github.io/2018/stl-benchmark-comparison-cpp-vs-julia/ and https://www.quora.com/Is-it-possible-for-Julia-to-perform-as-fast-as-C++-What-would-be-necessary-for-Julia-to-be-a-suitable-alternative-to-C++ states that Julia is just little slower than C++.

Model is based on Markov chain.

Our paper about model: https://www.sciencedirect.com/science/article/pii/S0006349516001600?via%3Dihub

to run - â€śrunMc36sm.jlâ€ť https://pastebin.com/LX5jAFV1

main model - â€śmc36sm_ss.jlâ€ť https://pastebin.com/QyaaW3Ti

commands to run:

```
using Revise, Plots;
plotlyjs();
includet("runMc36sm.jl"); includet("mc36sm_ss.jl"); includet("mc36sm.jl");
(vj, gj) = runGj();
```

also c++ version:

â€śmc36sm.cppâ€ť https://pastebin.com/cSgPUTjC

Kestutis