I am trying to fit my data to an inverse function from solving a DGL:
extf = 15
data = Import["http://pastebin.com/raw.php?i=mTQUNBua", "Table"];
ListPlot[data]
sol = DSolve[{a*L'[t] == b*L[t]^(-3/2) - c, L[0] == L0}, L, t]
Block[{b = 40000, c = 1.4 + extf, L0 = 10},
fit=FindFit[data, L[t] /. sol, {a}, t]]
Plot[fit,{t,0,100}]
Obviously this takes a huge amount of time. But after a while weird error messages start to show up. Is this in principle the right approach?
FindFit::nrlnum: "The function value {-47.5437+0.\ I,-32.8991+0.\ I,-22.3144+0.\ I,-13.9991+0.\ I,-7.09261+0.\ I,-1.16642+0.\ I,4.03917 +0.\ I,8.67588 +0.\ I,<<36>>,66.0943 +0.\ I,66.5301 +0.\ I,66.9401 +0.\ I,67.3269 +0.\ I,67.6895 +0.\ I,68.0316 +0.\ I,<<9950>>} is not a list of real numbers with dimensions {10000} at {a} = {1.}. "
$\endgroup$