Lotka Volterra fishing problem (Muscod): Difference between revisions
mNo edit summary |
m Moved C to Muscod |
||
| Line 1: | Line 1: | ||
The differential equations for the [[Lotka Volterra fishing problem]] in [[:Category: | The differential equations for the [[Lotka Volterra fishing problem]] in [[:Category:Muscod | Muscod code]] read as follows | ||
<source lang="cpp"> | <source lang="cpp"> | ||
| Line 13: | Line 13: | ||
</source> | </source> | ||
[[Category: | [[Category:Muscod]] | ||
Revision as of 09:18, 28 January 2016
The differential equations for the Lotka Volterra fishing problem in Muscod code read as follows
/* steady state with u == 0 */
double ref0 = 1, ref1 = 1;
/* Biomass of prey */
rhs[0] = xd[0] - xd[0]*xd[1] - p[0]*u[0]*xd[0];
/* Biomass of predator */
rhs[1] = - xd[1] + xd[0]*xd[1] - p[1]*u[0]*xd[1];
/* Deviation from reference trajectory */
rhs[2] = (xd[0]-ref0)*(xd[0]-ref0) + (xd[1]-ref1)*(xd[1]-ref1);