Lotka Experimental Design (VPLAN)

VPLAN

Differential equations:

c     RHS of the differential equations
 
      subroutine ffcn( t, x, f, p, q, rwh, iwh, iflag )
        implicit none
 
        real*8 x(*), f(*), p(*), q(*), rwh(*), t
        integer*4 iwh(*), iflag
 
        real*8 p1,p3,p5,p6, myu
 
c	fixed parameters
	p1 = 1.0
	p3 = 1.0
	p5 = 0.4
	p6 = 0.2
 
c       DISCRETIZE1( myu, rwh, iwh )
 
 
        f(1) = p1*x(1)        - p(1)*x(1)*x(2) - p5*myu*x(1)            
        f(2) = (-1.0)*p3*x(2) + p(2)*x(1)*x(2) - p6*myu*x(2)
 
      end

First Measurement function:

c     Messfunktion
 
      subroutine mess3( t, x, h, p, q, rwh, iwh, iflag )
        implicit none
 
        real*8 t, x(*), h, p(*), q(*), rwh(*)
        integer*4 iwh(*), iflag
 
        h = x(1)       
 
        iflag = 0
 
      end

Second Measurement function:

c     Messfunktion
 
      subroutine mess4( t, x, h, p, q, rwh, iwh, iflag )
        implicit none
 
        real*8 t, x(*), h, p(*), q(*), rwh(*)
        integer*4 iwh(*), iflag
 
        h = x(2)       
 
        iflag = 0
 
      end

Standard deviation of first measurement function

c     Standardabweichung der Messfunktion
 
      subroutine sigma3( t, x, s, p, q, rwh, iwh, iflag )
        implicit none
 
        real*8 rwh(*)
        integer*4 iwh(*), iflag
 
        real*8 t, x(*), p(*), q(*)
        real*8 s
        real*8 h
 
        s = 1.0d+0
 
        iflag = 0
 
      end

Standard deviation of second measurement function:

c     Standardabweichung der Messfunktion
 
      subroutine sigma4( t, x, s, p, q, rwh, iwh, iflag )
        implicit none
 
        real*8 rwh(*)
        integer*4 iwh(*), iflag
 
        real*8 t, x(*), p(*), q(*)
        real*8 s(*)
 
        s(1) = 1.0
 
        iflag = 0
 
      end