Hi all, I try to assess the parameters (K1,K2) of a model that describes the adsorption of a molecule onto on adsorbent.
equation: dq/dt = K1*C*(qm-q)-K2*q I know the value of 'qm' and I experimentally measure the variables 'q', 'C', and the time 't'. t C q 1 0 144.05047 0.0000000 2 565 99.71492 0.1105625 3 988 74.99426 0.1722100 4 1415 58.65572 0.2129545 5 1833 48.34586 0.2386649 6 2257 40.29413 0.2587440 7 2675 32.92470 0.2771216 8 3105 29.57162 0.2854834 9 3552 28.01424 0.2893672 10 3986 25.62167 0.2953337 11 4415 23.62612 0.3003101 12 4841 21.95523 0.3044769 13 5264 21.08464 0.3066480 14 5698 19.68040 0.3101498 15 6509 18.31788 0.3135476 16 6950 17.65868 0.3151915 17 7403 17.00206 0.3168290 18 8130 16.38856 0.3183589 19 9001 15.58544 0.3203617 20 9928 15.27882 0.3211263 21 11899 14.46415 0.3231579 22 16354 13.91779 0.3245204 23 18926 13.82630 0.3247485 24 21602 13.66776 0.3251439 25 24413 13.98560 0.3243513 26 27056 13.87143 0.3246360 27 29844 13.64881 0.3251912 It's a differential equation, thus I had a look on the command 'ode' from the deSolve package. I'm early stuck on the use of the function 'ode' cause I don't get how to define the function 'func' required by 'ode' Any help would be appreciated. Regards/Cordialement ------------- Benoit Boulinguiez Ph.D student Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences Chimiques de Rennes" Avenue du Général Leclerc CS 50837 35708 Rennes CEDEX 7 Tel 33 (0)2 23 23 80 83 Fax 33 (0)2 23 23 81 20 http://www.ensc-rennes.fr/ <http://www.ensc-rennes.fr/> ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.