It would be nice to allow sampling from the prior

Implementation should be straight-forward as a member function of Kalman