demo/dummy_model_for_testing.R

In randall-romero/CompEconR: An R version of Miranda & Fackler's CompEcon toolbox for MATLAB

# dummy model for testing the library

require("CompEconR")

params <- list(r=0.02)


myfunc <- function(flag,s,x,i,j,e,params){
  R <- 1+params$r

  ns <- nrow(s)
  ds <- ncol(s)
  dx <- if(hasArg(x)) ncol(x)

  out <- switch(
    flag,
    f = list(
      f   =  log(s-x/R),
      fx  =  1/(x-R*s),
      fxx = array(-1/(x-R*s)^2,c(ns,dx,dx))
    ),
    g = list(
      g   = x,
      gx  = array(1,c(ns,ds,dx)),
      gxx = array(0,c(ns,ds,dx,dx))
    ),
    b = list(lower=matrix(0,ns),upper=R*s)
  )
  return(out)
}


# mymodel <- new("dynamic.model",
#                horizon=Inf,
#                func=myfunc,
#                discount=0.95,
#                params = params,
#                X=matrix(seq(0,8,by=0.1)))


mymodel <- new("dynamic.model",
               horizon=Inf,
               func=myfunc,
               discount=0.95,
               params = params
               )



mybasis <- fundefn(type='cheb',n=7,a=0.5,b=2)
#mybasis2 <- fundefn(type='cheb',n=c(4,5),a=c(-1,0),b=c(1,3))

myoptions <- new("dpsolve.options",algorithm="newton")
#myoptions <- new("dpsolve.options",algorithm="newton")


#myfunc(flag='b',s=mybasis@nodes,i=1,j=1,params=params)
solution <- dpsolve(model=mymodel,basis=as.ref(mybasis),options=myoptions)




#basis.getPhi(as.ref(mybasis),order=2)
#print(basis.getPhi(as.ref(mybasis),order=-2))

#sprintf('%10s %2d\n','hola',6)