Dear all, I am successfully playing with sol = solve(ensembleprob,SRIW1(),EnsembleThreads(),trajectories=ϕ,saveat=δt); which generates a nice set of values. Among those values sol[1,:,:]; is a position, sol[2,:,:]; a velocity. unfortunately, I need an acceleration… how can I nicely proceed ? t…

sadly ended with a five-point stencil scheme : for i in 3:(lmnt-2) kt_a[i] = (-kt[i+2]+8*kt[i+1]-8*kt[i-1]+kt[i-2])/(12*δt) end; :roll_eyes:

Derivative for a set of ... points

rafael.guerra March 16, 2021, 2:36pm 3

@Frederic_Gilt, you may want to take a look at this discourse thread.

Topic		Replies	Views
New to julia: Convolution and juliadiff Numerics question	19	6283	August 21, 2017
Simple numerical differentiation New to Julia question , forwarddiff , finitediff	5	4595	December 23, 2020
How to find second and third derivative using ForwardDiff New to Julia	2	2202	February 19, 2019
How to take the derivative of a function? New to Julia	3	4574	March 14, 2022
Generating finite-difference stencils Numerics differentiation	19	3451	July 23, 2024