In[]:=
ClearAllend=6*model=
3
10
c-Kmeno3Kp+Kmeno3Kp+K3Kt-K3Kt-KpKt+KpKt
-Kpx-
K3Ktx
Kmeno3+Kt
Kpx
Kpx+
K3Ktx
Kmeno3+Kt
K3Ktx
Kmeno3+Kt
Kpx+
K3Ktx
Kmeno3+Kt
Kpx
Kpx+
K3Ktx
Kmeno3+Kt
Kmeno3Kp-K3Kt+KpKt
Out[]=
ClearAll
Out[]=
6000
Out[]=
c-Kmeno3Kp+Kmeno3Kp+K3Kt-K3Kt-KpKt+KpKt
-Kpx-
K3Ktx
Kmeno3+Kt
Kpx
Kpx+
K3Ktx
Kmeno3+Kt
K3Ktx
Kmeno3+Kt
Kpx+
K3Ktx
Kmeno3+Kt
Kpx
Kpx+
K3Ktx
Kmeno3+Kt
Kmeno3Kp-K3Kt+KpKt
In[]:=
data=ReadList["pori.dat",{Number,Number}]nlf=NonlinearModelFit[data,model,{{c,},{K3,},{Kmeno3,},{Kt,},{Kp,}},x]Show[ListPlot[data],Plot[nlf[x],{x,0,end},PlotStyleRed],FrameTrue]Normal[nlf]nlf["BestFitParameters"]
-5
10
-2
10
-3
10
-3
10
-3
10
Out[]=
{{0,0.},{30,2.0476×},{60,7.87794×},{90,1.68441×},{120,2.84323×},{150,4.22221×},{180,5.76087×},{210,7.43183×},{240,9.20596×},{270,1.10618×},{300,1.29784×},{330,1.49385×},{360,1.69257×},{390,1.8929×},{420,2.09375×},{450,2.29384×},{480,2.49238×},{510,2.6892×},{540,2.88367×},{570,3.07522×},{600,3.2638×},{630,3.44912×},{660,3.63087×},{690,3.80894×},{720,3.98323×},{750,4.15357×},{780,4.3199×},{810,4.48225×},{840,4.64057×},{870,4.79482×},{900,4.94507×},{930,5.0914×},{960,5.23383×},{990,5.3724×},{1020,5.50722×},{1050,5.63833×},{1080,5.76581×},{1110,5.88972×},{1140,6.01015×},{1170,6.12717×},{1200,6.24087×},{1230,6.3514×},{1260,6.45879×},{1290,6.56305×},{1320,6.66426×},{1350,6.76262×},{1380,6.85812×},{1410,6.95078×},{1440,7.04075×},{1470,7.12815×},{1500,7.21296×},{1530,7.29522×},{1560,7.37513×},{1590,7.45272×},{1620,7.528×},{1650,7.60101×},{1680,7.67193×},{1710,7.74078×},{1740,7.80756×},{1770,7.87234×},{1800,7.93527×},{1830,7.99635×},{1860,8.05558×},{1890,8.11299×},{1920,8.16861×},{1950,8.22252×},{1980,8.27478×},{2010,8.32544×},{2040,8.37457×},{2070,8.42222×},{2100,8.46846×},{2130,8.51334×},{2160,8.55693×},{2190,8.59929×},{2220,8.64047×},{2250,8.68054×},{2280,8.71947×},{2310,8.75723×},{2340,8.79385×},{2370,8.82937×},{2400,8.86381×},{2430,8.89722×},{2460,8.92964×},{2490,8.96109×},{2520,8.99162×},{2550,9.02126×},{2580,9.05005×},{2610,9.07801×},{2640,9.1052×},{2670,9.13163×},{2700,9.15732×},{2730,9.1823×},{2760,9.2066×},{2790,9.23018×},{2820,9.25305×},{2850,9.27523×},{2880,9.29674×},{2910,9.31758×},{2940,9.33778×},{2970,9.35737×},{3000,9.37637×},{3030,9.39482×},{3060,9.41274×},{3090,9.43014×},{3120,9.44702×},{3150,9.46339×},{3180,9.47929×},{3210,9.49471×},{3240,9.50966×},{3270,9.52417×},{3300,9.53826×},{3330,9.55192×},{3360,9.56518×},{3390,9.57804×},{3420,9.59053×},{3450,9.60264×},{3480,9.6144×},{3510,9.62581×},{3540,9.63687×},{3570,9.64762×},{3600,9.65804×},{3630,9.66816×},{3660,9.67799×},{3690,9.68757×},{3720,9.69689×},{3750,9.70597×},{3780,9.71479×},{3810,9.72335×},{3840,9.73166×},{3870,9.73972×},{3900,9.74752×},{3930,9.75507×},{3960,9.76237×},{3990,9.76941×},{4020,9.7762×},{4050,9.7828×},{4080,9.78923×},{4110,9.79548×},{4140,9.80157×},{4170,9.80749×},{4200,9.81323×},{4230,9.81881×},{4260,9.82421×},{4290,9.82945×},{4320,9.83451×},{4350,9.83941×},{4380,9.84413×},{4410,9.84869×},{4440,9.85312×},{4470,9.85744×},{4500,9.86165×},{4530,9.86575×},{4560,9.86973×},{4590,9.87359×},{4620,9.87735×},{4650,9.88099×},{4680,9.88452×},{4710,9.88793×},{4740,9.89123×},{4770,9.89442×},{4800,9.8975×},{4830,9.9005×},{4860,9.90342×},{4890,9.90626×},{4920,9.90903×},{4950,9.91172×},{4980,9.91434×},{5010,9.91687×},{5040,9.91933×},{5070,9.92172×},{5100,9.92402×},{5130,9.92625×},{5160,9.92841×},{5190,9.93049×},{5220,9.93252×},{5250,9.9345×},{5280,9.93642×},{5310,9.9383×},{5340,9.94012×},{5370,9.94189×},{5400,9.9436×},{5430,9.94527×},{5460,9.94688×},{5490,9.94844×},{5520,9.94995×},{5550,9.9514×},{5580,9.95281×},{5610,9.95419×},{5640,9.95553×},{5670,9.95683×},{5700,9.9581×},{5730,9.95934×},{5760,9.96053×},{5790,9.9617×},{5820,9.96283×},{5850,9.96392×},{5880,9.96497×},{5910,9.966×},{5940,9.96698×},{5970,9.96794×},{6000,9.96889×}}
-8
10
-8
10
-7
10
-7
10
-7
10
-7
10
-7
10
-7
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
-6
10
Out[]=
FittedModel
Out[]=
Out[]=
-1.05542(-2.36424×+0.0000118393-9.47504×)
-0.006005x
-6
10
0.000999561x
0.00500544x
-6
10
0.006005x
Out[]=
{c0.0000100002,K30.01134,Kmeno30.00132126,Kt0.00104403,Kp0.000999561}
In[]:=
fit=Plot[%%,{x,0,end}]
Out[]=