Different results when optimizing hyperparameter for a Gaussian process regression

Question

i\'m studying gaussian process regression, and i\'m trying to use the built-in functions from scikit-learn, and also trying to impement a custom function for doing so.

Th