So you want the code for Kriging using gstat. A simple Google search shows plenty of tutorials. You just need to apply the principles showed on these tutorials to your data.

I would advise you to start with Kriging and then move to co-Kriging.

Thank you for your reply dear Nikolaos Tziokas

As you said, there are many texts, but some of them are complicated for a beginner

I would like to take advice and guidance from colleagues who are more familiar with the use of the package

Best

Check the example in the link below:

https://github.com/r-spatial/gstat/blob/main/demo/examples.R

Thank you dear Muhammad Waqas Ahmed

I'll see it