HIGHER PRIORITY
-------------

Documentation
. document fsp
. proof read

momentDynamics():
. accept more general requestes for mu, e.g., include E[x] and E[x^2] 
  and all other monomials needed to close, up to order m
. in computation of Mu.sym and barMu.sym, the prod() operation should be done
  inside MuPAD to make things faster
. computation of .A and .B should eb sped up

readNet & Doc 
. "finishing without ';' termination does not yield error 
  and line is simply ignored, e.g.
  > derivatives:
  >	d/dt L += Y
  > EOF
. allow for keywords that reflect a more general use:
  e.g., 'species' can also be 'states', 'reactions' also 'jumps'
. allow for any order of sections (variable used in a switch statement 
  indicates with which section is being parsed, and keywords change
  this variable)

momentEquilibrium()
. is Abhi's LNA just lowDispersion?

support Monte Carlo simulations of nonzero derivatives and non-additive jumps.

Examples:
. Make all examples compatible with new syntax & momentEquilibrium()
.   gene/posfeed_nonlinear.*
.   gene/promoter_linear.*
.   gene/repressor_linear.*
.   bimodal/bistab2.*
.   bimodal/competence.*
.   bimodal/competence_simplest.*

.   PowerPoint/Biology/Banff07-stochasticity/matlab/promoter_linear.*
.   PowerPoint/Biology/Banff07-stochasticity/matlab/stochasticfocusing.*

.   research/biochemical/klavins/klavins.*
.   research/biochemical/division/division.*
.   research/biochemical/gene/matlab/

LOWER PRIORITY
-------------

Document
. momentDynamics()
. Document fullSSA(), 
. Document centerMoments

momentDynamics():
. Permit noise terms
. Permit normalized states (# molecules / typical # molecules)

.lnaClosure()
. speed up computation, perhaps using Jacobian?

readNet():
. Accept chemical reactions in the syntax of the .net file

readSBML():
. Understand volumes
. Understand units

Other
. Implement function that does (semi-automatic) expantion of 
  Michaelis-Menden kinetics into elementary reactions
