Mean-Field Variational Inference for Gradient Matching with Gaussian Processes