Difference between revisions of "Lotka Volterra fishing problem (Muscod)"

From mintOC
Jump to: navigation, search
m
m (Moved C to Muscod)
Line 1: Line 1:
The differential equations for the [[Lotka Volterra fishing problem]] in [[:Category:C | C code]] read as follows
+
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:C]]
+
[[Category:Muscod]]

Revision as of 10: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);