Discrete Kalman filter

The intended way is to provide the distribution of the initial state, you do this by providing the d0 argument here, for example

d0 = MvNormal(x0, R0)