Could anyone help with the code required for the McGinley Dynamic indicator?

MD = MD1 + (Price – MD1) / (k * N * (Price / MD1) ^ 4)

MD1 is the previous MD value calculated
k is a constant equal to 0.6
N is the moving average period
Price is the closing price

I had a look in indicators.c for some inspiration but trying to add the above is a bit beyond me at present.