Variable delay with shift index bigger than zero

New solver, works for time dependent single delay with shift index up to
three