# Weak inhibition

Left: the time-dependent probability density of a system
describing *weak inhibition*. Right: two trajectories sampled
from this model. Note how the system essentially defines a lower
dimensional manifold: a very narrow "bumerang" is thus
formed. This means that the sampling in the right picture only
slowly reveals the full behaviour of the system.

The density in the left picture has been obtained by solving an
*aggregated* form of the system using a very high order
spectral method.

(Plot info: each axis tick is 500 molecules and the contour levels
have been chosen by Matlab)

#### References

*S. Engblom: Spectral Approximation of Solutions to the Chemical
Master Equation*, in *J. Comput. Appl. Math.*
229(1):208--221,
2009: (doi)
*S. Engblom: Galerkin Spectral Method applied to the
Chemical Master Equation*, in *Commun. Comput. Phys.*
5(5):871--896,
2009: (abstract), (pdf).

Stefan Engblom
Last modified: Thu May 10 11:04:19 MEST 2007