Differential constraints for bounded recursive identification with multivariate splines