[LONGITUDINAL]
input = {C, k1, k2, lambda0, N0}

EQUATION:
t0   = 0
N1_0 = 0
N2_0 = 0
N3_0 = 0
Np_0 = N0
Nt = Np+N1+N2+N3
ddt_N1 = k2*C*Np - k1*N1
ddt_N2 = k1*N1 - k1*N2
ddt_N3 = k1*N2 - k1*N3
ddt_Np = lambda0*Np - k2*C*Np

shinyUI(fluidPage(
  tabsetPanel(
    tabPanel("Plot",
             fluidRow(
               column(2,
                      br(),
                      sliderInput("lconc", label = "log10(C):", min = -2, max = 2, value = 0, step = 0.2),
                      sliderInput("k1", label = "k1:", min = 0.01, max = 0.2, value = 0.07, step = 0.01),
                      sliderInput("k2", label = "k2:", min = 0.01, max = 0.2, value = 0.07, step = 0.01),
                      sliderInput("lambda0", label = "lambda0:", min = 0.005, max = 0.1, value = 0.03, step = 0.005),
                      sliderInput("N0", label = "N0:", min = 500, max = 5000, value = 2000, step = 250)
               ),
               column(10,
                      plotOutput("plot",  height="370px"))
             )),
    tabPanel("Mlxtran", pre(includeText("delbene.txt"))),
    tabPanel("ui.R", pre(includeText("ui.R"))),
    tabPanel("server.R", pre(includeText("server.R")))
  )
))

# source("../../initMlxR.R")
library(mlxR)

shinyServer(function(input, output) {

  r <- reactive({  
    param.value <-c(10^input$lconc, input$k1, input$k2, input$lambda0, input$N0)   
    f  <- list(name='Nt',time=seq(0,70))
    p   <- list(name=c('C','k1','k2','lambda0','N0'), value=param.value)
    
    res <- simulx( model     = 'delbene.txt',
                   parameter = p,
                   output    = f)
    return(res)
  })  
  output$plot <- renderPlot({
    r=r()
    pl=ggplotmlx(data=r$Nt, aes(x=time, y=Nt)) + geom_line(size=0.75)
    print(pl)
  })
})