Alpha C. Chiang. Elements Of Dynamic Optimization. 1992: Respostas do Capítulo 7 em R
Exemplo 1, página 171.
maxV=∫T0−(1−u2)1/2dtSujeito a:˙y=uy(0)=Ay(T) livreA, e T dados
O primeiro passo para a resolução do problema é instalar as bibliotecas necessárias e liberá-las para o uso. De início, considere usar a biblioteca deSolve, ideal para a resolução de equações diferenciais.
lapply(list("deSolve"), function(x){
if(x %in% installed.packages() == F){
install.