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