If I have the following data (**I have updated the data**):

      data={{{60, 1852.94}, {65, 178.035}, {70, 7.97143}, {75, 48.9479}, {80, 
      133.561}, {85, 8.65079}, {87, 1.78915}}

How can I fit it to the following equation?

    eq = 1/((x*c1)/((Exp[B1/(x - T0)])^0.75)*(1 - Exp[((Tm - x)*c2)/Tm^2])) /. {T0 -> 259.246,B1 -> 2595.89, Tm -> 88.2 + 273.15}

I am trying the following:


      fun[x_] = NonlinearModelFit[data,eq, {c1, c2}, x] // Normal (*c1=732975, c2=2.65721*10^6*)

Which gives me an erroneous `c1` and `c2` of `c1=732975` & `c2=2.65721*10^6`. When I do it in excel it gives me a `c1` of -65514626.34 and a `c2` of 68.4 which fits the data relative well. Why mathematica is not fitting the data or finding the correct `c1` and `c2` values?