demo 1 	
	#  
	# USAGE:  in  sm   (supermongo)
	#  : macro read gauss_test.sm       <<= loads macro into sm
	#  : noisedemo  1000                <<= runs macro with 1000 pts
	#
	define N $1
	set x=1,$N
	# dummy call to get ws dimensioned right
	get_hist x v ws -.5 .5 .02
	set ws=ws*0
	erase
	window 1 1 1 1

	set y=x*0
	limits 0 $N -.7 0.7
	box
	ptype 1 1
	# generate a lot of points to uniformly cover -0.5 to 0.5
	# -and- at the same time do averages of 10 samples => pseudo-gaussian
	do i=1,10 {
		set t=random($N)-0.5
		set y=y+t
		# keep plotting the uniform points
		points x t
		get_hist t v w -.5 .5 .02
		set ws=ws+w
	}
	echo  * Uniformly Distributed between -0.5 and 0.5 *
	toplabel  * Uniformly Distributed between -0.5 and 0.5 *
	# use call to "cursor" in order to pause
	cursor
	# prepare and plot the pseudo gaussian points
	set y=y/10
	lweight 3
	ctype red
	connect x y
	ctype default
	lweight 1
	define ave (sum(y)/$N)
	define rms (sqrt(sum(y*y)/$N))
	define err ($rms/sqrt($N))
	echo rms=sigma= $rms
	echo average= $ave +- $err
	echo
	echo * red: averages of 10 random points *

	print noise_file { y }

 	cursor

	erase
	ctype default
	limits x y
	box
	connect x y
	set s=abs(x-600)<50?.06:0
	ctype red
	cursor
	lweight 3
	connect x s
	lweight 1
	set ynew=y+s
	ctype 0
	cursor
        erase
	box
	limits x ynew
	erase
	box
	connect x ynew
	stats ynew ave sig kur
	echo ...
	echo stats ynew ave sig kur
	echo ave=$ave    sigma=$sig
	echo ...
	echo vecminmax ynew mn mx
	vecminmax ynew mn mx
	echo min=$mn   max=$mx
	echo
	cursor

	smooth ynew ysmo 20
	echo ******************* smooth ynew ysmo 20
	ctype red
	lweight 3
	connect x ysmo
	ctype 0
	lweight 1

	cursor
	erase
	box
	connect x ysmo
	stats ysmo ave sig kur
	echo ...
	echo stats ysmo ave sig kur
	echo ave=$ave    sigma=$sig
	echo ...
	echo vecminmax ysmo mn mx
	vecminmax ysmo mn mx
	echo min=$mn   max=$mx    max/sig=$($mx/$sig)
	echo
	cursor
	ctype red
	connect x s
	ctype 0
	cursor	

 	limits x ynew
	erase
	box
	connect x ynew

	smooth ynew ysmo 50
	echo ****************** smooth ynew ysmo 50
	ctype red
	lweight 3
	connect x ysmo
	ctype 0
	lweight 1

	cursor
	erase
	box
	connect x ysmo
	echo set region=ysmo if(x<500&&x>50)
	set region=ysmo if(x<500&&x>50)
	stats region ave sig kur
	echo ...
	echo stats region ave sig kur
	echo ave=$ave    sigma=$sig
	echo ...
	echo vecminmax ysmo mn mx
	vecminmax ysmo mn mx
	echo min=$mn   max=$mx    max/sig=$($mx/$sig)
	echo
	cursor
	ctype red
	connect x s
	ctype 0
	cursor	




 	limits x ynew
	erase
	box
	connect x ynew
	cursor

	smooth ynew ysmo 100
	echo ********************* smooth ynew ysmo 100
	ctype red
	lweight 3
	connect x ysmo
	ctype 0
	lweight 1

	cursor
	erase
	box
	connect x ysmo
	echo set region=ysmo if(x<500&&x>50)
	set region=ysmo if(x<500&&x>50)
	stats region ave sig kur
	echo ...
	echo stats region ave sig kur
	echo ave=$ave    sigma=$sig
	echo ...
	echo vecminmax ysmo mn mx
	vecminmax ysmo mn mx
	echo min=$mn   max=$mx    max/sig=$($mx/$sig)
	echo
	cursor
	ctype red
	connect x s
 	ctype 0
	
        cursor
		limits x ynew
	erase
	box
	connect x ynew
	cursor

	smooth ynew ysmo 200
	echo ***************** smooth ynew ysmo 200
	ctype red
	lweight 3
	connect x ysmo
	ctype 0
	lweight 1

	cursor
	erase
	box
	connect x ysmo
	echo set region=ysmo if(x<500&&x>50)
	set region=ysmo if(x<500&&x>50)
	stats region ave sig kur
	echo ...
	echo stats region ave sig kur
	echo ave=$ave    sigma=$sig
	echo ...
	echo vecminmax ysmo mn mx
	vecminmax ysmo mn mx
	echo min=$mn   max=$mx    max/sig=$($mx/$sig)
	echo
	cursor
	ctype red
	connect x s
 	ctype 0