* opt *
* set * 
h/del *
vec/del *
opt nstat 
set xmgl 2.5 
set pass 5 
set cshi .015 
set vfon -100 
set vsiz 0.3 
set xlab 1.2 
set ylab 0.8 
set xval 0.5 
set yval 0.3 
set cfon -30 
set tfon -30 
set asiz 0.2 
set tsiz 0.3 
set ygti 1.1 
opt nbox 
opt grid 
set ksiz 0.08 
*
*trace on
set plci 1
set lwid 1
set ltyp 1
set mtyp 8
set mscf .5
opt nsta
* AAA if you "set stat 111111" the ROUT file is cut away !!
*set stat 111111
set stat 1111
* AAA if you "opt liny" the ROUT file is cut away !!
* opt liny
*
vec/cre vinp(7)
*
zone 1 2
vec/inp vinp(1) 1
*
vec/cre veq(8) R -15. -133. -165. 85. -20. -145. -12. -65.
vec/cre veqs(8) R 53. 38. 40. 80. 60. 65. 52. 35.
vec/cre vhou(8) R 1.5 4.5 7.5 10.5 13.5 16.5 19.5 22.5
vec/cre vhous(8) R 8*0
null 0. 24. -250 250.
hplot/err vhou veq vhous veqs 8 24 0.28 '1'
set ltyp 4
line 0 0 24 0
set ltyp 1
vec/cre par(1) R
vec/cre step(1) R 
vec/cre pmin(1) R 
vec/cre pmax(1) R
vec/cre errpar(1) R
vec/fit vhou veq veqs p0 'S' 1 par step pmin pmax errpar
* result: -80.6+-16.8 Chi2=16.1 DOF:7, Chi2N=2.3
fun/plot x**2.5*exp(-x/2) 0 30
line 16.1 0 16.1 1.
*
h/cre/1d 100 ' ' 100 0. 50.
h/cre/1d 101 ' ' 300 -200 100.
call pval.for
*
opt stat
set stat 1111
h/plot 101
h/plot 100
*wait
opt liny
h/plot 100
line 16.1 0 16.1 1000.
opt nstat
opt logy
h/plot 100
opt liny
end:
*