make delay state varying

[AugustHao]

current delay ecdf input can be time-varying but not state varying, need to do this in a computationally efficient way, will interact with the convolve and get_convolution matrix functions