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)
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)