version 1.1.1.1, 2000/01/09 17:01:00 |
version 1.1.1.3, 2003/09/15 07:09:29 |
Line 11 set title 'data for first fit demo' |
|
Line 11 set title 'data for first fit demo' |
|
plot 'lcdemo.dat' |
plot 'lcdemo.dat' |
set xlabel "Temperature T [deg Cels.]" |
set xlabel "Temperature T [deg Cels.]" |
set ylabel "Density [g/cm3]" |
set ylabel "Density [g/cm3]" |
|
set key below |
|
|
print "now fitting a straight line to the data :-)" |
print "now fitting a straight line to the data :-)" |
print "only as a demo without physical meaning" |
print "only as a demo without physical meaning" |
Line 53 print "It's time now to try a more realistic model fun |
|
Line 54 print "It's time now to try a more realistic model fun |
|
load 'density.fnc' |
load 'density.fnc' |
show functions |
show functions |
print "density(x) is a function which shall fit the whole temperature" |
print "density(x) is a function which shall fit the whole temperature" |
print "range using a ?: expression. It contains 6 model parameters which |
print "range using a ?: expression. It contains 6 model parameters which" |
print "will all be varied. Now take the start parameters out of the" |
print "will all be varied. Now take the start parameters out of the" |
pause -1 "file 'start.par' and plot the function (-> return)" |
pause -1 "file 'start.par' and plot the function (-> return)" |
load 'start.par' |
load 'start.par' |
Line 106 print "\n\nNotice, however, that this would converge m |
|
Line 107 print "\n\nNotice, however, that this would converge m |
|
print "fitted in a more appropriate co-ordinate system:" |
print "fitted in a more appropriate co-ordinate system:" |
print "fit r 'hemisphr.dat' using 0:($1*$1+$2*$2+$3*$3) via r" |
print "fit r 'hemisphr.dat' using 0:($1*$1+$2*$2+$3*$3) via r" |
print "where we are fitting f(x)=r to the radii calculated as the data" |
print "where we are fitting f(x)=r to the radii calculated as the data" |
print "is read from the file. No x value is required in this case. |
print "is read from the file. No x value is required in this case." |
pause -1 "(This is left as an excercise for the user). (-> return)" |
pause -1 "(This is left as an excercise for the user). (-> return)" |
FIT_MAXITER=0 # no limit : we cannot delete the variable once set |
FIT_MAXITER=0 # no limit : we cannot delete the variable once set |
|
|
Line 195 set title 'initial parameters' |
|
Line 196 set title 'initial parameters' |
|
plot 'moli3.dat' w e, R(x) |
plot 'moli3.dat' w e, R(x) |
pause -1 "now start fitting... (-> return)" |
pause -1 "now start fitting... (-> return)" |
fit R(x) 'moli3.dat' u 1:2:3 via eta, tc |
fit R(x) 'moli3.dat' u 1:2:3 via eta, tc |
pause -1 "now look at the result (-> return)" |
pause -1 "Hit return to continue" |
set title 'fitted parameters' |
set title 'fitted parameters' |
replot |
replot |
|
|
Line 213 print "Remember that this file will always be appende |
|
Line 214 print "Remember that this file will always be appende |
|
print "from time to time!" |
print "from time to time!" |
print "" |
print "" |
pause -1 "Done with fitting demo (-> return)" |
pause -1 "Done with fitting demo (-> return)" |
|
reset |