def pvaj_to_pppp(p0, v0, a0, j0, delta_t):
    """
    Convert position, velocity, acceleration, and jerk at one timestep to positions for
    the next 4 timesteps.
    """
    a_arr = [a0, a0 + j0 * delta_t]
    v_arr = [v0]
    for a in a_arr:
        v_arr.append(v_arr[-1] + delta_t * a)
    p_arr = [p0]
    for v in v_arr:
        p_arr.append(p_arr[-1] + delta_t * v)
    return p_arr