diff --git a/core/model.jl b/core/model.jl
index 8aa7b2ee4ce1c35da59868397a5faf64dcef0cfa..6867206fd27004e059233dc1937cf67ea7e91d07 100644
--- a/core/model.jl
+++ b/core/model.jl
@@ -47,9 +47,10 @@ end
 
 function simulate(m::ContinuousTimeModel)
     # trajectory fields
-    full_values = zeros(0, m.d)
-    times = zeros(0)
-    transitions = Vector{Union{String,Nothing}}(undef,0)
+    full_values = zeros(1, m.d)
+    full_values[1,:] = m.x0
+    times = Float64[m.t0]
+    transitions = Union{String,Nothing}[nothing]
     # values at time n
     n = 0
     xn = m.x0