Ok, this is what I did then:

Using:

(changed D to d since D means derivative in Mathematica)

I wrote the following Mathematica code to try and fit the data to your expression. Unfortunately, Mathematica returns a negative value for d which quickly causes overflow with the sum (even with smaller upper limits on the sum overflow occurs):

Code:

vals={{0,11.58},{30,72.61},{60,97.89},
{90,106.05},{120,109.63},{150,117.02},
{180,119.57},{210,125.99},{240,129.56},
{270,132.85}};
lp = ListPlot[vals, PlotRange -> {{0, 200}, {0, 300}}]
coef = FindFit[vals, M[t], {d, k}, t]
p1 = Plot[f[t] /. coef, {t, 0, 3}]

Mathematica returns:

{d -> -1.84547*10^8, k -> 122.325}

Oh yeah, if the equation is still not correct, click on it. A small window will appear with the LaTex code I wrote to generate the math code. You can understand it. \frac means fraction in french brackets, power is ^, so forth. with the [tex] brackets starting and ending it. Try and cut and paste it into a new post and edit the code as you wish even if it comes out messy. It's a start and you should try posting in Latex so everything is neat, well defined and explicit so people will follow it easily.