By downloading this code you agree to the
Source Code Use License (PDF). |
{ version = 1.04; (* of fdr.p 1995 Jan 14 *)}
(* begin module describe.fdr *)
(*
name
fdr: computes best values of sigma/D for f_D(r) fit of real data.
synopsis
fdr(histog: in, fdrp:in, outfile: out, output: out)
files
histog: output file from genhis.
fdrp: parameter file for fdr.
sigma,sigmastep (reals): initial value and incremental increase in sigma.
dim,dimstep (reals): initial value and incremental increase in D.
temperature (real): Kelvin temperature of system.
outfile: writes best value of sigma/D and error in f_D(r) fit.
output: sends messages to the user.
description
D-1 2 2
r exp(-r /2sigma )
f (r) = -----------------------
D D D/2-1
Gamma(D/2) sigma 2
We compute the values of D, sigma that minimze the least-squares error
between f_D(r) and the actual data read in from histog.
examples
See fdr.eg (not yet written )
documentation
See Appendix 3 to TDS's CCMM paper.
see also
genhis.p
author
V. Jejjala (vishnu@ncifcrf.gov)
bugs
Unknown but existent (existential bugs, as it were... :-) ). If you find a
bug, please forward your discovery to the author.
disclaimer
The author of this code is in no way responsible for any havoc that this
program may cause to your computer system, your files, or your reputation.
Use this program and the data it generates at your own risk.
technical notes
Use large data sets in histog when possible.
*)
(* end module describe.fdr *)
{This manual page was created by makman 1.45}
{created by htmlink 1.62}