; this one test linefitting commands
disp1d 1 disp2d 0 cdisp2d 0
read fit_test
evaln 33 142 addbase %% max
zoom 1 636 807
minimax (10*$noise) % peak showpeaks
report ($npk1d == 21) "peak"

zoom 1 643 687 
showpeaks
iter 10
linefit lorentz
set a = (abs($chi2-36.4)<1)
peak showpeaks
linefit gauss
report ($a &  (abs($chi2-10.2)<1)) "Linefit"

zoom 0

dim 1 chsize 50
one tm 50 50  put tab                ; create a dump TAB[]
one mult 1.5 specw 1000 em 100       ; create a dump data-set
iter 10   miniter 10
addnoise .02 123

set p1 := 1
set p2 := 0.1

FITGENE '$p1*exp(-$x/$p2)' 2
; ../macro/showexp '$p1*exp(-$x/$p2)'
set a = (abs($p1-1.75)<0.05)
set a = ($a & (($p2-0.1273)<0.05))
set a = ($a & ($chi2 < 1e-3))
report $a "fitgene"

set p1 := -1
set p2 := 1
MINIMIZE 'sin($p1)*sin($p2)' 2
set a = (abs($p1+2*atan(1))<0.05)
set a = ($a & (abs($p2-2*atan(1))<0.05))
set a = ($a & (abs($chi2+1) < 1e-3))
report $a "minimize"

unset p1  unset p2

