Using Mata function inside `nlcom`?

Question

I'm interested in a function that does not admit a closed-form expression, so I approximate it using a Mata function.

How can I call this Mata function from nlcom?

I perfectly understand that I could calculate numerical derivatives etc. on my own, but feels a bit unnecessary given that everything is already in nlcom.