search for: nonstiff

..., I'm struggling to solve this ODE using R, vdpol <- function (h, v, t) ( list(c ( -0.1*v/(pi*(2*10*h-h^2)), (v = (-0.1*v/(pi*(2*10*h-h^2))^2) + 2*9.81*h)) )) library(deSolve) yini <- (c(h = 20, v=0)) nonstiff <- ode(y = yini, func = vdpol, times= seq(0, 30, by = 0.01), parms = 1) It says that The number of derivatives returned by func() (4) must equal the length of the initial conditions vector (2) I'm not sure how when i only h...

...mponent in error vector 0 8 The length of the real work array actually required: 1732 9 The length of the integer work array actually required: 110 14 The number of Jacobian evaluations and LU decompositions so far: 4834 15 The method indicator for the last succesful step, 1=adams (nonstiff), 2= bdf (stiff): 2 16 The current method indicator to be attempted on the next step, 1=adams (nonstiff), 2= bdf (stiff): 2 -------------------- RSTATE values -------------------- 1 The step size in t last used (successfully): 2.189586e-06 2 The step size to be attempted on the ne...