Simulation of a multineuron interaction model

We simulate the dynamics of a multineuron model (RS model) with an energy function given by the product between the squared distances in phase space between the state of the net and the stored patterns. We obtain the relative frequency $f(m_0)$ that an arbitrary pattern is retrieved from an initial overlap $m_0$ and estimate the size of the basins of attraction for different activities $a$. Two limit cases are taken into account : when patterns and antipatterns are stored (p.a.s.) and when only the patterns are stored (o.p.s.). For the $a=0.5$, p.a.s. nets a limit for the load parameter was not found, but for the other cases ($a\neq 0.5$ or o.p.s. configuration) the relative size of the basins of attraction may become too small