Torkel Loman

There are methods (like in SteadyStateProblem?) that are dependent on an initial condition one give, right? In these cases conversation law constraints should not be necessary (since these are implied...

Not sure exactly what you are looking for, but here are some ideas: ```julia f(du,u,p,t) = du[1] = 5 + 18*u[1]^4/(u[1]^4+11.9^4) - u[1] function g(du,u,p,t) du[1,1] = .1*sqrt(5 + 18*u[1]^4/(u[1]^4+11.9^4))...

The first one is just a very simple bistable system (X activating its own production and getting degraded). I guess people come across it quite often but I don't think...

You mean that it should rather be ```julia function g(du,u,p,t) ... du[1,3] = -0.15*2.5*sqrt(u[1]) ... du[2,3] = 0.15*2.5*sqrt(u[1]) ... end ``` I might be wrong but I think that is...

On the topic of spatial things, I've always imagined that the natural way to go is to have one `ReactionSystem` in each compartment (In most cases probably the same `ReactionSystem`)....

Seems we've been burned on this a few times already now... Is there an automatic-ish way of doing this, or would one just copy the code that's in the docs...

I have started working on this, but there are some issues over at https://github.com/JuliaDocs/Documenter.jl/issues which probably need to be sorted out before it would be possible to make sensible doctests....

While I agree ``` rx = @reaction k, 2A --> B ``` is much better than ``` rx = @reaction k, A+A --> B ``` is there a strong reason...

I am not sure I understand fully, but: - What do you want to do with your model? If you want to simulate time development from an initial condition, using...

I have to admit I'm not fully sure exactly what you request, but then I'm not a chemist, which seems what this mainly targets?