ample = 1.05295;
sigma = 0.509317;
model = a*x^c
data = Table[{x,mample*(1/(sigma*Sqrt[2*Pi]))*Exp[-(1/2)*(((x)/sigma)^(2))]}, {x,0, 2.1, 0.05}];
nlm = NonlinearModelFit[data, model, {a, c}, x]
Plot[nlm[x], {x, 0, 2.1}, Epilog -> Point[m], PlotRange -> All]