Skip to main content
33 votes
Accepted

Turing patterns

I developed a reaction-diffusion-advection model of pattern formation in semi-arid vegetation (tiger bush) 20 years ago, which shows a type of Turing instability. Plants ($n$) consume water ($w$) and ...
Chris K's user avatar
  • 20.4k
22 votes

Turing patterns

I did some work with the Brusselator some time ago. This is the reaction-diffusion equations which generate Turing patterns. There are some things you need to know: (1) The non-linear PDEs have ...
Dominic's user avatar
  • 2,914
12 votes
Accepted

Solving ODEs vs VoronoiMesh

The following example makes use of my package, IGraph/M. Needs["IGraphM`"] First we generate some random points to use as Voronoi centres. We apply a few steps ...
Szabolcs's user avatar
  • 236k
9 votes
Accepted

How to use Manipulate for protein sequence alignment?

I'm using as examples the mouse and human sequences for Rhodopsin (the protein in retinal rod cells that detects photons). ...
MelaGo's user avatar
  • 9,111
6 votes

Performing maximum likelihood estimation for SIR Model

ParametricNDSolve and NonlinearModelFit work for this cases. ...
cvgmt's user avatar
  • 84.1k
6 votes

How do I plot this system of ODEs using Manipulate to change the b and k variables?

...
MarcoB's user avatar
  • 67.7k
6 votes
Accepted

Plotting a 4-dimensional graph with the fourth dimension depicted using a color scheme

...
Bob Hanlon's user avatar
  • 163k
4 votes
Accepted

Avoiding Concave Polygons

The logic might need some tweaking depending on the order of the points, but for these examples, it works. Essentially, replace the computation of the two new endpoints ...
Lukas Lang's user avatar
  • 34.4k
4 votes
Accepted

Modeling Infection Spread among Different Age Groups with Contact Matrix

I really enjoy your epidemic series :-) It is always better to decouple the weedy details of a specific model from the key question. Key question should be purified to bare minimum. Perhaps helpful, --...
Vitaliy Kaurov's user avatar
3 votes
Accepted

Plotting Dengue Infection Prevalence Sensitivity to Biting Rate Variations

...
Alex Trounev's user avatar
  • 48.8k
3 votes

Implementing Izhikevich's spiking neuron model

We can use input data from the page. Tonic Spiking model can be implemented as follows ...
Alex Trounev's user avatar
  • 48.8k
3 votes
Accepted

Plotting several curves of the same variable for different values of a parameter

ParametricNDSolve is the tool for the job, combined with Table: ...
MarcoB's user avatar
  • 67.7k
3 votes
Accepted

How to plot contour plots of R0 against two parameters?

...
Bob Hanlon's user avatar
  • 163k
3 votes
Accepted

Insert, remove or shift two random time points of data

...
kglr's user avatar
  • 401k
3 votes

Bioinformatics: RNA-seq and Gene Set Analysis Pipelines

Gene List Pre-Processing ABOUT: In this example we have two gene lists, both from mice. One gene list contains columns of features, the other has condition replicates. GOAL: The goal is to Convert ...
SumNeuron's user avatar
  • 5,442
2 votes

Difficulties in compiling the steady state of a biological model

Comment Maple gives multiple solutions, ...
zhk's user avatar
  • 12k
2 votes
Accepted

Birth Death Process with Delay

I figure it out how to simulate birth-death process with one species. Black curve is DDE (delayed DE). ...
OkkesDulgerci's user avatar
2 votes
Accepted

SIR Parameter Estimation in Mathematica

It looks like we try to estimate parameter $\beta$ with known $\gamma, \mu$ using one point only reportedData = n*0.08. In this case we can compute $\beta$ as ...
Alex Trounev's user avatar
  • 48.8k
2 votes
Accepted

Performing maximum likelihood estimation for SIR Model

...
Athanasios Paraskevopoulos's user avatar
1 vote
Accepted

Writing a sum-of-squares function or SIR model fits

Just changing the last little bit: ...
ydd's user avatar
  • 6,944
1 vote

Non Linear System of Diferential Equations with periodic function

Your main problem was not to define all paramters. A good way to look whether you got all variables is Variables[Level[eqs, -1]] Two posibilities of solution: a) ...
Akku14's user avatar
  • 17.4k
1 vote

Insert two random time points to data

Perhaps you can make use of StepFunction from my answer to question How can the behavior of InterpolationOrder->0 be controlled?: ...
Carl Woll's user avatar
  • 132k
1 vote
Accepted

Insert two random time points to data

...
kglr's user avatar
  • 401k
1 vote

Difficulties in compiling the steady state of a biological model

Here's how I would do it: ...
Chris K's user avatar
  • 20.4k

Only top scored, non community-wiki answers of a minimum length are eligible