vectorising all forloops (setting dP, enforcing cutoff)
deleted dynplot