Hi there,

Is there a way in which i can solve the following optimal control problem numerically with Maple?

dH/dt=λ-µ H-(1-u1)β H V+δ I,

dI/dt=(1-u1)β H V-σ I,

dV/dt=(1-u2)k I-γV,

dλ1/dt=-1+λ1µ+βV(1-u1)(λ1-λ2),

dλ2/dt=λ1δ+λ2σ-λ3(1-u2)k,

dλ3/dt=β H(λ1-λ2)(1-u1)+λ3γ.

where

u1=βHV(λ1-λ2)/A1,

u2=-λ3kI/A2,

σ = α + δ,

and

H(0)=1.7*10^8, I(0)=0, V(0)=400.

λ=5*10^5, µ=0.003, β=4*10^(-10), α=0.043, δ= 0.2 k=6.24, γ=0.65.

A1=900, A2=1000.

Answers and advice are very appreciated.

Thank you all for reading.

Aylin