c for fitting by x** function @ ANAPAW
REAL FUNCTION FUNC(X,Y) !for a 2-Dim histogram
COMMON/PAWPAR/PAR(2)
FUNC=-Y+PAR(1)*X+PAR(2)
END
こういうテキストファイルを作り、名前をfunc.fとする。
anapaw(paw)上で2次元の絵(histID=100)を出したら、
ANAPAW> His/fit 100 func.f ! 2 par
と打つと、以下のようにフィッティング結果が出る。
*** Unknown vector
PAR
1
**********************************************
* *
* Function minimization by SUBROUTINE HFITH *
* Variable-metric method *
* ID = 105 CHOPT = *
* *
**********************************************
Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01
FCN= 58617.65 FROM MIGRAD STATUS=CONVERGED 20 CALLS 21 TOTAL
EDM= 0.42E-08 STRATEGY=1 ERROR MATRIX UNCERTAINTY= 0.3%
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 P1 -0.23138E-01 0.74368E-03 -0.96424E-07 -0.13032
2 P2 139.58 0.67839E-01 0.10123E-05 -0.52769E-03
CHISQUARE = 0.4187E+02 NPFIT = 1402
0 件のコメント:
コメントを投稿