Again in the book "chemical kinetics with maple and mathcad", chapter 4 introduces different ways for rate constants determination.
Take figure 4.8 as an example (corresponding attached file is "Fig_4_08_PineneIsomerization_InvTask_MC14" ), numerical integration method is used searching rate constants:
but i think that solving ODEs is more direct and general than integration method, so I solve the problem with ODEsolve method:
Both methods looks good enough, but differential method results are closer to that of DYANFIT than for integral method.
But differential method seems much more initial-condition sensitive than integral method, e.g. if guess values of k1--k4 are all set to 1, the calculation will be very time-consuming. and I dont know how to improve this situation.
Does anybody give any suggestions on "lowering/ reducing" the initial - condition sensitivity of differential method?