Rev 3482 Rev 3483
Line 16... Line 16...
16 h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1) 16 h4(x)=(T0/tau)*((p0/f4(x))**(R*tau)-1)
17 h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1) 17 h5(x)=(T0/tau)*((p0/f5(x))**(R*tau)-1)
18   18  
19 #plot "mereni_tlaku_16122013.log" using ($10):($17*100000) with linespoints title "P1", f1(x), "" using ($11):($17*100000) axes x1y1 with linespoints title "P2", f2(x), "" using ($12):($17*100000) axes x1y1 with linespoints title "P3", f3(x), "" using ($14):($17*100000) axes x1y1 with linespoints title "P4", f4(x), "" using ($15):($17*100000) axes x1y1 with linespoints title "P5", f5(x) 19 #plot "mereni_tlaku_16122013.log" using ($10):($17*100000) with linespoints title "P1", f1(x), "" using ($11):($17*100000) axes x1y1 with linespoints title "P2", f2(x), "" using ($12):($17*100000) axes x1y1 with linespoints title "P3", f3(x), "" using ($14):($17*100000) axes x1y1 with linespoints title "P4", f4(x), "" using ($15):($17*100000) axes x1y1 with linespoints title "P5", f5(x)
20   20  
-   21 a=1
-   22 b=1
21   23 c=1
22   24  
23 f1(x)=a*x+b*x**2+c 25 f1(x)=a*x+b*x**2+c
24 fit f1(x) "mereni_tlaku_16122013.log" using ($10/1000):(h1($10)) via a, b, c 26 fit f1(x) "mereni_tlaku_16122013.log" using ($10/1000):(h1($10)) via a, b, c
25 #fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2 27 #fit f2(x) "mereni_tlaku_16122013.log" using ($11):($17*100000) via k2, q2
26 #fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3 28 #fit f3(x) "mereni_tlaku_16122013.log" using ($12):($17*100000) via k3, q3