version 1.1.1.2, 2000/01/22 14:16:06 |
version 1.1.1.3, 2003/09/15 07:09:29 |
Line 54 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 107 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 196 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 |
|
|