Hi Rik,
Numerical derivatives are inherently noisy, and subject to spikes due to time step differences, jumps in the solution, etc. You may want to filter the signal being differentiated. In Simulinik, a straightforward way to do a first-order filter is to use a transfer function block with 1/(tau*s + 1) where tau is the filtering time constant.
Martin
