My problem is that given a dataset, I want to program fitting a gamma distribution on this data by estimating the two parameters(shape and the scale parameters) using Maximum Likelihood Estimation. I checked wikipedia article on MLP estimation of Gamma Distribution but it is not obvious there how we can have two equations where we can have the MLP estimate of the shape and the scale parameters and more importantly how to program them?
The answer is on the Wikipedia page. There is not a closed form solution so you have to use an iterative method like the one they have provided. 

