Ok, I have the following equations: c'[t] = a[t]*b[t]*k1 - c[t]*k2, a'[t] = (-a[t])*b[t]*k1 + c[t]*k2, b'[t] = (-a[t])*b[t]*k1 + c[t]*k2, c[0] =0, a[0] =10, b[0] =10. I have data for concentration of "c" at all time intervals from 0-16. Here is the data: {0, 0}, {1, 4.98112}, {2, 6.61986}, {3, 7.42423}, {4,7.89529}, {5, 8.20004}, {6, 8.41008}, {7, 8.56123}, {8,8.67342}, {9, 8.7586}, {10, 8.82442}, {11, 8.87595}, {12, 8.91671}, {13, 8.94923}, {14, 8.97535}, {15, 8.99643}, {16, 9.01351}. I need Maple to Solve the ODE's and optimize the values for k1 and k2 to fit the data usins a least-square method or something along the lins of that. Can someone send me the code for that?