I have created a dataset from the following:

set obs 10
gen goals=_n-1
gen alpha=1.4
gen beta=1

I want to use a gamma distribution function to estimate the likelihood for each observation in goals given the alpha and beta parameters defined. I believe i need to do this in mata but have never used this programming language and am at a loss of to create the function in mata and then access that function to generate a new variable

Giulio DiDiodato